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Abstract 



Inflationary cosmology is the leading explanation of the very early universe. Many 
different models of inflation have been constructed which fit current observational data. 
In this work theoretical and numerical methods for constraining the parameter space 
of a wide class of such models are described. 

First, string-theoretic models with large non-Gaussian signatures are investigated. 
An upper bound is placed on the amplitude of primordial gravitational waves produced 
by ultra-violet Dirac-Born-Infeld inflation. In all but the most finely tuned cases, this 
bound is incompatible with a lower bound derived for inflationary models which exhibit 
a red spectrum and detectable non-Gaussianity. 

By analysing general non-canonical actions, a class of models is found which can 
evade the upper bound when the phase speed of perturbations is small. The multi- 
coincident brane scenario with a finite number of branes is one such model. For models 
with a potentially observable gravitational wave spectrum the number of coincident 
branes is shown to take only small values. 

The second method of constraining inflationary models is the numerical calculation of 
second order perturbations for a general class of single field models. The Klein-Gordon 
equation at second order, written in terms of scalar field variations only, is numerically 
solved. The slow roll version of the second order source term is used and the method 
is shown to be extendable to the full equation. This procedure allows the evolution 
of second order perturbations in general and the calculation of the non-Gaussianity 
parameter in cases where there is no analytical solution available. 
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1. Introduction 



In the past cosmology was a speculative science. The scarcity of observational data 
meant that many conflicting theories for the evolution of the universe were entertained, 
with nothing but personal opinion to differentiate between them. The explosion in 
the quantity and quality of observational data in recent years has led to a much more 
competitive marketplace of ideas about the physical beginning of the universe. 

The Big Bang scenario has emerged as a cohesive framework for the evolution of the 
universe from very early times. The observation of the Cosmic Microwave Background 
(CMB) provided much supporting evidence for this scenario ["' "]. This relic radia- 
tion, emitted 300,000 years after the Big Bang, continues to be our primary source of 
information about the early universe. 

The inflationary scenario is an attempt to solve problems with the standard Big Bang 
picture and provide an origin for the fluctuations in energy that seeded the growth of 
structure in the universe [5, 78, 119, 190, 191]. These fluctuations link the classical 
scales of relativistic gravity with the quantum scales of Planck level physics. There are 
many possible realisations of inflation and there has been an explosion in the number 
of theoretical models which agree with current observational limits (for reviews see, for 
example, Refs. [1, 18, 111]). 

Ground and space-based observations have significantly challenged theoretical cos- 
mological models with a wealth of new data. The Wilkinson Microwave Anisotropy 
Probe (WMAP) mission [ ], in conjunction with supernova surveys and other evi- 
dence, have shown that the fluctuations in the temperature of the CMB are 10^ times 
smaller than the background value and that the magnitude of the fluctuations is roughly 
independent of the angular scales at which they are measured. This is in agreement 
with the predictions of inflationary models and has led to other scenarios being ruled 
out. An upper bound has been placed on the amplitude of gravitational wave pertur- 
bations and bounds have also been placed on the deviation of the fluctuations from a 
purely random Gaussian distribution. 
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Constraining the parameter space of inflationary models is an important step towards 
limiting the number of observationally viable models, and ultimately towards identifying 
one such model as the best candidate to describe the physics of the early universe. 

The goal of this thesis is to constrain inflationary models in two very different ways: 
by deriving analytic limits on their parameter spaces, and by demonstrating a numerical 
calculation which will allow the investigation of higher order perturbations. Both these 
methods have the potential to limit the parameter space of the models investigated and 
possibly to rule them out. 

In Chapter 2 the foundations are laid for these investigations. The geometry and 
physics of the Friedmann-Robertson- Walker universe are presented and inflationary cos- 
mology is introduced to alleviate problems with the standard Big Bang scenario. Slow 
roll conditions are then defined to ensure an adequate duration of inflation. Despite 
its elegance, this homogeneous cosmology does not provide an adequate description of 
our universe. To understand the inhomogeneities that are present in reality, first or- 
der cosmological perturbation theory is employed. Models with non-canonical actions 
can also be considered. The relationships between observable quantities and the model 
parameters are altered in this case, meaning these models could be distinguished from 
those with canonical actions. The departure of primordial perturbations from a Gaus- 
sian random distribution could also reveal significant information about the underlying 
physics at work. 

In Part I of this thesis, analytical bounds are placed on a class of non-canonical 
inflationary models. These models illustrate the dynamics of extended objects called 
branes in superstring theory and are considered to be some of the most promising 
candidates for achieving inflation using string theory. 

Chapter 3 outlines the Dirac-Born-Infeld (DBI) scenario in terms of the string the- 
oretic background and how it applies in four dimensions as a realisation of inflation. 
The six extra dimensions required by string theory play an integral role in this scenario. 
These are compactified into a complex manifold whose geometry allows extended regions 
called throats to exist. DBI inflation consists of a brane moving in one of the throats. 
The inflaton field is the radial distance of the brane from the tip of the throat. Trans- 
lating the higher-dimensional motion into four dimensions introduces a non-canonical 
term into the effective action. The real nature of the action then enforces an upper 
bound on the kinetic energy of the inflaton, allowing a sufficiently long period of in- 
flation. The total inflaton field variation is directly linked to the amplitude of tensor 
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modes which can be produced. 

In Chapter 4 the repercussions of this relationship between the change in the field 
value and the tensor mode amplitude are explored further. In the DBI scenario, Bau- 
mann & McAllister [19] placed a conservative upper bound on the total production of 
tensor modes during infiation, by assuming the brane does not propagate further than 
the length of the throat. By considering only the period of observable infiation, which 
takes place over a much smaller region of the throat, we have derived a new bound 
which is considerably stronger. In the generic case, the ratio of the amplitudes of the 
tensor and scalar perturbations must be less than 10~^. This is below even the most 
optimistic forecasts for the sensitivity of future observational experiments. 

If attention is limited to brane motion down the throat, another complementary 
bound on the tensor modes can be derived, which depends on the non-Gaussianity 
of the scalar modes produced during infiation. The DBI scenario is inherently non- 
Gaussian in nature, but, even assuming the largest levels allowed by observations, the 
tensor-scalar ratio must exceed 0.005. These two bounds are clearly incompatible in 
the generic case and only a very fine-tuned selection of model parameters allows the 
standard DBI scenario to survive. By taking a more phenomenological approach and 
allowing the other parameters to vary, conditions are found under which the bounds 
can be relaxed. 

A more general class of models which evade the upper bound are identified in Chap- 
ter 5. The DBI scenario is characterised by a simple algebraic relation, in which the 
sound speed of fiuctuations is inversely proportional to the contribution to the non- 
Gaussianity. By allowing the proportionality constant to vary, a new family of actions 
is derived for which the bound on the tensor-scalar ratio can be relaxed. 

Instead of considering a single brane moving in the throat, a more natural scenario 
might involve multiple branes. These could be created from the energy released by 
a brane/anti-brane annihilation and could move up the throat away from the tip. In 
Ref. [] Thomas & Ward described the case when these branes are coincident. When 
a large number of branes coincide, the resultant action is similar to the single brane 
action and is restricted by the bounds on the tensor-scalar ratio. For a small, finite 
number of branes, however, the action is non-Abelian in nature and is one of the 
family of "bound-relaxing" actions described above. Nevertheless, this model is still 
constrained by observations and, if a detectable tensor signal is required, only two or 
three coincident branes are allowed. This limit on the number of branes is strongly 
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dependent on the non-Gaussianity and a tightening of the observational bounds could 
rule out the possibility of an observable tensor signal from this model. 

In Part II, the focus of the thesis moves from analytical to numerical techniques. Sec- 
ond order cosmological perturbations are numerically calculated for single field canon- 
ical inflationary models. 

In Chapter 6, the system of equations for the numerical calculation is developed. In 
order to understand non-linear perturbative effects, it is necessary to examine models 
using perturbation theory beyond first order. The gauge transformation for second 
order perturbations is outlined and the effect on scalar quantities is considered in the 
uniform curvature gauge. In Ref. [ ' the Klein-Gordon equation for second order 
perturbations was written in terms of the field perturbations alone. This forms the 
basis of the numerical calculation once it is transformed into Fourier space. As the 
original equation involves terms quadratic in the first order perturbations, the Fourier 
transformed equation contains convolutions of these perturbations. As a first step 
towards demonstrating the calculation for the full equation, the slow roll version of the 
source term is considered in the second order equation. The second order perturbations 
can be linked to observable quantities including the curvature perturbation and the 
non-Gaussianity parameter. 

The Klein-Gordon equations are the central governing equations of the calculation 
described in Chapter 7. They must first be rewritten in a form more suitable for nu- 
merical work. This involves changing the time coordinate to the number of elapsed 
e-foldings and writing the convolution terms in spherical polar coordinates. Four dif- 
ferent potentials will be investigated, each of which has a single field which is slowly 
rolling. The parameters for these models are set by comparing the calculated power 
spectrum of first order scalar perturbations with the latest WMAP data. The initial 
conditions for the background field and perturbations must also be specified. The sec- 
ond order perturbations are initially set to zero, to highlight the creation of second 
order effects. As this is a novel procedure, a thorough description of the implementa- 
tion of the calculation is given. Where an analytic solution for the convolution terms 
is possible, this is compared with the calculated value. Numerical parameters are set 
by minimising the relative error in the calculation of one of the terms. 

The results of the numerical calculation are presented in Chapter 8. Three different 
ranges of the discretised momenta are considered and general results presented for the 
quadratic potential. As expected for a single field, slow roll model, the second order 
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perturbations are highly suppressed compared to the first order ones. The source term 
of the second order perturbation equation is similar in form to the power spectrum of the 
first order perturbations. It decreases rapidly until horizon crossing after which a more 
steady amplitude is maintained. The results for all four potentials are also compared. 
Differences are apparent in the behaviour of the models after horizon crossing. This 
calculation represents only the first step towards a full numerical integration of the 
second order Klein-Gordon equation. The next stages towards this goal are outlined. 
The second order equation for single field models without the slow roll assumption is 
written in the correct form for numerical use and the second order equations for the 
two field case are presented in vector form. 

In Chapter 9 the results of the thesis are discussed and some final conclusions are 
presented. 

Conventions 

Throughout this thesis units are chosen such that Mpl = (^^ixGY^I'^ = 2.4 x 10^^ GeV 
defines the reduced Planck mass and c = h = 1. 

An overdot (') is used for differentiation with respect to proper time t and a prime 
(') for differentiation with respect to conformal time r]. From Chapter 7 onwards, the 
dagger symbol (^) denotes differentiation with respect to the number of e-foldings Af. 
A subscripted comma denotes partial differentiation by the symbol it precedes, e.g. 

^''^ dip- 

The (++-1-) convention in the notation of Misner et al. [195] is used throughout. 



2. Inflationary Cosmology 



In this chapter the foundations of inflationary cosmology are described. In Section 2.1 
the physics of an isotropic and homogeneous universe is reviewed. The inflationary 
scenario is introduced in Section 2.2. First order cosmological perturbation theory is 
presented in Section 2.3 and inflationary models with non-canonical actions are de- 
scribed in Section 2.5. The current observational limits on inflationary models are 
outlined in Section 2.4 and departures from Gaussian statistics are parametrised in 
Section 2.6. 

2.1. The Friedmann- Robertson- Walker Universe 

The cosmological principle is central to the Friedmann- Robertson- Walker (FRW^) Uni- 
verse. According to this postulate, there is no privileged place in the universe and no 
privileged direction in which to make observations. These assertions are formalised by 
assuming that the universe is homogeneous and isotropic at every point. This clearly 
conflicts with the highly inhomogeneous nature of matter on planetary and solar system 
scales, but is assumed to hold as larger and larger scales are considered. Surveys of 
the observable universe indicate that this assumption is valid up to the largest scales 
observed [51, 207]. Historically, homogeneity and isotropy were assumed primarily for 
simplicity. Many alternative approaches can be taken. Violating these assumptions 
can be done, for example, by specifying a preferred direction or supposing that the 
universe is formed by a series of voids connected by filaments. Although many of these 
approaches have been disregarded due to lack of evidence, some are still allowed by 
observations [6, 8, 68, 82]. 

This section outlines the dynamics of the standard Big Bang scenario. By assuming 
homogeneity and isotropy, the equations of motion of a fluid-filled universe can be 

^Lemaitre is sometimes also included in this group to give FLRW. 
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derived. What follows here is a standard exposition of well-known physics and has 
been the subject of numerous reviews including Refs. [103, 114, 195]. 

By imposing both homogeneity and isotropy on a general 4-dimensional metric, the 
line element ds^ of the FRW universe with coordinates (t, r, 6', u) is obtained: 

ds^ = -dt^ + a\t) (^YZ^ + (^^' + sm'(^)dc^')) > (2-1) 

where K = +1,0 or —1 depending on whether the universe is closed, flat or open 
respectively. The time-like coordinate in the metric is t, known as proper time. The 
spatial part of the FRW metric is multiplied by the scale factor a{t). This characterises 
the size of space-like hypersurfaces at different times. In an expanding universe, a grows 
with increasing t and d > 0. The definition of the Hubble parameter, H, captures this 
expansion: 

H=-. (2.2) 

a 

The Einstein equations can be derived by the variational principle from the action S, 
where S = Seh + Sm- This is the sum of the Einstein-Hilbert (S'eh) and matter (S'm) 
actions which are defined as 



^— j d^a;^(/? + 2A,) , (2.3) 
Sm= I d\^\CM. (2.4) 



Here g is the determinant of the metric (yf^jy, R is the Ricci scalar, G is Newton's 
gravitational constant, Ac is a cosmological constant term and Cm is the sum of the 
Lagrangian densities for all the matter fields. Changing either the matter or gravity 
actions will affect the resultant physics. In this work we focus our attention only on 
the matter Lagrangian and will use the standard Einstein-Hilbert action throughout. 
We can now write down the Einstein equations for a general matter Lagrangian: 

Ri,u - ]^Rgt,u = SnGTf,^ + AcQ^^u , (2.5) 

where T^,y is the stress energy tensor obtained by the variation of the matter Lagrangian. 
In the definitions above we have included a cosmological constant term, A^,, for com- 
pleteness. In the early universe this term is sub-dominant and will be negligible until 
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much later [I ! ' ] . From now on we will disregard the contribution of such a term in the 
early universe. 

We concentrate now on the case of a universe filled with a perfect fluid. Suppose u'^ 
is the 4- velocity of this fluid with m^m^ = — 1. The stress-energy tensor of the fluid is 

Tt = iE + P)u^u, + P5l^ , (2.6) 

where E is the matter energy density and P is the isotropic pressure. The trace of T 
is given by 

T'; = -E + 3P. (2.7) 

The Einstein equations and the stress-energy tensor of the perfect fluid can now be 
used to derive the equations of motion of the fluid. From the metric in Eq. (2.1), the 
00 and ij components of the Ricci tensor can be found: 

Roo = —3- , (2.8) 
a 

Rij = -fij [2(2^ + ad + 2K] , (2.9) 

where is the time independent spatial part of the metric in Eq. (2.1). The Friedmann 
equations are then determined from the Einstein equations (2.5). The 00 equation gives 

while the trace of the Einstein equations gives the Raychaudhuri or acceleration equa- 
tion 

a AttG , , , , 

- = —{E + 3P). 2.11 

a 6 

By combining these two equations we can determine a continuity equation for the energy 
density: 

E + 3H{E + P) = 0. (2.12) 

The last three equations, (2.10), (2.11) and (2.12), will determine the evolution of the 
perfect fluid. Two important solutions of these equations are the radiation and matter 
dominated universes. In the standard Big Bang scenario the universe is dominated by 
radiation to a good approximation until matter becomes dominant at later times [103]. 
These different components change the rate of expansion of the universe. For relativistic 
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radiation Prad = -f'rad/S and integrating the continuity equation (2.12) gives Ej-m oc a~^. 
Matter conversely is taken to be dust-like with zero pressure and so -Ematter cl^^- The 
dependence of a on t can then be found from Eq. (2.10), giving a oc t^/^ and a oc 

^2/3 

for the radiation and matter eras respectively. 

Instead of using proper time as above we could bring the scale factor outside the 
whole metric and use conformal time rj defined by 

r dt , , 

V= -■ 2.13 
J a 

The metric written in conformal time is then 

ds^ = a\ii) (^-dr/2 + + (d^^ + sin\e)duj'')^ , (2.14) 

As all the coordinates in the line element are now scaled by a(?7), we have defined 
a coordinate grid which does not change as the universe expands. These "comoving" 
coordinates allow distances to be compared at different eras with ease. A comoving 
distance x can be translated into a physical distance d by 

d = ax. (2.15) 

The physical distance changes as the universe expands but the comoving distance will 
remain fixed. 

One particularly important distance is the maximum distance light could have prop- 
agated from some initial time ti to a later time t. From Eq. (2.16), this is simply the 
conformal time integrated from the initial time and is called the comoving or particle 
horizon. If the initial time is restricted to being at some finite time in the past, as in 
the Big Bang scenario, then the particle horizon will be finite. Two points which are 
further apart than this finite distance could never have been in causal contact. This is 
the origin of one of the major problems with the standard Big Bang scenario and will 
be discussed in the next section. Rewriting the comoving horizon as 

rda' 1 

shows that it is also the logarithmic integral of the comoving Hubble radius 1 / aH. This 
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distance is how far particles can travel in one "e-folding" , the time for a to expand by 
one exponential factor. The number of e-foldings between two measurements of the 
scale factor, and aj, is given by 



Particles that are separated by more than the Hubble radius cannot be in causal 
contact now. Particles separated by more than the comoving horizon, however, could 
never have been in causal contact. In addition to the Hubble parameter H, it will be 
useful to define the parameter "H = aH = a' /a. The comoving Hubble radius is then 



In this section we introduce the inflationary scenario. First we briefly describe how it 
solves two major problems with the standard Big Bang picture: the flatness problem 
and the horizon problem [ ]. We go on to describe canonical slow roll inflation, the 
generation of perturbations from quantum fluctuations and inflation from non-canonical 
actions. 

2.2.1. Problems with the Big Bang Scenario 

Although remarkably successful in describing the evolution of the universe from very 
early in its history, the standard Big Bang scenario suffers from a number of serious 
problems. Two of the main problems are described in this section. 

Flatness Problem 

The Friedmann equation (2.10) can be re- written as 



(2.17) 



1/n. 
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n{t) - 1 



K 




(2.18) 



(a/7)2 



where Q{t) = E{t) /Ecrit and the critical density E'crit = 3H'^/8ttG. If a > then Q 
approaches the critical value = 1 over time, whereas if a < it diverges from this 
value. The flat universe, = 0, is an unstable flxed point in the parameter space. 
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Current observations confirm that Q = 1 within about 2%, at a 95% confidence level 
[104]. During the radiation and matter dominated eras aH is decreasing with time, so 
that fl diverges away from 1. If the measured value is now very close to 1 then in the 
past it must have been even closer. The fine-tuning in the initial conditions required 
for this proximity to f2 = 1 is known as the flatness problem. 

Horizon Problem 

The particle horizon, also known as the comoving horizon, defines the maximum sep- 
aration between two points that have been in causal contact sometime in the past. 
During the radiation and matter eras, this comoving horizon increases monotonically 
and so length scales which are now entering the horizon would have been far outside 
it in the past. The CMB as observed by the WMAP satellite is extremely smooth on 
scales that would have been far outside the horizon at the time of last scattering [104]. 
These regions of space have very similar energies and yet according to the Big Bang 
scenario they could never have been in causal contact. 

2.2.2. Inflation and Canonical Slow Roll 

Inflation is a period of accelerated expansion in the size of the universe which took 
place just after the Big Bang [5, 78, 119, 190, 191]. During this expansion phase the 
comoving Hubble radius {aH)~^ decreases and the isotropic pressure of the universe is 
negative [18, 114]: 




a>0 E + 3P <0. 



(2.19) 




(2.20) 





a 



This parametrisation illustrates that inflation only occurs when Eh < ^- In this subsec- 
tion we describe briefly how inflation solves the problems outlined above and outline 
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the inflationary dynamics of single scalar field models. 

Both the horizon and flatness problems described above are statements about our 
reluctance to impose fine-tuned initial conditions. Inflation removes the need to fix 
these conditions at the start of the Big Bang. A period of decreasing Hubble radius 
before the radiation period could explain the homogeneity of temperatures in the CMB 
at large scales. Comoving scales that entered the horizon recently, such as those we 
observe in the CMB, would have been within the horizon previously. During this period, 
the energy density could reach an equilibrium value. Figure 2.1 shows how, by extending 
the era of inflation far enough into the past, any comoving length could previously have 
been inside the horizon. Observations require that inflation lasted at least long enough 
that all the scales we measure today were previously inside the horizon. 



log(fc-i) 




Large scale 



Small scale 



Inflation 



Radiation 



Matter 



e-foldings 



Figure 2.1.: Comoving scales that have recently entered the horizon would previously 
have been inside the horizon, if the inflationary period extended far enough 
into the past. 
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Now consider the time derivative of |f2 — 1| as defined in Eq. (2.18): 

If the universe is not flat to begin with, a period of inflation of sufficient duration will 
push it towards = 1, solving the flatness problem. Instead of an unstable point in 
the parameter space, f2 = 1 is an attractor during the inflationary phase. 

To solve the horizon and flatness problems the duration of the inflationary era must be 
sufficiently long. Approximately 50-70 e- foldings is considered standard [ ]. Inflation 
can last longer than this but only the last 50-70 e-foldings will be important for the 
length scales of our observable universe. 

During inflation the universe is filled by material exhibiting negative isotropic pres- 
sure. Therefore, whatever drives inflation cannot be matter or radiation in their usual 
forms. The simplest proposal is to flll the universe with a single scalar fleld (p. The 
canonical action for this field is 



(2.22) 



Cm = P{^,X)=X-V{^), (2.23) 

where X = —}^g^^d^(pd^ip denotes the kinetic energy of v?, V{ip) is the potential and 
P{ip,X) is called the kinetic function. In Section 2.5 we will consider other choices for 
P. 

The equation of motion for (f, for the canonical action P = X — V^, is 

AT/ 

¥3 + 3Hip - VV + = . (2.24) 

If we now restrict ourselves to considering the homogeneous part of the field, if = 
(p{t), the V'^ip term disappears and the functional derivative of V becomes a standard 
derivative V^^. With these choices we have the following relations for the matter energy- 
density and isotropic pressure: 

E=^^' + Vicp) , (2.25) 
P=l^'- V{^) . (2.26) 

Under these conditions the kinetic function P{f, X) can be identified as the isotropic 
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pressure. The dynamics of the field are governed by the potential V{(f). Infiation 
requires P < —E/3, so from Eqs. (2.25) and (2.26) infiation can also be thought of as 
a period when the potential energy dominates over the kinetic energy. 

Infiation needs to last long enough to solve the problems described above. A generic 
potential is not likely to satisfy these requirements without fine-tuning. One approach is 
to enforce conditions on the potential under which the infiationary period is necessarily 
long. We have seen that for infiation to occur the potential needs to dominate over 
the kinetic energy. From Eq. (2.25), this occurs in the limit P — ?■ —E or equivalently 
Eh ^- However, for this to remain the case for a sufficiently long period the second 
derivative of (f must be small. If we define another parameter 

- ^^^^ ^ ^oo7^ 
dma H(f ^Heh 

then taking ItjhI <^ 1 ensures that ip and Eh change slowly. This allows an infiationary 
phase of sufficient duration to occur. 

The approximations ^ 1 and |?7| <^ 1 are known as the slow roll conditions because 
they force the infiaton field if to roll down the potential V slowly. The parameters Eh 
and rjH are the slow roll parameters. Setting if to be small is equivalent to making 
the friction Hif in Eq. (2.24) dominant. With these approximations the equations of 
motion for a slowly rolling field become 

(2.28) 

^ ^-^Viip) . (2.29) 



2.3. Perturbations 

We considered a homogeneous scalar field in the analysis of Section 2.2. Such a field, 
however, will lead only to a homogeneous universe later. How does the myriad structure 
that we see around us form? From stars to galaxies to clusters, the gravitational 
force has concentrated energy density over the history of the universe, but some initial 
fluctuation must have been present to begin this process. One of the main achievements 
of infiation is to provide a physical origin for such initial fluctuations. In Chapter 6 
we will formally develop cosmological perturbation theory up to second order. In this 
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section we review first order perturbation theory and introduce the observable quantities 
important for inflation. 

Suppose that a full inhomogeneous scalar field is split into a homogeneous back- 
ground field (fo, as described above, and an inhomogeneous perturbation 6ip{ri,x^) for 
i = 1,2,3. For the following analysis to be applicable the perturbation must be much 
smaller than the background field value. From the amplitude of perturbations in the 
CMB this approximation can be seen to be valid [iOl]. For single field models no 
mixing of adiabatic and non-adiabatic modes occurs [203]. Therefore, throughout this 
thesis we will only consider adiabatic perturbations and ignore any isocurvature mode 
present. 

If we suppose that e is a small quantity then the split in can be written as [136] 

^{r],x') = Mv) + eS^{^,x'). (2.30) 

The perturbation 5(p{ri,x^) can be further expanded in powers of e. We will follow the 
custom of not explicitly writing e, instead relying on the order of the perturbation, 
denoted by a subscript, to keep track. If we expand in a Taylor series then up to second 
order (i.e. including terms up to e^) we have: 

^{t], x') = Mv) + SMV, x') + ^Sifi2{v, x') . (2.31) 

There is some freedom in how the split of the perturbations into different orders is made. 
We will suppose that the first order perturbation 6(pi contains only linear contributions 
and the higher order terms contain non-linear terms. 

Instead of working in coordinate space, we can also consider the perturbation in 
Fourier space using the definition 

Sy^iv, x') = j d^kS^ik') expithx') , (2.32) 

where k^ are the components of the comoving wavenumber vector k. The amplitude 
of this vector A; = |k| identifies whether a particular mode is inside or outside the 
comoving horizon. Wavemodes inside the comoving horizon are identified by A; > aH, 
while k < aH for those outside the horizon. 

We must also consider perturbations in the metric tensor (yf^j,. If the background 
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metric is the FRW one described in Section 2.1 then the metric can be written with 
perturbations, up to first order, as follows: 

fi-oo = -a^{l + 20i) , 

g,, = a" {5i, + 2Cuj) ■ (2.33) 

The 0-2 and i-j perturbations can be decomposed into scalar, vector and tensor parts 
[136]: 

Bii = Bi^i — Sii , 

Clij = -i^iSij + Ei^ij + Fi(ij) + ^hiij , (2.34) 

where = ^(-^lij + -^ij,i)- The vectors Su and Fu are divergence free and the 

tensor part huj is divergence free and traceless: 

Stk = 0, F,\, = 0; /if,, = 0, hU = 0. (2.35) 

In the previous equations is the lapse function, is the curvature perturbation, Bi 
and El are the scalar part of the shear, Su, and Fu are the vector parts of the shear, 
and hiij is the tensor perturbation describing gravitational waves. 

Splitting an inhomogeneous spacetime into background and perturbation is not a 
covariant operation. This leads to an ambiguity in the choice of coordinates which 
must be rectified by choosing a gauge. Gauge transformations relate physical results in 
one gauge to those in another. To choose a gauge one must specify how spacetime is 
foliated, i.e., a shcing, and how coordinates in one spatial hypersurface are related to 
those in another, i.e., a threading [136]. We will employ the uniform curvature gauge 
in which spatial hypersurfaces are fiat. This is also known as the fiat gauge. 

The gauge transformation vector at first order, can be split into scalar and vector 
parts 

er = (ai,/3i; + 7i), (2.36) 
where the vector part obeys '^i^ = 0. A scalar quantity such as the infiaton perturba- 
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tion will transform as [134, 136] 

6^pi = 6ipi + if'^ai , (2.37) 

where a tilde (~) denotes a transformed quantity. For the metric perturbations the 
transformations for the scalars are 

01 = 01 + -Hai + a[ , (2.38) 

^1 = ^1 - -Hai , (2.39) 

B,=B,-a, + P[, (2.40) 

^i=^i + /3i, (2.41) 

and for the vectors 

S^' = S,' - 7f , (2.42) 

F\' = F^' + Yi. (2.43) 

The tensor perturbation huj does not change under transformation at first order, but 
does at subsequent orders. The fiat gauge, which we will use, is the one in which spatial 
hypersurfaces are not perturbed by scalar or vector perturbations, so ipi = Ei = and 
Fii = 0. This is equivalent to the transformation 

«i = ^, fii = -Ei, l\ = -Fl. (2.44) 

A gauge invariant infiaton perturbation variable is the Sasaki-Mukhanov variable 
[145, 147, 170]: 

^^1 = 5^1 + ^'^^, (2.45) 

In the flat gauge this is just 5ipi. We will work in flat gauge from now on and so will 
drop the tildes on quantities in that gauge. 

Another very important gauge invariant quantity is the comoving curvature pertur- 
bation TZ. At first order in the flat gauge TZ is related to the infiaton perturbation by 
[136] 

7^ = ^5(^1 . (2.46) 
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We are interested in the power spectrum of the curvature perturbation as this is directly 
related to the temperature fluctuations that we can observe in the CMB. 

The action (2.23), including perturbations of and g^i, up to first order, is varied 
to get the equation of motion of 6ipi. In the flat gauge the equation can be rewritten 
in terms of the inflaton field values only by eliminating the metric perturbations using 
Eq. (2.45). In Fourier space and in terms of the conformal time rj, the closed form of 
the first order perturbation equation of motion is [136] 



n 



n 



Sifii{k') = 0, (2.47) 



where Vq is the background value of the potential V{ip). Substituting u = aSipi gives 
the Mukhanov equation [147] 



u"{k') + 



k' 



z 
z 



u{k') = 



(2.48) 



where z = aipQ/Ti. 



2.3.1. Quantum Perturbations 

So far we have considered classical perturbations. However, the generation of fluc- 
tuations is a quantum effect and we need to consider the perturbations as quantum 
operators in some vacuum. 

In Minkowski space the quantisation of u{k^) is straightforward. The perturbation 
modes can be written in terms of quantum operators as 

u{k') u{k') = w{k')a{k') + w*{-k')a\-k') . (2.49) 

The mode function w{k^) obeys the same equation of motion as u{k^): 

w{k')=0. (2.50) 

The operators and d are the usual creation and annihilation operators. They act on 
quantum states by adding or removing particles. The zero particle vacuum state, |0), 



w"{k') + 



k'- 
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is such that 

at|0) = |l), a|0)=0. (2.51) 
In Minkowski space these operators have the usual commutation relations 

[a(k), at(k')] = (27r)35(k - k') (2.52) 

and 

[at(k), at(k')] = [a(k), a(k')] = . (2.53) 
The w modes are normalised by the condition [144] 

w\k')w\k') - w*\k')w{k') = i . (2.54) 

In the expanding FRW background the choice of vacuum is not straightforward. 
Suppose one observer selects a zero particle state as the vacuum. Another observer 
accelerating with respect to the first will see particles being created in this "vacuum" 
state due to the Unruh effect [9(i, 19<S]. In selecting the vacuum we must choose one 
of the many equivalent options. To do this we consider the far past where ?7 — )■ — oo. 
The wavelengths of all the modes are then much smaller than the Hubble radius and 
curvature scale. The modes are therefore assumed to evolve in flat space. This suggests 
the Minkowski vacuum as the most natural vacuum state to select and this choice of 
vacuum at early times is known as the Bunch-Davies vacuum. In the limit 7] — —oo 
(or equivalently k/aH — > oo), the mode equation (2.50) becomes 

w"{k') + k^w{k') = 0, (2.55) 

which has the plane wave solution 

w{k') = ^e-'^"^ . (2.56) 



This is the initial condition for modes which are well inside the horizon. 

Now consider the de Sitter limit in which eh ^ and H is constant. We have 
z"/z = a" /a = 2/rf so the mode equation is [18] 



w"{k') + 



w{k') = Q. (2.57) 
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A full general solution for w is 

w{k') = A^— (l - ^] + b""— [l + ^] . (2.58) 

Taking the condition (2.54) along with the solution for subhorizon modes in Eq. (2.56) 
we find that A = 1 and 5 = 0. Thus the full solution in de Sitter space is [114] 




(2.59) 



Inflation in spacetimes that are close to de Sitter will contain perturbations with a 
spectrum defined by Eq. (2.59). The slow roll approximation is enough to ensure that 
inflation occurs in a quasi-de Sitter spacetime. However, the initial conditions for 
Fourier modes in Eq. (2.56) apply to non slow roll models so long as they are applied 
well before horizon crossing. 



2.3.2. Power Spectra and Spectral Indices 

The power spectrum of the inflaton perturbation 6ipi = u/a can now be deflned as 

(5v.i(kOM(k2)) = (27r)3<5(ki + k,)Pl{h) = {2n)'6ik, + ks)^^^^^ , (2.60) 

where (...) denotes the ensemble average. If taken over a large enough volume, the 
ensemble average and spatial average are equivalent [125]. The power spectrum 
depends only on the magnitude of the wavenumber vector, k = |k|, but has dimensions 
of k~^. A dimensionless power spectrum can be deflned as 

Vl=Al^—Plik). (2.61) 

In a similar way we can deflne the power spectrum of the comoving curvature per- 
turbation TZ = H6{pi/(po: 



(7^(kl)7^(k2)) = (27r)=^(5(ki + k2)Pl{k,) , 



(2.62) 
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and the dimensionless power spectrum 

Vl = ^]^^^2Pn{k). (2.63) 

A slow roll inflation model in a quasi-de Sitter spacetime will have the Fourier mode 
solution given in Eq. (2.59). After horizon crossing, when k <^ aif, this gives \w\'^ = 
1/(2A;V) so 



2 



for the scalar perturbation spectrum and 



Kik)=(^X(^] , (2.65) 



\ipoJ V27r 

for the comoving curvature perturbation spectrum. Models that are not slowly rolling 
usually require their more complicated mode equations to be numerically solved. 

We have discussed in depth the scalar perturbations but tensor perturbations can 
also be produced. The tensor perturbation hij has two polarisations, hg for s = +,x. 
The amplitude of each can be thought of as a separate scalar field. The analysis for 
each field is similar to that above with the substitution hg = 26ipi/Mpi^. After horizon 
crossing in a quasi-de Sitter space the spectrum for each polarisation is 



-(-)' 



'Plik) = ^{-] , (2.66) 



and the overall tensor perturbation spectrum is 

(2.67) 

The ratio of the tensor to curvature perturbations (tensor-scalar ratio) r is defined as 

r=:^, (2.68) 
where r is usually quoted at a particular k but could in principle depend on k. The 
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tensor-scalar ratio can also be written in terms of Eh- 

r = IGeh . (2.69) 

As ^ 1 for slow roll models of inflation the amplitude of tensor perturbations that 
these models produce is much smaller than the amplitude of curvature perturbations. 

If the curvature perturbation power spectrum, V^{k), is independent of wavenumber 
k, it is said to be scale invariant. The spectral index Ug is a measure of the deviation 
from scale invariance: 

where Ug = 1 denotes a scale invariant spectrum. The spectral index of the tensor 
power spectrum can be similarly defined: 

although this definition means that the spectrum is scale invariant if ut = 0. The 
spectral indices and indeed the spectra themselves are usually calculated at an arbitrary 
pivot scale. The WMAP results for V-^ and outlined in Section 2.4 are quoted at 
the scale k = 0.002Mpc"^ 

If there is a non-trivial dependence of or on k then higher order derivatives 
can be taken to give the running of the quantities. The runnings of the spectral indices 
are 

d In d In ut ,^ 

"^=dhr¥' "^ = ^hrF- (2-^'^ 

In the slow roll approximation Ug and Ut can be written in terms of the slow roll 
parameters en and rju, evaluated at k = aH using dln(aif) ~ Hdt: 

ng-l = -4eH + 2r]H, (2.73) 
nT = -2eH- (2.74) 

Combining Eq. (2.74) and Eq. (2.69) gives a powerful consistency condition for slow 
roll inflation: 

r = -8nT. (2.75) 
For the slow roll approximation to be valid for single field canonical inflation mod- 
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els, Eq. (2.75) must hold. Current observations are not accurate enough to test this 
condition but it is hoped that this will be possible in the future. 

2.4. Current Observations 

There have been rapid improvements in the quantity and quality of cosmological data 
sources in the last twenty years. From the launch of the COBE satellite in 1989 [25, 26], 
through the currently ongoing WMAP mission [104, 189], to the recent launch of the 
Planck satellite [158], space based observations have been at the forefront of the effort 
to collect data. Complementing these have been ground and balloon based missions 



including CBI [138, 182, 183], VSA [51], ACBAR [107, 108] and BOOMERANG [143, 



Major recent data releases have provided significant confirmation of the FRW model 
of the universe. The Hubble parameter today has been measured as Hq = 72 ± 
8 km/s/Mpc by the Hubble Key Project [66]. The WMAP 5- Year data release (WMAP5) 
[104] quotes their results combined with data from Baryon Acoustic Oscillations in 
galaxy distributions (BAO) [156] and supernova surveys (SN) by the Hubble Space 
Telescope and others [11, 162, 163, 205]. This combined data constrains the universe 
to within two percent of the flat Q = 1, K = case outlined in Section 2.1. 

The amplitude of the scalar curvature perturbations V-^ was first measured accurately 
by the COBE satellite [25, 26]. The WMAP5 normalisation is taken at a different scale 
to the COBE result, measuring 



where the pivot scale fcwMAP = 0.002Mpc~ ~ 5.25 x 10~^°MpL. The spectral index of 
scalar perturbations for models with tensor-scalar ratio r 7^ is given by the combined 
WMAP5+BA0+SN measurement as 



157, 166]. 



P^(A;wMAp) = 2.457 X 10-^ 



(2.76) 



n, = 0.968 ±0.015. 



(2.77) 



The detection of B-mode polarisation would provide definitive proof of the existence 
of primordial gravitational modes and much observational effort is being expended in the 
attempt to achieve such a detection [21, 46, 157, 176, 183, 199]. The observational bound 
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on r from WMAP5 using only the B-mode power spectrum is weak with r < 4.7 at the 
95% confidence level, when is fixed at the best fit value. Including other polarisation 
data from the E-mode and TE power spectra reduces this bound to r < 1.6, again with 
Hs fixed. A stronger bound has been obtained with the B-mode power spectrum by the 
BICEP experiment, giving r < 0.73 [ !'•]. The strongest bound to date on the tensor to 
scalar ratio is given when the temperature power spectrum data is also included in the 
WMAP analysis. For the pure WMAP5 data without any restriction on but with 
no spectral running the bound is r < 0.43. When BAO and SN data is combined with 
the WMAP5 data the bound on the tensor to scalar ratio becomes 

r< 0.20, (2.78) 

at the 95% confidence level. 



2.5. Non-Canonical Inflation 

In the previous section we considered the dynamics of a scalar field with a canonical 
action P((y9,X) = X — V{ip), where X = —^g'^'^d^ipduip is the kinetic energy. In this 
section we will generalise that analysis to include non-canonical actions. Non-canonical 
scalar field actions appear frequently in string theory derived inflationary models. In 
Chapters 3, 4 and 5 there are explicit examples of non-canonical scenarios. 
We will consider an action of the same form as before 

5* = y d^xy^]^ 

with minimal coupling to the gravitational sector. Varying this action gives the stress- 
energy tensor in Eq. (2.6) where = df^(p/\/2X. The energy density E is defined 

as 

E = 2XPx - P , (2.80) 

and for a homogeneous scalar field the kinetic term Pi^p, X) is the isotropic pressure. 
It proves convenient to define two parameters in terms of the kinetic function P and 



R + P{^,X) 



(2.79) 
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its derivatives [116, 174]: 



The first parameter, Cg, is called the sound speed of the fluctuations in the inflaton 
field. This can be significantly less than unity for non-canonical actions, in contrast to 
slow roll inflation driven by a canonical field such that Cg = P^x = 1- Christopherson 
& Malik showed in Ref. [17] that Cg is in fact the phase speed of the fluctuations and 
not the sound speed which is defined as P/E. However, in common with the rest of 
the literature, we will continue to use Cg as defined in Eq. (2.81). 

The generation of quantum perturbations in the non-canonical case is similar to the 
canonical one, but now includes contributions from Cg. Letting u = a6ipi, the Mukhanov 
equation for the Fourier modes, Eq. (2.48), becomes [69, 144]: 



where z has been redefined as 



z" 



z 



u{k')=0, (2.83) 



WE + P aJ2XPx 

^ - ^ ' . (2.84) 



We quantise the u modes using the Bunch-Davies vacuum as above and work with 
the amplitude w{k'^). Instead of considering whether modes are inside the comoving 
horizon, it is now important to distinguish between modes inside and outside the sound 
horizon, defined by kc^ = aH. Far inside the sound horizon, where kcg ^ aH, the 
mode solution takes a similar asymptotic form to Eq. (2.56): 

w(F) = — ^6-^^"=^ (2.85) 
v2fccg 

Following the same analysis as above, the amplitude of the curvature perturbations 
generated during inflation can be found and is given by [69] 



(2.86) 
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This expression is only valid after exit from the sound horizon. In contrast the tensor 
perturbations are not affected by the change in the action. The expression for the power 
spectrum in Eq. (2.67) is still valid. This should be evaluated after the modes have 
exited the normal horizon, i.e., when k < aH. The consistency relation (2.75) is now 
defined as [69] 



Hence, a sound speed different to unity leads to a violation of the standard inflationary 
consistency equation, which might be detectable in the foreseeable future [116, 117]. 

2.6. Non-Gaussianity 

The initial fluctuations described above have Gaussian statistics, with no correlations 
between modes on different scales. All the information about the perturbations can be 
obtained from the two-point function or power spectrum as deflned in Eq. (2.62). For a 
Gaussian random fleld all higher point functions are either zero or combinations of the 
two-point function. In particular the three-point function of 71, {7l{xi)lZ{x2)'R-{x^)) , 
will be zero for purely Gaussian TZ. We can write the Fourier transform of the three 
point function in terms of the bispectrum B [17]: 



where translation invariance imposes the conservation of the k vectors and the bis- 
pectrum depends only on the magnitude of each wavenumber. Any deviation from 
Gaussianity will result in a non-zero bispectrum value. Because of the delta-function, 
the wavevectors form triangles in Fourier space and 5 is a function of only two vari- 
ables. The bispectrum generated by inflationary models takes two main triangular 
shapes, squeezed and equilateral [13]. 

The flrst parametrisation of non-Gaussianity was deflned in real space in terms of 
the Gaussian part of the perturbation. As the non-linearity is localised in real space 
the parameter is known as the local non-Gaussianity f^^: 



r = IGcsEh = —ScsUt . 



(2.87) 



(7^(kl)7^(k2)7^(k2)) = {2nfS%k^ + k2 + ks)Bih, h, h) 



(2.88) 



(2.89) 



Here the quadratic component represents a convolution and IZq denotes the Gaussian 
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contribution [130]. We use the WMAP sign convention for /nl throughout. This is 
the opposite of the Maldacena convention: f^j^^^ = _ y^jiaWacena^ consequence of 
this choice of sign is that positive /nl imphes a decrease in temperature in the CMB 
compared to the Gaussian case. This can be seen by noting that at hnear order the 
temperature anisotropy in the CMB can be related to the curvature perturbation by 
7^G ^ -5AT/T. 

The local non-Gaussian parameter f^{] can be related to the bispectrum by: 

B{h, h, h) = ^-Jl^l [PUk,)Pl{k,) + Pl{k,)PUh) + Pl{k,)Pl{k,)\ . (2.90) 

If P^{k) is approximately scale invariant, P^{k) = ck~^, then the bispectrum becomes 
[18] 



B{k,,k„ks) = p'^lc' 



+ 7^ + 



(2.91) 



7^3 7^3 2^3 7^3 Z^3 //i3 
1 2 2 3 3 1 

This expression is maximised if one of the ki is much smaller than the other two. 
Momentum conservation then requires that they are approximately equal. This con- 
figuration is a squeezed triangle in momentum space where, for example, k^ -C ^1,^2. 
In single field inflation f^^ is proportional to the slow roll parameters and therefore 
expected to be small. Non-linear contributions from the coupling of the gravitational 
potential to the curvature perturbation are expected to produce f^l of order one which 
would be much larger than the 0{6h) contributions from single field, slow roll inflation 
[17, 104]. Any detection of /^l at greater than 0(1) would present a challenge to such 
single field slow roll models. The current bounds on the non-Gaussianity parameter are 
not strong but have been steadily tightening. The WMAP5 bound on the local form of 
/nl is 

-9</^°£<lll. (2.92) 

This observational limit still includes /^l = at the 95% confidence level. 

The other important case is where the three momenta have equal magnitude, which 
corresponds to the equilateral triangle limit. Non-Gaussianity of this shape is chiefly 
produced by models with non-canonical kinetic terms as defined in Section 2.5. The 
equilateral non-Gaussianity parameter can be evaluated in terms of the sound speed 
Cs and the A parameter defined in Eq. (2.82). The leading-order contribution to the 
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non-linearity parameter is given by [43, 174] 

f^l = f i_ _ iV 1 f i_ _ 1 _ 2aV (2.93) 

■'^^ 108 J 81 J ^ ^ 

Data from WMAP3 imposed the bound < 300 on this parameter [189]. The more 
recent WMAP5 data set improves on this bound somewhat [101], and also indicates 
that it is distinctly asymmetric. At the 95% confidence level, the current bound on the 
equilateral triangle is 

-151</^1< 253. (2.94) 

The main difference between the local and equilateral types of non-Gaussianity are 
the eras and methods of production. Local non-Gaussianity parametrises non-linear 
correlations which are local in real space. Non-linear processes taking place outside 
the horizon are the cause of these correlations. This is Production of this type of 
non-Gaussianity occurs irrespective of whether the perturbations are Gaussian when 
they cross the horizon. For single field models the magnitude of f^{] is proportional 
to the deviation of the scalar curvature power spectrum from scale invariance and is 
therefore expected to be small. On the other hand, models with multiple fields can 
produce a large amount of local non-gaussianity by the evolution of a non-infiaton 
field outside the horizon and the subsequent transfer of fluctuations in this field into 
curvature perturbations. A detection of non-negligible f^l would therefore be a very 
strong indication that multiple degrees of freedom are present in the early universe. 

In contrast, equilateral type non-Gaussianity is peaked when the momenta of the 
three modes are very similar and is generated by higher order derivative terms. Both 
the time and space derivatives become negligible once the modes have left the horizon 
and therefore any contribution to the bispectrum peaked in the equilateral shape takes 
place when the modes are inside the horizon. The extra derivative terms required are 
found generally in non-canonical models which were discussed in Section 2.5. In this 
case the amplitude of is proportional to the inverse of the sound speed squared and 
can be large. 

In the case of single field DBI infiation, discussed in Part I of this thesis, the non- 
canonical action in Eq. (3.12) contains a non-linear function of d^j^ip in the square-root 
term. These higher derivative terms are related to the magnitude of the equilateral 
type through Eq. (2.93). In the relativistic limit in which the sound speed is small, 
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can become arbitrarily large. Indeed the current observational limit on restricts the 
degree to which the relativistic limit can be reached and tighter bounds on could 
make such a limit inconsistent. 

In summary there are two main types of non-Gaussianity, which are produced in 
very different fashions^. Local non-Gaussianity is produced outside the horizon and is 
comprised of correlations which are local in real space. Equilateral non-Gaussianity is 
produced by higher derivative terms when similar modes are inside the horizon. It is 
generated by models which have non-canonical actions. 

2.7. Discussion 

In this chapter the physics of the FRW universe has been described. Inflation has been 
introduced to solve problems with the standard Big Bang scenario. To solve these 
problems the inflationary period must be of sufficient duration. This can be ensured 
by using models which comply with certain slow roll conditions. 

To explain inhomogeneities in the early universe, cosmological perturbation theory 
was presented up to first order. The power spectrum of scalar perturbations, V^, the 
spectral index of this spectrum, Ug, and the ratio of tensor-scalar perturbations, r, are 
the main observable quantities against which models can be tested. Slow roll models 
must also satisfy a consistency relation between the tilt of the tensor spectrum and r. 
Current observations favour an almost scale invariant red spectrum [ug < 1) with a 
low level of tensor signal. The accuracy of the current data is not yet good enough to 
meaningfully evaluate the slow roll consistency relation. 

As well as the standard models, one can also construct inflationary models in which 
the action takes a non-canonical form. In these models the sound speed of scalar 
fluctuations, Cg, plays a pivotal role. The predictions for scalar perturbations are altered 
by a factor of Cs, as is the slow roll consistency relation. Non-canonical models also often 
exhibit strong non-linear effects which can be parametrised using the non-Gaussianity 
parameter /nl- Canonical single field slow roll models do not predict large amounts of 
non-Gaussianity. 

In this thesis, inflation is taken to be the mechanism by which inhomogeneities in 
matter are seeded and the horizon and flatness problems of the Big Bang are solved. 

^Not all non-linear processes fit into these two categories and other types have been proposed including 
one "orthogonal" to the equilateral type [] 77] . 
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However, the inflationary paradigm is not without its own challenges. 

Chief amongst these is the lack of a unique underlying theory. Many high energy 
theories have been shown to produce an inflationary phase. Often, however, these 
require a great deal of flne-tuning in order to produce a sufficient number of e-foldings 
of inflation. Lack of knowledge about the governing physics at high energy scales 
hampers our understanding of the cause of inflation and undermines any analysis of the 
generic nature of the initial conditions required. 

The overall duration of inflation is also unknown. Observations only require that 
currently observable scales were previously inside the horizon. Thus the onset of infla- 
tion is not constrained and could occur far in the past. However, allowing such a long 
inflationary period typically increases the fine-tuning necessary and can lead to other 
issues. 

There are further problems with the infiationary paradigm, including the lack of an 
explanation for how energy in the infiaton field is transferred to the other constituent 
parts of the universe, and indeed the fact that no scalar field has yet been directly 
observed. We will continue to employ the infiationary paradigm in this thesis but it is 
important to acknowledge that some challenges remain to be overcome. 

This chapter laid the foundations for the two main parts of this work in which first 
analytic and then numerical techniques are used to constrain infiationary models. In 
the next chapter the DBI brane infiation scenario is presented. 
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3. Introduction to 

Dirac-Born-Infeld Inflation 



3.1. Introduction 

The inflationary scenario provides the theoretical framework for the early history of 
the universe. It has now been successfully tested by observations, including the five 
year data from WMAP [104]. Despite this success, however, the high energy physics 
that resulted in a phase of accelerated expansion is still not well understood. String 
and M-theory attempt to unify the fundamental interactions including gravity. The 
early universe provides a unique window into high energy physics at scales currently 
unreachable by particle accelerators. It is therefore important to develop inflationary 
models within string theory and to confront them with cosmological observations. 

One class of string theory models that has received considerable attention is D-brane 
inflation [7, 30, 33, 34, 38-40, 44, 45, 50, 57, 58, 65, 67, 85, 88, 89, 91, 93, 148, 179- 
181, 185, 200]. (For some recent reviews, see [22, 23, 49, 80, 122, 139]). The Dirac- 
Born-Infeld (DBI) scenario of the compactified type IIB theory is a well-motivated 
model [7, 185], in which inflation is driven by one or more D-branes propagating in 
a warped "throat" background. In the simplest version of the scenario, the inflaton 
parametrises the radial position in the throat of a single D3-brane. The brane dynamics 
are determined by the DBI action in such a way that the inflaton's kinetic energy is 
bounded from above by the warped brane tension. The regime where this bound is 
nearly saturated is known as the "relativistic" limit. 

In Part I of this thesis we will explore the observational consequences of DBI inflation. 
In general, primordial gravitational wave fluctuations and non-Gaussian statistics in 
the curvature perturbation provide two powerful discriminants of inflationary scenarios. 
The nature of the DBI action is such that the sound speed of fluctuations in the inflaton 
can be much less than the speed of light. This induces a large and potentially detectable 
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non-Gaussian signal in the density perturbations [7, 43, 174, 185]. 

In this chapter we introduce string theory, warped compactifications and DBI infla- 
tion. In Chapter 4 we will derive upper and lower limits on the amplitude of the tensor 
perturbations. We will explore how these bounds may be relaxed in Chapter 5 and 
discuss multi-brane scenarios which permit observable tensor signals. 

3.2. String Theory and Extra Dimensions 

The desire to unify seemingly disparate theories has been a driving force in theoretical 
physics for more than a hundred years. This effort has produced the Standard Model 
(SM) of particle physics which unifles three of the four fundamental forces in a robust 
theoretical framework. Since the realisation of the SM, a clear goal of theoretical 
physics has been the uniflcation of the fourth force — gravity — into this framework. 
String theory is one of the leading contenders for achieving this uniflcation. In this 
section we will introduce some of the string theory concepts that will be required later 
to understand DBI inflation. Many review articles and text books have been written 
about string theory and its application to cosmology and a short list of recent works 
includes Refs. [20, 49, 86, 91, 121, 139]. 

In string theory there are two main types of strings, referred to as closed and open. 
These are distinguished by the fact that closed strings form a continuous loop while 
open strings have two unconnected ends. There are several different string theories 
which are linked in pairs by a process called duality. Physical descriptions in one 
theory can be translated into a dual description in the other. The dual version often 
exhibits properties that are useful for solving problems in the original setting. We will 
work in the framework of the Type IIB theory since this has proven to be the most 
fruitful for generating models of cosmological inflation [49, 121]. 

3.2.1. Extra Dimensions 

String theory predicts that the one time-like and three spatial dimensions that consti- 
tute the observable universe do not represent the complete spacetime manifold. Instead, 
our universe is a 10 or 11 dimensional spacetime and physical theories in 3+1 dimensions 
must therefore be able to explain why the other 6 or 7 dimensions are unobservable. 
One of the challenges of string theory is how to "hide" these extra dimensions in such 
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a way as to recover the standard four- dimensional cosmology at low energies. 

The early work of Kaluza and Klein (KK) in formulating higher dimensional theories 
laid the groundwork for the current treatment of extra dimensions in string theory 
[92, 99]. By compactifying an extra dimension onto a circle of finite radius an infinite 
tower of extra fields are introduced into the lower dimensional theory. The mass of 
these fields is inversely proportional to the size of the extra dimension. The appearance 
of these unobserved massive fields is avoided by taking the radius to be extremely 
small. This leaves a massless degree of freedom which must be accounted for in the 
four-dimensional effective action. 

A similar procedure is undertaken when compactifying string theory from a 10 or 11 
dimensional description down to four dimensions (for reviews see Refs. [56, 73]). In ten 
dimensional type IIB theory the six extra dimensions are compactified into a Ricci flat 
Calabi-Yau (CY) manifold which can be described by three complex coordinates [206]. 
Because any Ricci flat metric can be rescaled onto another Ricci fiat metric, there is 
no unique solution for the metric on the CY manifold. Instead a family of solutions 
exists with many free parameters. These parameters remain after the compactifica- 
tion, in analogy to the size of the extra dimension in KK compactification, and can 
depend on position in the four-dimensional spacetime. They appear as fields in the 
four- dimensional theory and are known as moduli. These fields are not subject to any 
symmetry and so their individual values at different spacetime points can affect the 
physics at those points. 

3.2.2. T-duality 

In string theory an extra space time symmetry is present which relates physical prop- 
erties in theories with large compactification radius with those in theories with small 
radius. Suppose we have a string theory compactified on a circle of radius L. The "T- 
duality" transformation which relates two physical theories with this one compactified 
dimension is 

~ a' 

L^L = ^. (3.1) 

Now consider what effect this transformation will have on the momentum of a closed 
string. Instead of being a continuum, the momentum takes discrete values j/L for 
j G Z. This is a KK tower of massive states. As we complete a circuit around the 
compact dimension, the value of the coordinate function embedding the string in the 
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background will increase by 2t:wL for w G Z. This w, called the winding number of 
the string, can only be non-zero for closed strings, which can be wrapped around the 
periodic dimension. 

The total mass of the string contains terms with both the KK tower of states and 
the new tower of winding states: 

M^ = f^ + ^ + ---, (3.2) 

where the string parameter a' is related to the string mass scale by a' = If L 

is taken to infinity, the w states become infinitely massive and only the w = 
state is left with a continuum of momentum values. Thus, the uncompactified result 
is recovered. However, if L — )■ 0, the j ^ states are now infinitely massive as in the 
standard KK picture. Unlike the standard case, there is now a continuum of winding 
states with w ^ 0, again giving an uncompactified dimension. This major departure 
from the standard compactification result is a purely stringy effect. 

The formula for the mass spectrum, Eq. (3.2), is invariant when j and w are exchanged 
given the transformation in Eq. (3.1). Writing the equations of motion in terms of L, 
having interchanged j and w, gives a new theory which is compactified on a circle of 
radius L. This is known as the T-dual theory [94, 167]. The two theories are physically 
identical since T-duality is an exact symmetry of string theory for closed strings. The 
T-duality applies to all physics in the theory and in particular also affects open string 
modes. These behave in a different way under T-duality to closed strings as will be 
described below. 



3.2.3. D-Branes 

The dynamics of extended objects known as branes are particularly important for build- 
ing inflationary models. As mentioned in Section 3.2.2, string theories are linked by 
T-duality. Fundamental parameters such as the size of the extra-dimensions, the string 
coupling and the coordinate solutions of the strings are related by such a symmetry. 

We introduced T-duality by explaining its effects on closed strings. But what happens 
to the open strings in a T-dualised theory? Open strings, as their name suggests, have 
two open ends and consequently cannot have a conserved winding number such as w. 
Suppose once more that one of the D dimensions is compactified. As L — )■ 0, the non- 
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zero momentum states become infinitely massive, but in contrast to the closed case 
there is now no continuum of winding states. Thus, the open string now lives in D — 1 
dimensions similar to the result of standard KK compactification [ ]. The endpoints 
of the open strings then observe Dirichlet boundary conditions, taking fixed values in 
the compactified direction. There are still closed strings in this theory, however, and 
these continue to move in the full D dimensions after being T-dualised. 

The result is similar if more than one coordinate is made periodic. If D — p — 1 
spatial dimensions are compactified, for some p, then the ends of the open strings can 
still move freely in the other p spatial dimensions on a p + 1 dimensional hypersurface. 
This hypersurface is called a Dirichlet brane or Dp-brane. The closed string modes 
move in the full D dimensions. In Type IIB theory with supersymmetric strings, an 
extra condition implies that only Dp-branes with p = 1,3, 5, 7, 9 are stable^. 

Dp-branes can be considered as dynamical objects in their own right with a tension 
given by [bb]^ 

Tp = — , (3.3) 

where Qs is the string coupling and rUs is the string mass scale. Their dynamics is 
governed by the action 



pjd^^'^^/^g, (3.4) 



where cjab is the induced metric on the brane with internal coordinates for a = 
0,...,p, given by [86] 

Eq. (3.4) is the general form of the DBI action which will be used later. 

In the simplest versions of slow roll inflation, only a single scalar field with a suffi- 
ciently fiat potential is required to satisfy the slow roll conditions outlined in Section 
2.2.2. Since D-branes are charged (with Ramond-Ramond charge), a D-brane and an 
anti-D-brane (D) separated by some distance will be attracted to each other. The sep- 
aration distance can be identified as a scalar degree of freedom and under appropriate 
conditions could play the role of the infiaton field [30, 33, 34, 57, 58, 170]. 

As described above, compactifying dimensions introduces scalar fields known as mod- 

^ There is also a p = — 1 D-instanton in which the time direction along with all spatial directions is 

subject to Dirichlet boundary conditions [74-76]. 
^Tp here is Tp in Ref. 
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uli. These fields must be accounted for in the dynamics unless some way can be found 
to stabilise them by fixing their masses to be large. Initial efforts to induce infiation 
using D-branes ignored the issue of moduli stabilisation. Instead, it was assumed that 
whatever stabilisation mechanism was used would have no discernible effect on the infla- 
tionary physics. Kachru et al. [ ] recognised that in fact stabilisation will be important 
and must be taken into account. 

3.2.4. Warped Throats 

The moduli must be stabilised so that they do not appear in the effective action as 
massless flelds. This can be achieved by switching on background fluxes in the com- 
pactified space. These fiuxes are analogous to magnetic fields in the higher dimensional 
space. By Gauss' theorem the compact space will now have a quantised non-zero total 
charge. In the presence of fiuxes, a general form for the ten dimensional metric is [20]: 

ds^ = e^^^y^^Ax'^dx" + e-^^^y^g^ndy'^dy^ , (3.6) 

where the function Aijj) varies across the compact dimensions y*". Compactifications 
in which A varies significantly with y are called warped compactifications and e"^'-^-' is 
referred to as the warp factor. These warped compactifications are qualitatively similar 
to the Randall-Sundrum scenario [31, 160]. 

Flux compactification of type IIB string theory to four dimensions results in such a 
warped geometry, where the six-dimensional CY manifold contains one or more throats 
[56, 71, 73]. The metric inside a throat takes the same form as in Eq. (3.6): 

dsl, = h\p)dsl + h'\p) (dp2 + pM4J , (3.7) 

where the warp factor h{p) is a function of the radial coordinate p along the throat and 
X5 is a Sasaki-Einstein five-manifold. 

In many cases, the ten-dimensional metric (3.7) can be approximated locally by the 
geometry AdS^ x X5, where the warp factor is given by h = p/L and the radius of 
curvature of the AdS^ space is defined by 



47r\N 
Vol(X5)m4 ' 



(3.8) 
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such that Vol(X5) is the dimensionless volume of X5 with unit radius and N is the D3 
charge of the throat. 

In the Klebanov-Strassler (KS) background [07], the throat is a warped deformed 
conifold and corresponds to a cone over the manifold X5 = T^'^ = SU(2) ® SU(2)/U(1) 
in the UV limit (p — )■ 00). This has a volume Vol(T^'^) = 167r'^/27 and topology x S^, 
where the S'^ is fibred over the S^. 



Figure 3.1.: A conifold can be deformed to remove the singularity at the tip. 

There are two 3-cycles in the warped throat. The first is the subspace and is known 
as the A-cycle. The second, called the B-cycle, is the 5*^ times a circle extended in the 
direction of the throat radius. The three-form fluxes F3 and H^, aligned with these 
cycles, are turned on to make the warped throat a solution of the Einstein equations 
[49]. The cycles are threaded with quantised units of flux M and K given by: 



where M, K ^ X. The D3 charge of the throat, N, is related to the quantised fluxes 
hy N = MK. The wrapping of the fluxes along the cycles of the conifold smooths out 
the conical singularity at the tip of the throat with an cap [97, 98], as shown in 
Figure 3.1, and the warp factor asymptotes to a constant value in this region. 

In this section we have summarised the concepts that will be required to discuss 
DBI inflation. The compactified warped throat described here will provide the setting 





(3.10) 



(3.9) 
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for tliis string tlieoretic realisation of inflation. In the next section we connect the 
geometry and physics of the string compactification with inflationary cosmology and 
establish the observational parameters that will directly enable concrete constraints to 
be formulated. 



The DBI scenario is based on the compactification of type IIB string theory on a Calabi- 
Yau (CY) three-fold, where the form-field fluxes generate locally warped regions known 
as throats, as described above. The propagation of a D3-brane in such a region can 
drive inflation, where the inflaton field is identified with the radial position of the brane 
along the throat. Inflation can occur whether the brane moves towards or away from 
the tip of the throat. Since the radial distance is an open string mode, the field equation 
for the inflaton is determined by a DBI action. 

In general, the low-energy world-volume dynamics of a probe D3-brane in a warped 
throat is determined by an effective, four-dimensional DBI action, as described in Sec- 
tion 3.2.4. The inflaton field is related to the radial position of the brane hj if = y/TsP, 
where T3 is the brane tension defined in Eq. (3.3). The action is then given by [185] 



where R is the Ricci curvature scalar and T{ip) = T^h'^^ip) defines the warped brane 
tension. As in Section 2.5 we refer to P{if,X) as the kinetic function for the inflaton, 
X = —\g^^^V is the kinetic energy of the inflaton and V{ip) denotes the field's 
interaction potential. Typically in warped compactifications of IIB supergravity, this 
potential is determined by the relevant fluxes and brane interaction terms. We will 
ignore the precise origin and form of this potential, but simply note that it is highly 
sensitive to the string theoretic construction. For the purpose of this thesis we will 
simply treat it as an arbitrary function of the inflaton field. (See, for example, Ref. [88] 
for a discussion on the precise form that the inflaton potential may take.) 

We consider a spatially flat and isotropic cosmology sourced by a homogeneous scalar 
field. In this case, the Friedmann equations for a monotonically varying inflaton can 
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P(^, X) = -T{ip) v/l-2T-i(^)X + T(v9) - V{ip) , 



(3.12) 
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Figure 3.2.: A representation of tlic Calabi-Yau manifold in the 6 compactificd dimen- 
sions. Throats are connected to the main bulk. D3-branes appear as dots. 
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be expressed in the form [185] 



?,Ml^H\^) = V{^) - T{^) [l - ^1 + AM^^T-^Hl 



(3.13) 



(3.14) 



^1 + AM^^T-^H%' 



In Section 2.5 we introduced the speed of sound of inflaton fluctuations. For the 
kinetic function in Eq. (3.12), we find from Eq. (2.81) that 



The condition that the sound speed be real imposes an upper bound on the kinetic 
energy of the inflaton, 0^ < T{ip), which is independent of the steepness of the potential. 
The motion of the brane is said to be relativistic when this bound is close to saturation. 
We will assume throughout Part I of this thesis that motion takes place in the relativistic 
limit in which Cg ^ 1. 

We now define the epoch that is directly accessible to cosmological observations 
as "observable inflation". We will assume that this phase occurred when the brane 
was located within a throat region and moving towards the tip of the throat. We 
denote the parameter values evaluated during observable inflation by a subscript star 
( =„). Observable inflation corresponds to about 4 e-foldings of inflationary expansion, 
ATV* ~ 4, and occurred somewhere between 30 to 60 e-foldings before the end of 
inflation. 

The definitions of the slow roll parameters defined in Section 2.2.2 change when Cg is 
not equal to unity and we will include a third parameter, s, which quantifies the rate 
of change of Cg. The inflationary dynamics during this phase can be quantified in terms 
of these three parameters: 



1 



Vi-2r-ix. 



(3.15) 
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where 7 = l/cg. We will assume that the quasi-de Sitter conditions {eh, \vh\, \ s\} ^ 1 
apply during observable inflation. In this regime, the amplitudes and spectral indices 
of the two-point functions for the scalar and tensor perturbations are given by 



2 1 

r, 2 

^T = -,J7^. 3.20 

=4eH-2r7H + 2s, (3.21) 
nt = -2eH , (3.22) 

respectively. and rit are evaluated when k = aH but and are evaluated when 
the scale with wavenumber k crosses the sound horizon fccg = aH. 

A further important consequence of a small sound speed is that departures from 
purely Gaussian statistics may be large [7, 43, 174, 185]. DBI inflation produces non- 
Gaussianity maximised in the equilateral configuration and the leading contribution is 
in the form of Eq. (2.93). When c^P^x = 1 the second term in Eq. (2.93) is identically 
zero and becomes [43, 117] 

a- 4(^-1). (3.23) 

When Cs <^ 1 a significant level of non-Gaussianity is produced. For a homogeneous 
field 2X = 0^, so from Eq. (3.15) we find that 

02 = T(^)(l-a. (3.24) 

Eqs. (3.19), (3.23) and (3.24) may then be combined to provide a relation for the 
warped brane tension: 



T(^) vr^ 2 



16 '^V 3/, 



'NL 



^K{^-VJ^] ■ (3.25) 



3.4. The Lyth Bound 

In the next two chapters we will use a powerful result due to Lyth [123]. This links the 
change in value of the infiaton field during inflation to the production of tensor modes. 
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This relation was originally derived for canonical actions but can be straightforwardly 
extended to the case of non-canonical actions such as the DBI action. 

Eqs. (2.86) and (2.67) imply that the variation of the inflaton field during inflation 
is related to the tensor-scalar ratio by [19, 123] 



1 /d(^ 



(3.26) 



where A/" denotes the number of e-foldings as defined in Eq. (2.17). The total variation 
in the inflaton field between the epoch of observable inflation and the end of inflation 
is then given by 



MpL \8CsP,X 



1/2 

A/'cff, (3.27) 



where 



If r/{csP^x) varies sufficiently slowly during observable inflation, the corresponding 
change in the value of the inflaton field is given approximately by [19, 123] 



Av^V (AAC: 



This equality links the total variation of the inflaton during observable inflation with 
the tensor-scalar ratio, i.e., the amplitude of gravitational waves produced during that 
period. In Chapter 4 we will show how the dynamics of the DBI scenario allow an 
upper limit to be imposed on r using this relation. 

In deriving Eq. (3.29) we have assumed that r/cgPx varies slowly during observable 
inflation. For the DBI case, CgP^x = 1 and the change in r can be related to the change 
in en and Cg through Eq. (2.87). As we have taken en, \vh\, \s\ ^ 1 the tensor-scalar 
ratio will indeed vary slowly over the observable epoch. For more general models where 
CsP,x 7^ 1 we have that 

16^ i2eH - 2r]H) • (3.30) 



dAf 
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Therefore rjc^P^x varies slowly as long as Px is not too small, i.e., close to O^ejj). 
This will not be the case in the models studied in Chapters 4 and 5. 
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3.5. Discussion 

In this chapter we have introduced the Dirac-Born-Infeld inflationary scenario. Many 
attempts have been made to provide an inflationary expansion phase in the early uni- 
verse using string theory. In compactifying from ten dimensions down to four, compli- 
cated geometries and additional fluxes must be used to stabilise the remaining moduli 
fields. 

The DBI scenario is a particular example of the non-canonical inflationary para- 
digm described in Section 2.5. The radial position of a D3-brane in a warped throat 
is identified as the inflaton field. While the brane propagates up or down the throat, 
the kinetic energy of the inflaton is bounded above by requiring the sound speed of 
fluctuations to be real. This bound holds no matter how steep the potential of the 
field. The relativistic limit takes the bound to be close to saturation and the sound 
speed to be small. In the case of DBI inflation the speed of sound parameter takes 
the simple form Cg = 1/-P.X- The previously derived results for and n^, as well as 
the redefined slow roll parameters (3.16)-(3.18) can then be expressed in terms of this 
parameter. 

Significant non-Gaussianity in the density perturbations spectrum can be generated 
due to the small sound speed of the inflaton fluctuations. This non-Gaussianity can 
be related to the brane tension and tensor-scalar ratio through Eq. (3.25). The tensor- 
scalar ratio can also by related to the variation in the inflaton field by the Lyth bound 
(3.26). This relation can be refined by focusing only on the period of observable infla- 
tion. In the next chapter we will derive two competing bounds on r which will strongly 
constrain the parameter space for DBI models. 



4. Observational Bounds on DBI 
Inflation 



4.1. Introduction 

In this chapter two bounds on the amphtude of primordial gravitational waves will be 
derived, which severely challenge the standard DBI inflationary scenario. By consider- 
ing the field range of observable inflation inside a warped throat, the tensor-scalar ratio 
r will be constrained to be less than 10~^. In contrast a lower bound of r > 0.005 will 
be derived when the power spectrum of scalar perturbations has a red spectral index. 
These clearly incompatible bounds can be relaxed by using a more general form of the 
DBI action. 

The gravitational wave background generated in DBI inflation was initially investi- 
gated by Baumann & McAllister (BM) [19]. By exploiting a relationship due originally 
to Lyth [123], these authors derived a fleld-theoretic upper limit to the tensor ampli- 
tude and concluded that rather stringent conditions would need to be satisfled for these 
perturbations to be detectable. Moreover, the special case of DBI inflation driven by 
a quadratic potential is incompatible with the WMAP3 data when this constraint is 
imposed [, ' =]. 

Our aim in this chapter is to derive observational constraints on DBI inflation that are 
insensitive to the details of the throat geometry and the inflaton potential. In general, 
there are two realisations of the scenario, which are referred to as the ultra-violet (UV) 
and infra-red (IR) versions. These are characterised respectively by whether the brane 
is moving towards or away from the tip of the throat. We focus initially on the UV 
scenario and derive an upper bound on the gravitational wave amplitude in terms of 
observable parameters. This limit arises by considering the variation of the inflaton 
fleld during the era when observable scales cross the Hubble radius, and we flnd in 
general that the tensor-scalar ratio must satisfy r < 10~^. This is below the projected 
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sensitivity of future CMB polarisation experiments [21, 199]. 

On the other hand, the WMAP5 data favours a red perturbation spectrum, with 
Ug < 1, when the scalar spectral index is effectively constant [lUI]. For models which 
generate such a spectrum, we identify a corresponding lower limit on the tensor modes 
such that r > 0.1(1 — n^). This is incompatible with the upper bound on r when 
1 — Us — 0.03, as inferred by the observations. 

Therefore a reconciliation between theory and observation requires either a relaxation 
of the upper limit on r or a blue spectral index [us > 1). The DEI scenario would need 
to be generalised in a suitable way for the upper bound on r to be weakened. Necessary 
conditions are identified that a generalised action must satisfy for the BM constraint 
and our newly derived bound to be relaxed. Such conditions are shown in Chapter 5 
to be realised in a recently proposed IR version of DEI inflation driven by multiple 
coincident branes [19 1]. 

4.2. An Upper Bound on the Primordial 
Gravitational Waves 

In Ref. [19] Baumann & McAllister derived a field-theoretic upper bound on the tensor- 
scalar ratio. They achieved this by noting that the four-dimensional Planck mass is 
related to the volume of the compactified CY manifold, Vg, such that M|l = Vqk^q, 
where k^q = ^{2iTy g'^m~^ = vr/T| for a D3-brane^. In general, the compactified volume 
is comprised of bulk and throat contributions, Vq = Ve bulk + Ve throat- The latter is given 



where puv denotes the radial coordinate at the edge of the throat (defined as the region 
where h{puv) is of order unity). The geometry of the throat is shown in Figure 4.1. 

If one assumes that the bulk volume is non-negligible relative to that of the throat 
(^6 throat < Ve), it follows that M|l > Vq 

throat 'tio^- ^ warpcd AdS^ x geometry, 
this leads to an upper limit on the total variation of the infiaton field in the throat 



by 




(4.1) 



^We parametrise the Planck scale in terms of the D3-brane tension out of convenience, and note that 
there is no physical relationship between the two. 
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Figure 4.1.: Geometry of the warped throat. The radial coordinate p is measured from 
the tip of the throat to puv at the join with the bulk manifold. 
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region in terms of the D3 charge: 



< ^ ■ (4-2) 



Condition (4.2) may be converted into a corresponding limit on the tensor-scalar 
ratio by noting from the definition (3.16) that 0^/M|l = '^BhH^ / P,x- This implies 
that the variation of the inflaton field is given by the Lyth bound (3.26) [19, 123]: 



(4.3) 



where Af is the number of e- foldings as defined in Eq. (2.17). Since v?*, the field value 
during observable inflation, is less than (puv, this results in an upper bound on the 
observable tensor-scalar ratio [19]: 

32 , , 

n < ,o • (4.4) 



iV(A4ff, 



The effective number of e-foldings, A/'efr, defined in Eq. (3.28), is a model-dependent 
parameter that quantifies how r varies during the final stages of inflation. Since CgPx = 
1 in the standard DBI model, it follows that A/'efr = A/'end if ^ is constant during inflation, 
where Mend is the total number of e-foldings from the epoch of observable inflation until 
inflation ends. 

Typically, one expects 30 < A/'cs < 60, although smaller values may be possible if 
the slow roll conditions are violated after observable scales have crossed the horizon. 
Furthermore, A^ ^ 1 is necessary for backreaction effects to be negligible [19]. Hence, 
the constraint (4.4) imposes a strong restriction on DBI inflationary models. On the 
other hand, the numerical value of A^fr is uncertain. Our aim here is to focus on the 
range of values covered by the inflaton field during the observable stages of inflation. 
This will result in a constraint on the tensor modes that can be expressed in terms of 
observable parameters. 

To proceed, we denote the change in the value of the inflaton field over observable 
scales by Ayj,, = a/7^^P*- Since the brane moves towards the tip of the throat in UV 
DBI inflation, it follows that > Pend > 0, which implies that 

P*>|Ap,|. (4.5) 



4.2: An Upper Bound on the Primordial Gravitational Waves 



59 



This change in the inflaton value will correspond to a fraction of the throat volume, 
lAVg*! < Vg throat ^ Ve, where equality in the second limit arises if the bulk volume is 
negligible. Hence, |Ay9^,| is bounded such that 

^ < • (4-6) 



The observations of the CMB that directly constrain the primordial tensor pertur- 
bations only cover multipole values in the range 2 < Z < 100. This is equivalent to 
AAC ~ 4 e-foldings of inflationary expansion and, in general, corresponds to a narrow 
range of inflaton values. To a first approximation, therefore, the fraction of the throat 
volume (4.1) that is accessible to cosmological observation can be estimated to be 

|AV6.|^Vol(X5)^^. (4.7) 
Combining the inequality (4.5) with Eq. (4.7) then implies that 

|AV6*|>Vol(X5)^^. (4.8) 
Substituting the condition (4.8) into the bound (4.6) gives 



VMpL^ Vol(X5) VMpL 
and using T((^) = Tjjv" and Eq. (3.25) yields the upper limit 

Hence, employing the Lyth bound (4.3) in the form (A(y9^,/MpL)^ ~ r(AAC)^/8 results 
in a very general upper limit on the tensor-scalar ratio: 

327r3 , / 1 \ , , 



(AA4)^voi(X5) V ^ir 



Condition (4.11) is only weakly dependent on the level of non-Gaussianity when 
"/nl ^ ^ may therefore neglect the factor involving this parameter. Substitut- 
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ing the WMAP5 normalisation ~ 2.5 x 10 then implies that 

2.5 X 10-^ 

Furthermore, the most optimistic estimate for the minimum number of e-foldings that 
could be probed by observation is AA/'^, ~ 1, whereas a generic compactification arises 
when the volume of the Einstein five-manifold is Vol(X5) ~ 0{7r^) [97]. This yields 
a model-independent upper bound on the tensor-scalar ratio for standard UV DBI 
inflation: 

< 10~^ (4.13) 

The bound (4.13) is the main result of this section. This value of r is significantly below 
the sensitivity of future CMB polarisation experiments, which will measure r > 10~^ 
[21, 199]. If CMB observations are able to span the full range of e-foldings such that 
AA/"* ~ 4, this constraint is strengthened to r,, < 2 x 10^"'^"'^. 

Before concluding this section, we should explicitly outline all the assumptions that 
have lead to Eq. (4.13). First, we are considering the relativistic limit where Cg <^ 1. 
We are also restricting ourselves to considering the UV scenario where a brane moves 
towards the tip of the throat. This ensures that Eq. (4.5) is satisfied. For the Lyth 
bound to take the form in Eq. (3.29), we have assumed that r varies slowly during the 
observable period of inflation. This is justified as the change in r can be written in terms 
of the quasi deSitter parameters eH,VH and s and we have assumed their magnitudes 
are much less than unity. 

The estimate (4.7) was derived under the assumption that the integrand in Eq. (4.1) 
is constant. This inevitably introduces errors into the bound (4.11). However, the two 
limiting cases of interest in KS-type geometries arise when the warp factor scales either 
as /i oc p or as ~ constant [97, 98]. In both cases the integral (4.1) can be performed 
analytically. Indeed, if we specify h cc p'^ for some constant a, evaluate the integral 
from to + Ap*, and expand to second-order in a Taylor series, we find that 



AFe.-Vol(X,)^(Ap 



P* I (5-4c .)(Ap,) 



P* 



(4.14) 



This implies that the error in Eq. (4.7) is no greater than about 3{Ap^,/p^,) if < a < 1. 
More generally, it follows that a similar error will arise for any warp factor h oc p°'^''\ 
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where the function a{p) satisfies < a(p) < 1 over observable scales. We conclude, 
therefore, that Eq. (4.7) provides a sufficiently good estimate of the volume element for 
a generic warp factor^. 

In order to neglect the /^l term in Eq. (4.10) we have assumed that — /nl > 5. 
As Cg has been taken to be small this is expected to be the case. The volume of 
the Sasaki-Einstein manifold is taken to be 0{t[^) in keeping with the values for 
known solutions. The WMAP5 normalisation of the scalar perturbation power spec- 
trum has also been used. Finally, in going from Eq. (4.12) to the final numerical figure 
in Eq. (4.13) the most "optimistic" value, AA/'^, ~ 1, has been chosen as this leads to 
the least restrictive bound on r. As described above a more realistic value of 4 would 
severely constrain r due to the strong dependence of Eq. (4.12) on AA/"*. 

4.3. A Lower Bound on the Primordial 
Gravitational Waves 

The analysis of the previous section indicates that standard versions of UV DBI infiation 
generate a tensor spectrum that is unobservably small. Therefore, r = can be assumed 
as a prior when discussing the WMAP5 data. However, in this case the data disfavours 
a scale-invariant density spectrum at close to the 3(T level (2.78cr) when the running 
in the spectral index, = dris/dlnfc, is negligible [101]. Furthermore, a blue spectral 
index is only marginally consistent with the data when r 7^ and as = 0. (The inferred 
upper limit is < 1.018.) Although the results from WMAP5 do allow for a blue 
spectrum if there is significant negative running in the spectral index, we will focus 
in this section on models that generate a red spectral index < 1, since these are 
preferred by the current data. 

In general, the spectral index may be related to the tensor-scalar ratio. After differ- 
entiating Eq. (3.15) with respect to coordinate time, and employing Eqs. (3.14) and 
(3.21), we find that^ 



1 - = Aeh + 



1 — 7' 



2s 



2M^^TJH, 
7 TH 



(4.15) 



^As we shall see in the following section, even an order of magnitude error will make little difference 
to our final conclusions. 
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where the minus (plus) sign corresponds to a brane moving down (up) the warped 
throat. The second term in Eq. (4.15) can be converted into observable parameters 
by defining the 'tilt' of the non-linearity parameter [10]: 

d In I f*""* I 

This implies that s = 3f^n'^^/[2{l — 3/^l)] and substitution of Eqs. (3.21)-(3.23) into 
Eq. (4.15) then yields 

eq 



l — ^yr^+^T/^^MP.)^. (4.17) 

In Ref. [117], brane infiation near the tip of a KS-type throat was considered, where 
the warped brane tension asymptotes to a constant value. In this regime, Eq. (4.17) 
reduces to the condition r ~ 2.3(1 — Us)/ a/— /^l when is sufficiently large to be 
detectable by Planck, i.e., 1/^11 > 5. It then follows from the WMAP5 best-fit value 
Ug ~ 0.968 and lower limit > —151 [ ] that the gravitational wave amplitude is 
bounded both from above and below such that 0.001 < r < 0.01. These bounds follow 
from current WMAP5 limits on the spectral index and the non-linearity parameter, but 
do not take into account the field-theoretic upper bound that must be imposed on the 
variation of the infiaton field during infiation. 

More generally, in UV DBI infiation where the brane moves towards the tip of the 
throat, it is reasonable to assume that the warp factor decreases monotonically with 
the radial coordinate over the observable range of infiaton values, i.e., dh/dp > 0. 
This condition is satisfied for AdS^ x X5 compactifications and KS-type solutions. 
Consequently, the third term in Eq. (4.17) will be semi-negative definite, which implies 
that 



^ ^' - pcq I ^NL 

4V* "^^^^1-3/^L 



3/Nl + r^^>l-^.- (4.18) 



Condition (4.18) is a consistency relation on UV DBI infiation in terms of observable 
parameters and it may be combined with the upper bound (4.11) to confront the sce- 
nario with observations. Firstly, let us assume that the tensor-scalar ratio is negligible. 

^The relationship between ip and H^^p is defined in Eq. (3.14). When 1^ < and the brane moves 
down the throat (UV case) H ^p > 0. Alternatively in the IR case when cp > we have H < 0. In 
order to remove the ambiguity we rewrite Eq. (3.14) using —H,^ = ^\H,^\ where the minus (plus) 
sign corresponds to the UV (IR) case. 
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The WMAP5 data implies that 1 — > 0.026 at Icr, and this is only compatible with 
condition (4.18) if 

^ -2s > -3(1 - n,)f^l > -0.078/^1 . (4.19) 

However, when — /^l ^ this would violate the slow roll conditions that must be 
satisfied for a consistent derivation of the perturbation spectra (3.19). For example, 
the conservative bound \s\ < 0.1 with 1 — ~ 0.05 is violated if — /nl > ^{^)- 

In view of this, let us consider the case where the tensor perturbations are non- 
negligible. The magnitude of the second term in condition (4.18) is suppressed by a 
factor of (— /nl)'^''^ ^ 1 relative to the first. This is expected to be a significant effect 
in DBI infiation. Consequently, by saturating the WMAP5 limit > —151, we arrive 
at a lower bound on the tensor-scalar ratio which applies to any model for which the 
ratio |nNl//Nll is negligible: 

n > > ^ . (4.20) 

This second bound requires r > 0.005 for the WMAP5 best-fit value 1 — ~ 0.032, 
which is incompatible with the upper limit (4.13). 

In general, therefore, it is difficult to simultaneously satisfy the bounds on r with the 
WMAP5 data in standard UV DBI infiation. There is a small observational window 
where a blue spectrum is consistent with the data, in which case the lower limit (4.20) 
does not apply. However, if the tensor modes are negligible, as implied by the inequality 
(4.13), the data strongly favours a red spectral index with Ug < 0.974, and this violates 
the condition (4.20). A significant detection of a red spectral index requires either a 
violation of the slow roll conditions or a sufficiently small value for the volume of X5. 
In particular, combining the limits (4.12) and (4.20) results in the condition 

9x1 0~^ 

and we find that Vol(X5) < 10^'^ is required for typical values 1 — Ug — 0.05 and 
AA/"* ~ 4. This is comparable to the limit on the volume derived for the special case of 
a quadratic inflaton potential [19]. 

As noted in Refs. [19] and [23], condition (4.21) may be achieved if X5 corresponds 
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to a Y^''^ space. Previously only two five- dimensional Sasaki-Einstein metrics were 
explicitly known, and T^'^ on x S^. The Y^''^ metrics described in Ref. [TD] are 
a countably infinite number of Sasaki-Einstein metrics on S"^ x 5*^. The metrics are 
parametrised by the two topological numbers p and q, which are coprime when the 
y^''^ is topologically S"^ x S^. The volume of one of these manifolds is proportional 
to 1/p. Hence by setting q = 1 and letting p become large, this volume can be made 
arbitrarily small On the other hand, the largest volume occurs for p = 2, q = 1 

giving Vol(y^'^) 0.2971'^. Small volumes could also be realised by orbifolding the S'^ 
symmetry of a KS-type throat. 

On the other hand, the upper limit (4.12) on the gravitational waves follows as a 
consequence of assuming the constraint (4.5). This could be violated in IR versions of 
the scenario, where observable scales crossed the Hubble radius when the brane was 
near the tip of the throat and (f ^ MpL [39, 10]. Nonetheless, we emphasise that the 
upper bound (4.12) on the tensor modes will also apply to any IR DBI model for which 
|A(y9^,| In view of the above discussion, we will proceed in the following section 

to discuss a framework for generalising the DBI scenario so that the constraints on the 
tensor modes can be satisfied. 

4.4. Relaxing the Upper Bounds 

In this section, we take a phenomenological approach and consider the following kinetic 
function which has a more general form than the DBI one but still contains a square 
root term: 

P = -f^(^)^l-f^(^)X- fciv) , (4.22) 

where /i(v?) are unspecified functions of the infiaton field. We will assume implicitly 
that these functions have a suitable form for generating a successful phase of infiation. 
A direct comparison with Eq. (3.12) indicates that the standard DBI action can be 
recovered by setting /a/b = 2. This implies that CgP^x = 1 and greatly simplifies the 
form of Eq. (4.3). Another important property in the DBI case is that the warp factor 
uniquely determines the kinetic structure of the action, i.e., /i^ oc oc /^^ In view 
of this, it is interesting to consider whether the gravitational wave constraints could be 
weakened by relaxing one or both of these conditions. 
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We can differentiate P(X, ip) in Eq. (4.22) to find: 

The sound speed of fluctuations in the inflaton, defined in Eq. (2.81), is then given by 

Cs = ^/l^hX = ^^^. (4.25) 

and the scalar power spectrum (2.86) by 

o 1 H'^ 

However, the consistency equation (2.87) and non-Gaussianity constraint (3.23) remain 
unaltered for this more general class of models [117]. It follows, therefore, that the 
CMB normalisation condition (3.25): 

generalises to a constraint on the value of /^(yj,,,): 



_fA_\ _vr^ 



^ _^2p2 M _ __ 1 (4_28) 



Finally, the expression for the scalar spectral index follows by generalising the derivation 
of Eq. (4.17) given by 



It is straightforward to show that for the more general kinetic function this expression 
becomes 

i-„..^;rr^.^,y^(^Mr.) . ,4.30, 
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4.4.1. A More General BM Bound 

The BM bound (4.4) restricts the maximal variation of the scalar field ip in the full 
throat region for DBI infiation. This is determined by expression (4.2) for generic 
warped geometries that are asymptotically AdS^ x X5 away from the tip of the throat. 
However, in Section 3.4 the Lyth bound was also defined for general non-canonical 
actions. For the more general kinetic function (4.22), the BM bound becomes 

32 16 

To use this bound we must be able to calculate A/'eg over the full range of e-foldings of 
inflation. This requires knowledge of the behaviour of and fs over that range. 

A more cautious approach would be to restrict our attention to the observable stage 
of inflation. Assuming that the variation of /a/b = 2csPx is negligible during that 
epoch, we can use Eq. (3.29) which states that 



Mpl;, 8 \c,P,xJ, 4 V/a/ 



Tb 



(4.32) 



In addition, if observable scales leave the horizon while the brane is inside the throat, 
the change in the fleld value must satisfy |A(y9*| < ipuv- It follows from Eqs. (4.31) and 
(4.32), therefore, that 

32 16 

< iv(w^'^^'^^* = iv(AAg^^^^^^^* • ^^-^^^ 

Condition (4.33) will be referred to as the generalised BM bound. We have been 
conservative by restricting our discussion to the observable phase of inflation. A stronger 
condition is obtained by using Eq. (4.31), which is equivalent to substituting AAA* — )■ 
ACflf. If /a/b remains nearly constant over the last M e-foldings of inflation, then Mes 
may be as large as 60 and the right hand side of Eq. (4.33) will be reduced by a factor 
of 225. Thus, the generalised bound (4.33) should be regarded as a necessary (but not 
sufficient) condition to be satisfied by the tensor modes. 

Given expressions (4.28) and (4.33) we can either constrain r using a specified value 
for fsi or find a necessary condition on the value of fsif*) for the generalised BM 
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bound to be satisfied using r and V^: 



In UV models, identical arguments that led to the lower limit (4.20) on the tensor- 
scalar ratio will also apply in this more general context if, as expected, /a(v^) is a 
monotonically increasing function. 

A necessary condition for the lower and upper limits (4.20) and (4.33) to be compat- 
ible, therefore, is that 

In IR scenarios, however, the positive sign will apply in the last term of the right- 
hand side of Eq. (4.30). Hence, assuming > and neglecting the term proportional 
to ^^l/ /nl yields another upper limit on the tensor-scalar ratio: 

4(1-^2.) 

r, < \ . 4.36 

Combining conditions (4.34) and (4.36) therefore leads to a constraint on fsi'^*) for 
the generalised BM bound to be satisfied in IR inflation: 



fsM^^ {AK? V-sg 

To summarise, for the more general kinetic function in Eq. (4.22), the parameters 
/a and /b must satisfy Eqs. (4.28) and (4.34), where we have restricted our interest 
to the era of observable inflation. For UV models the lower and upper bounds on the 
tensor-scalar ratio will be compatible if Eq. (4.35) is satisfied. For IR models two upper 
bounds on r have been found, which when combined constrain fs as in Eq. (4.37). This 
constraint will prove useful in Section 5.4.1. 



4.4.2. The New Upper Bound for General Models 

The newly derived upper bound on r for DBI models, Eq. (4.11), arises because the warp 
factor in standard DBI models completely specifies the kinetic energy of the infiaton 
field. Deriving a corresponding bound for the generalised model (4.22) would be more 
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involved, since the CMB normalisation (4.28) only directly constrains the function 
/^((/j) and this may not necessarily depend on the warp factor. Instead, the constraint 
(4.9), given by 

can be combined with Eq. (4.32) to derive a limit on the tensor-scalar ratio in terms of 
the warp factor and P{X, ip). We find that 



<^7::T7^U7^ , (4-38) 



1/3 4/3 

This bound is valid for any P{X, ip) in the warped throat including the generalised 
DEI function given in Eq. (4.22). For a specific model where the warp factor and the 
functions fi{^) are determined by particle physics considerations, condition (4.39) may 
be interpreted as a bound that relates the tensor modes directly to the value of the 
inflaton field during observable inflation. This constraint provides a consistency check 
that any given model must satisfy irrespective of the form of the inflaton potential. 

It is worthwhile to compare Eq. (4.39) with the BM bound for the full evolution given 
in Eq. (4.31). To evaluate the BM bound requires knowledge of and over the 
whole of the inflationary era. In contrast, using Eq. (4.39) only requires values during 
observable inflation. However CgPx = /a/b/2 must be assumed to be slowly varying 
for Eq. (4.39) to be valid. As discussed in Section 3.4 this is a reasonable assumption 
for models in which Px is larger than 0{e'\j) during observable inflation. 

As the two bounds provide upper limits on r their relative strength can be compared. 
The bound (4.39) is stronger than the full throat BM bound if 

KI'N < 20 {yo\{X,)g,f' . (4.40) 

For typical field-theoretic values Vol(X5) c:^ 0{tx^), ~ O.lMpL and — 10~^, this 
implies that 

hf^N <?,m^-^^. (4.41) 
If the more conservative approach outlined above is taken, the BM bound for observable 
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inflation, Eq. (4.33), is weaker than Eq. (4.39) when 

h^J^N < 300 . (4.42) 

To summarise this section, new bounds have been derived which generalise those 
described in Section 4.2 to the case of the kinetic function in Eq. (4.22). The expressions 
(4.31), (4.33) and (4.39) imply that the bounds on r could be relaxed if 2csPx = 
/a/s ^ 1 on observable scales. It is therefore important to develop string-inspired 
models where this condition arises naturally. We will explore this possibility in the 
next chapter. 
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Figure 4.2.: A schematic of recent DBI inspired models 

We have found that the standard DBI model appears to be in conflict with observa- 
tions. Many attempts have since been made to modify the original scenario in order 
to evade the bounds derived above. These new models can be classified according to 
whether they involve single or multiple fields and single or multiple branes. Figure 4.2 
lists some of the models in each category. 

The most straightforward extensions of the DBI model are single field, single brane 
models. These have a single degree of freedom, as in the D3-brane model, but they 
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rely on other physical mechanisms to ease the bounds on r. A natural extension to 
the single D3-brane model is to consider a Dp-brane wrapped around a (j9 — 3)-cycle 
of the internal space. This leads to a change in the relationship between p and from 
that defined in Section 3.3 [21, 102, 202]. For example, Becker et al. [1 ] have proposed 
a model in which inflation is driven by a wrapped D5-brane. In this case, the range 
of allowed values for the inflaton becomes independent of the throat charge, A^, which 
weakens the upper bound on the tensor-scalar ratio to r < 0.04. (Strictly speaking 
there is a weak dependence on the charge since Ayj ~ N"^/^.) However, in arriving at 
this bound, it was assumed that backreaction effects of any fluxes in the throat were 
negligible. Kobayashi et al. [\ 02] considered both D5 and D7 wrapped brane models, 
but concluded that the former case required an excessively large background charge in 
order to relax the bounds on r. This requirement is highly constraining, but is still not 
as restrictive as the value of the charge required by the single brane scenario, which 
effectively rules this model out. Thus, a wrapped brane conflguration is preferable to 
the single D3-brane model, but the parameter space of the former is still severely limited 
by the WMAP5 observations [ ]. 

Another interesting proposal is warped Wilson line DBI. In this scenario, moduli flelds 
associated with Wilson lines play the role of the inflaton [12]. This scenario is T-dual to 
the standard DBI model with non-parallel branes. In general, the model describes the 
physics of a single brane with multiple position flelds and multiple Wilson line flelds. In 
Ref. [12], observational predictions were derived for the case when the brane position 
is flxed and only one Wilson line degree of freedom is used. This implementation 
is therefore a single brane, single fleld model. By following the method outlined in 
Section 4.2 for this single fleld model, a lower bound on r was derived, instead of the 
upper bound (4.12) [I 2]. The lower bound (4.20) remains valid for this scenario. There 
are, therefore, two lower bounds on r and the inconsistency of the standard DBI model 
is not replicated. 

Changing the physical setting can also allow larger fleld ranges, which in turn can 
relax the bounds on r. One such example is the case of a D4-brane in compactifled 
manifolds containing monodromies. The large fleld variations in this single brane, single 
fleld model lead to possibly observable tensor modes [I'^fi]. Although formulated in Type 
IIA string theory, the monodromy scenario has a simple inflationary interpretation as 
a large fleld, slow roll model with a potential V{lp) oc y)^/^. 

The tensor-scalar ratio and other observable quantities are signiflcantly altered if the 
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throat geometry is not of the AdS^ type, even in the case of the standard D3-brane 
model [ ■]. In Ref. [■''>■], a one parameter family of solutions was found, which interpo- 
lates between the Klebanov-Strassler (KS) ['J/] and Maldacena-Nuiiez [ ] throats. As 
the throat geometry moves away from KS, more non-Gaussianity is produced whereas 
the tensor-scalar ratio is reduced. The choice of throat geometry, therefore, could affect 
the bounds on r and must be considered when models are compared. 

The next class of models that can be investigated are the single brane, multi-field 
configurations. The warped throat is six-dimensional, so it is natural to consider cases 
where the D3-brane is not restricted to a radial trajectory. This was investigated 
in Refs. [" ] and [ ]. Increasing the degrees of freedom in this way introduces the 
possibility of entropy mode production. There are also changes in the predictions for 
the amount and type of non-Gaussianity produced and the constraints on r can be 
eased [10, 109-111, 141, 142, 161]. The bounds on r could also be affected if a non- 
negligible part of the curvature perturbation was produced by a curvaton field [129]. 
Curvaton fields arise generically in scenarios containing warped throats, particularly 
for propagation near the tip [101, 113]. 

Investigations have also been made into models with multiple branes, each of which 
has a single dynamical field. In Refs. [-i?] and [ ' ], no interactions between branes were 
considered, and the branes could conceivably propagate in different throats. The action 
for n decoupled branes in the relativistic limit is the sum of n copies of the DBI action 
(3.11). The power spectrum of curvature perturbations is enhanced by a factor of n^^"^ 
with respect to the single brane case. Consequently, the value of the tensor-scalar ratio 
will be reduced. 

We have not yet addressed models with multiple branes but only one effective degree 
of freedom. Multiple M5-branes in M-theory act with an effective single degree of 
freedom, but the Lyth bound is now significantly weakened. Large field ranges and 
an observable tensor signal are therefore possible [106]. Another proposal is that of n 
D3-branes which are coincident and propagating in a warped throat [28, 83, 194, 202]. 
The non-Abelian nature of the interactions between the branes differentiates this model 
from other multi-brane models and the model is also known as "Matrix Inflation" . In 
Chapter 5 we investigate this model in the relativistic limit for both large and small n, 
and show how the constraints derived in this chapter can be applied. 
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4.6. Discussion 

In this chapter, we have derived an upper hmit on the amphtude of the primordial 
gravitational wave spectrum generated during UV DBI inflation. We considered the 
maximal inflaton field variation that can occur during the observable stages of inflation 
and assumed only that the brane was propagating inside the throat during that epoch. 
The bound (4.12) is valid for an arbitrary inflaton potential and warp factor (modulo 
some weak caveats) and can be expressed entirely in terms of observable parameters, 
once the volume of the five-dimensional sub-manifold of the throat has been specified. 
The inferred upper limit on r is surprisingly strong. We find that the standard UV 
scenario predicts tensor perturbations that are undetectably small, at a level < 10^^. 

The current WMAP5 data favours models that generate a red spectral index, Ug < 
1, when both the gravitational waves and running in the scalar spectral index are 
negligible. For UV versions of the scenario, we have identified a corresponding lower 
limit on r which applies in this region of parameter space, > 0.1(1 — n^). It is 
clear that the standard scenario cannot satisfy both the upper and lower bounds on the 
tensor modes for the observationally favoured value 1 — ~ 0.03. 

The generality of our analysis implies that modifying either the inflaton potential or 
the form of the warp factor is unlikely to resolve this discrepancy. On the other hand, 
there are a number of possible ways of reconciling theory with observation. In general, 
either the upper or lower limit on r needs to be relaxed. Weakening the latter would 
require a violation of the slow roll conditions or a blue spectral index. A value of > 1 
is compatible with WMAP5 if the running of the spectral index is sufficiently negative, 
but is only marginally consistent if just the tensor modes are non-negligible. The 
upper limit on r can be weakened by reducing the volume of or by generalising the 
DBI action. Furthermore, it need not necessarily apply in IR versions of the scenario, 
although the BM bound will still hold in such cases. We considered a generalised version 
of the DBI action and identified a necessary condition on the form of such an action 
for the BM bound to be relaxed. 

In conclusion therefore, we have shown that primordial gravitational wave constraints 
combined with cosmological observations of the density perturbation spectrum act as 
a powerful discriminant of DBI inflationary models. They also serve as an important 
observational guide for identifying viable generalisations of the scenario. In Chapter 5 
we will explore one particular generalisation, the multi-coincident brane scenario intro- 



4.6: Discussion 73 



duced in Ref. [194]. 



5. Multi-Coincident Brane Inflation 



5.1. Introduction 

We have seen in Chapter 4 that the standard DBI inflationary model is severely con- 
strained by current observational data. The amplitude of tensor perturbations is 
bounded from above by r < 10^. When the brane is moving towards the tip of the 
throat, a complementary lower bound on r can be derived such that the two bounds 
are incompatible using current observational data. In this chapter we will explore how 
to evade, reconcile and weaken these bounds by considering a more general class of 
models that exhibit properties similar to the standard DBI scenario. 

In Section 5.2 we consider the special algebraic properties of the DBI action. We 
identify a general class of non-canonical inflationary models where the leading-order 
contribution to the non-Gaussianity of the curvature perturbation is determined entirely 
by the speed of sound of the inflaton fluctuations. In these models, the bounds on r 
can be relaxed if significant non-Gaussianities are generated. 

As reviewed in Section 4.5, many alternative ways to relax these bounds have been 
proposed, including models based upon multiple fields, the addition of angular momen- 
tum as another degree of freedom and using different throat geometries. However, in 
most cases the extra degrees of freedom introduced in these models do not solve the 
problem [1]. The bounds are relaxed only by a small fraction, and therefore these mod- 
els should still be regarded as unsatisfactory since they require an extreme amount of 
fine tuning in order to satisfy the observational constraints. 

One alternative possibility is to consider multiple brane configurations^. One scenario 
in which multiple branes are expected is after brane flux annihilation, in which branes 
travelling down the throat annihilate with the trapped flux, creating new branes [53, 
90, 194]. These are then attracted by other branes and fluxes in other throats and 

^In certain limits this approach is actually dual to considering wrapped branes [202]. 
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propagate toward the bulk. In Ref. [191] Thomas & Ward argue that it is unlikely that 
only a single brane is left after the flux annihilation process, due to the large amount of 
fine tuning necessary to achieve this. Instead it is more likely that a number of branes 
remain. 

In the case where n branes are localised initially at equal distances, I > Is, and 
subsequently follow the same trajectory, the effective theory is equivalent to that of n 
copies of the action for a single brane. A more general initial condition, particularly 
for branes created in the IR region of the throat [38, 53, 90], is that the branes should 
be separated over a range of scales, with a subset being coincident and the remainder 
being widely separated. 

In Section 5.3 we introduce the multi-brane model with n coincident branes described 
in Ref. [I') I]. We will consider two limiting cases of this model. The large n case is 
similar in form to the original DBI model, and we will show in Section 5.4 that it can be 
constrained using the formalism derived in Chapter 4. In contrast, the effective action 
in the relativistic limit for a finite number of branes is shown to be in the class of actions 
for which the bounds on r can be relaxed. In Section 5.5 we find that such models can in 
principle lead to a detectable gravitational wave background if the number of coincident 
branes is sufficiently small. 

5.2. Relaxing the Upper Bounds on the 
Tensor-Scalar Ratio 

In Chapter 3 we described the standard DBI scenario, in which the kinetic function 
P{(f,X) takes the form given in Eq. (3.12): 



where T{(p) = T^h^{ip) is the warped brane tension and V{(p) is the inflaton potential. 
The standard DBI scenario is algebraically special, in the sense that the kinetic function 
satisfies the constraints 



P(^, X) = -T(v9) v/l-2T-i(^)X + T(v9) - V{^) , 



(5.1) 




(5.2) 



5.2: Relaxing the Upper Bounds on the Tensor-Scalar Ratio 



76 



We saw in Section 4.4 that the bounds (4.33) and (4.39) on the tensor-scalar ratio could 
in principle be significantly relaxed in models where (cgP^x)* ^ 1. In view of the second 
relation in Eq. (5.2), it is of interest to begin by taking a phenomenological approach 
and to consider the more general class of models that satisfy the relation 

\-l = aA, (5.3) 

for some positive constant a. 

5.2.1. Approximate Solution 

A large non-Gaussian signature in the curvature perturbation is typically generated in 
models where the sound speed of fiuctuations is small. By substituting Eq. (5.3) into 
the definition of in Eq. (2.93), we can see that for this case 

/n1 (X A . (5.4) 

Recall that the definitions of and A in Eqs. (2.81) and (2.82) are in terms of P and 
its derivatives with respect to X. We can require that the magnitude of is large by 
considering scenarios in the limit where Cg is small and A is large, or equivalently: 

X^P,xxx > XP^xx > P,x . (5.5) 

Having taken this limit, the constraint (5.3) reduces to the third-order, non-linear, 
partial differential equation 

Pxx = ^P,xPxxx . (5.6) 
Changing the dependent variable to T = P^xx/P,x reduces Eq. (5.6) to 

aT^x = {Q-a)T\ (5.7) 

and it is straightforward to integrate Eq. (5.7) exactly. The remaining integrations can 
also be performed analytically and the general solution to Eq. (5.6) for a 7^ 6 is given 
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by2 

Picp,X) = -h{^) [1 - h{^)X]' - Mcp) , (5.8) 
where fi{(p) are arbitrary functions of the scalar field and 

(5.9) 

a — b 

The inequalities (5.5) are satisfied in the relativistic limit, where X ~ I//2, justifying 
their use. We consider the inflationary dynamics in the relativistic limit in what follows 
but for completeness we show in Appendix A.l that Eq. (5.3) can be solved analytically 
without this approximation. 

5.2.2. Consequences 

The standard DBI scenario is recovered from Eq. (5.8) for Z = 1/2, a = 2 (or s = — 1 in 
the exact case following a redefinition of fi). More generally, however, Eq. (5.8) implies 
that 

c.Px^^^{l-hXr-''/\ (5.10) 



.2 .1-/2^ 



2(1-0 



(5.11) 



when X ^ 1//2- Self-consistency therefore requires / < 1. Moreover we find from 
Eq. (2.93) that 

_ 5(59 - 55/) 

486 ' ' 
/eq ^ 5 /59 - 55/\ 

"" 972 i^3rJ- ^^-^^^ 

Hence substituting Eqs. (5.10) and (5.12) into the BM bound (4.33) and the bound 
(4.39) implies that 

..<-S.^(-#V'""" (5.14, 



^The special case a = 6 results in an exponential dependence of the kinetic function on X . However, 
we do not consider this model further, since it does not lead to a weakening of the gravitational 
wave constraints. 
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and 

respectively. 

We conclude, therefore, that the upper limit on the tensor-scalar ratio could be sig- 
nificantly relaxed if / < 1/2, since the non-linearity parameter is at present only weakly 
constrained at > —151. Although it is possible to phenomenologically construct 
a model which has a value of I in this range, it is clearly preferable to identify UV 
complete models that satisfy this requirement within a string theory context. Unfortu- 
nately this is quite difficult to achieve since the infiaton will either be associated with 
an open or closed string mode. The open strings are governed by relativistic actions of 
the DBI form, whilst closed strings arise from compactification of Einstein gravity and 
are typically put into a canonical form. However, there do exist classes of open string 
models which satisfy the above requirement, namely those associated with multiple 
coincident branes. 

More specifically, if the branes are spatially separated, the effective action is alge- 
braically equivalent to that of a single brane. It will therefore not satisfy the bound 
on l^. In the remainder of this chapter we will examine the case of n coincident branes 
as described in Ref. [191]. The large n limit of this configuration will also fall into the 
class of models with 1 = 1/2, which is equivalent to the single brane case. On the other 
hand, if it is assumed that n is finite, the special properties associated with the matrix 
degrees of freedom become important and this results in a kinetic function satisfying 
/ < 1/2. 

5.3. The Multi-Coincident Brane Model 

We have seen how the form of the kinetic function P can significantly change the 
strength of the bound (4.39) on the tensor-scalar ratio, depending on its explicit form. 
One model in which a suitable form for P is realised is the multiple coincident brane 
model as outlined by Thomas & Ward in Ref. [19 1]. In this model, the fiux annihilation 
process generates n coincident branes that are initially located at the bottom of a throat 

^In this discussion, we are ignoring the non-trivial backreaction of these branes on the background, 
and therefore one should be careful about the range of validity of the effective action. 
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region. The dynamics of this configuration is determined by the non-Abelian world- 
volume theory [150, 151]. This theory exhibits extra stringy degrees of freedom which 
arise due to the fuzzy nature of the geometry. 

In general the open string degrees of freedom for n coincident branes combine to 
fill out representations of U{n), as opposed to f/(l)" in the case of separated branes. 
This introduces a non-Abelian structure into the theory. In the single brane case, the 
fluctuations of the brane are characterised by induced scalar fields on the world-volume. 
However, for multiple branes these scalars must be promoted to matrix representations 
of some gauge group. Typically the transverse space of any given compactification will 
always admit an S0{3) isometry. Scalars can therefore be chosen to transform under 
representations of the algebra of 5*0(3) ~ SU{2) by making the identifications 

Lp' = Rma' i = 1,2,3, (5.16) 

where Rm is some scale with canonical mass dimension, and the a* are specified to be 
the irreducible generators satisfying the commutator 

[a\a^] = 2ieijka^ , (5.17) 

and the conditions 

-Tria'a^) = C6'^ = (n^ - 1)6'^ , (5.18) 
n 

where C is the quadratic Casimir of the gauge group. The irreducibility condition 
corresponds to the configuration being in the lowest energy state. It is therefore an 
additional fine-tuning of the initial conditions. 

The Myers prescription requires a symmetrised trace (denoted STr) to be made over 
the gauge group. This implies that the symmetric averaging must be taken over all the 
group dependence before taking the trace: 

STr{Ai ...A,) = -Tr{Ai ... A, + all permutations) . (5.19) 

S • 

For a large number of branes, n ^ 1, the symmetric trace can be approximated with a 
trace, which results in the usual DBI action multiplied by a potential term (as described 
in Refs. [90, 194]). However for finite n, the symmetrisation becomes more important 
and it is essential to have some means of performing this operation. A prescription for 
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treating tlie symmetric trace at finite n was proposed in Refs. [150] and [140], using 
highest weight methods and chord diagrams. The result is that the STr acts on different 
spin representations of SU{2) in the following manner: 

ra/2 

5Tr(a^a^)'' = 2(2g + l) J](22-l)29, neven, (5.20) 

i=l 

(n-l)/2 

STr{a'ay = 2{2q + 1) ^ {2if'', nodd. (5.21) 

In order for the solution to converge in this prescription, it is also necessary to modify 
the definition of the radius of the SU{2) sphere. In the large n limit, this is given by 

p2 = X^Rl^-Tr(a'a') = X^BLC , (5.22) 
n 

where A = 2nl1 = 27rm~^, whereas for finite n, it becomes 

= A^fiiLim,^^ ( ^|:^P' ) = A^fi^(n - 1)^ . (5.23) 

This converges to the large n result in the appropriate limit. This point is important, 
since the warp factor of the four- dimensional theory is typically of the form h = h{p). 
The next two sections will examine this coincident brane model in both the large and 
finite n limits. 



5.4. Coincident Brane Inflation with a Large 
Number of Branes 

Taking the limit of a large number of coincident branes significantly simplifies the 
non-Abelian action. The symmetrised trace can now be replaced with a normal trace 
operator and the expression for p takes the form in Eq. (5.22). For the case where a 
fuzzy two-sphere is embedded in a three-cycle in the X5 manifold, the kinetic structure 
of the action is given in the large n limit by [194] 



p = -nn 



(5.24) 
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wliere 

W{ip) = + C- 1/1-4 (5,25) 

defines the so-called 'fuzzy' potential, C = n^CT^ /m'^ is a model-dependent constant 
and (7 ^ in the large n limit. 

The kinetic term (5.24) is clearly of the same form as the single-brane DBI term, 
with Z = 1/2 in the scheme outlined in Section 5.2. We can therefore apply the analysis 
of Section 4.4 to investigate whether the bounds described in Chapter 4 can be relaxed. 
Comparison with Eq. (4.22), 

P = -f^(^)^l-f^(^)X~ fciv) , (5.26) 

implies that /a/b = 2nW and /b = 2/{T^h^). Hence, the new features of this model 
relative to the standard single-brane scenario are parametrised in terms of the fuzzy 
potential W{ip). This configuration is conjectured to be dual to a D5-brane which is 
wrapped around a two-cycle of the throat [14, 64, 193]. 



5.4.1. Bound on n during IR Propagation 

The regime ^ 1 is of interest for relaxing the gravitational wave constraints^. The 
generalised BM bound for IR models, with branes propagating towards the bulk, is 
given by Eq. (4.37): 

N 47r2 (l-n,)P|- ^ ' 



n 

As we know Jb, this may be expressed as a limit on the value of the warp factor h{(fj,) 
on CMB scales: 

.rJ^V<iteM,. ,5.8, 



MplJ v^=3a(AAC)2 

We now consider whether this limit can be satisfied for reasonable choices of param- 
eters when the warped compactification corresponds to an AdSs^ or KS throat, respec- 
tively. Recall that the warp factor for the AdS^ throat is given hj h = ip/{^/T^L). 
Condition (5.28) therefore reduces to a constraint on the value of the infiaton during 



■^Note that the case n ^ 1 and W ^ 1 will not significantly relax the BM bound, since we require 
n ^ N for backreaction effects to be negligible. 



5.4: Coincident Brane Inflation with a Large Number of Branes 



82 



observable inflation: 

However, non-perturbative string effects are expected to become important below a 
cutoff scale, (/^cut = ^cutV^-^, where /icut is the value of the warp factor at that scale. 
For consistency, therefore, one requires > fcut, so that 

which implies an upper limit on the D3-brane charge: 

, 64r^,.(i - ^ Mp, y 

Assuming the typical values rus ~ 10~^MpL, AA/"* ~ 4 and h^ut ~ 10~^ implies 

N < 1.76 X 10«(1 - n,)(-/^l)-i/2 , (5.32) 

and for 1 — < 0.05 and — /^l > 5 the inequality becomes 

AT < 4 X 10^ (5.33) 

For an Ac/S's throat, the fuzzy potential Ty is a constant, and the condition that 
W ^ 1 becomes 

Hence, combining inequalities (5.31) and (5.34) implies that 

C< ^'^:'ff,,li i^^V. ,5.35, 



-3/^l(AAr,)2 VoKXs) V^cutm, 
and specifying Qs ~ 10~^ and Vol(X5) ~ vr^ then yields the limit 

C < 2.25 X 10^(1 - ns){-f^lr^/^ < 5 x 10^ (5.36) 



or equivalent ly. 



n < 225. 



(5.37) 
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In deriving the action (5.24) the number of coincident branes was assumed to be large. 
However we have now found that for the case of branes propagating towards the bulk, 
the number of such branes is bounded from above. Furthermore, since /a/b — constant 
for the AdS^ throat, the stronger form of the inequality (4.33) may be used. The 
right hand side of inequality (5.35) would be reduced by a factor of (A/'es/^VV*)^ by 
substituting AA/"* — ?■ Afcs- This ratio could be as high as (60/4)^ ~ 200, leading to n 
being less than 15. In this case the assumption of large n would clearly be inconsistent 
and the model would be ruled out. 



5.4.2. Bound on D3 Charge at the Tip of the Throat 

Since the branes are initially located at the tip of the throat, another case of interest is 
the IR limit of the KS geometry, where the warp factor asymptotes to a constant value 
[71]: 

/itip = exp ■ (5.38) 



3Mg., 

In this case, the generalised BM bound (5.28) becomes 

8ttK , . f \ f 647r5(l - n^Vl , 



The radius of the three-sphere at the tip of the KS throat is of the order (g^My^'^ 
in string units and this must be large (and at the very least should exceed unity) 
for the supergravity approximation to be reliable. Substituting this requirement into 
expression (5.39) results in a necessary (but not sufficient) condition on the D3-brane 
charge for the generalised BM bound to be satisfied: 

Recall that a necessary condition for the backreaction of the branes to be negligible 
is iV ^ n and this implies that the exponential term in (5.40) will dominate unless 
the string coupling constant is extremely small. Hence, for the parameter estimations 
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quoted above, we deduce tlie lower limit 

iV- 121niV > -6.8 + 121n I — (^5 4^^ 




which becomes > 10^ for 1 — ~ 0.05 and — /^l > 5. 

In general, however, the K and M units of flux are not independent. F-theory com- 
pactification on Calabi-Yau four-folds provides a geometric way of parametrising type 
IIB string compactifications [52, 71, 77, 100, 178, 204]. Global tadpole cancellation con- 
strains the topology of the four-fold and this restricts the brane and flux conflgurations. 
When the KS system is embedded into F-theory, the constraint is given by [71] 

||=n + MK, (5.42) 

where x is the Euler characteristic of the four- fold. Hence, = MK < x/24 and 
together with condition (5.41), this implies that 

X > 2400 , (5.43) 

for > 10^. It is known that the Euler number for four- folds corresponding to hyper- 
surfaces in weighted projective spaces can be as high as x < 1,820,448 [100], so there 
are many compactiflcations that could in principle satisfy the generalised BM bound. 
On the other hand, the above limit on the Euler characteristic does allow us to gain 
some insight into the topology of the extra dimensions, since compactiflcations which 
result in a small Euler characteristic would be incompatible with the generalised BM 
bound. 



5.5. Coincident Brane Inflation with a Finite 
Number of Branes 

5.5.1. The Finite n Model 

The coincident brane model outlined in Section 5.3 takes a signiflcantly different form 
if, instead of assuming a large number of coincident branes, there are now only a small 
flnite number. The prescription for the symmetrised trace given in Eqs. (5.20) and 
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(5.21) must be used, where p is determined from Eq. (5.23). The resulting kinetic 
function for n coincident branes in the finite n limit is therefore given by 



P = -TsSTr (^h\p) (-^0'=r^(a^ar+^(^f ) + V{p) - h\p) 



(5.44) 



where 



Z^\^h-\p), Y^A\^Rth-\p), , (5.45) 



q J r(3/2-g)r(l + g)- 

Note that the second and third terms in Eq. (5.44) are singlets under the STr and 
therefore contribute terms proportional to n. The physics of these branes away from 
the large n limit is particularly interesting as discussed further in Refs. [19 1] and [202]. 

It was shown in Ref. [n.~] that the functional forms of the kinetic function, P, and 
corresponding energy density, E, for all the solutions with n > 2 can be derived re- 
cursively from the n = 2 solution. We will use the notation P„ and En to denote the 
non-singlet sector of the kinetic function and energy density for the n-brane solutions. 
The full pressure and energy densities are then given by P = P„ — nT^lV — h'^) and 
E = En + nT^CV — h'^), respectively. Using Eq. (5.20) and the expressions 

f:A'Ci^)=Vl + A, (5.46) 



A:=0 

oo 



k=0 

then implies that the terms P2 and E2 can be derived 




1 + 2r - (2 + 3Y)ZRl 



VT+Y^l - ZRl 
1 + 2Y - YZRj) 

Vl + Y{l - zRi)y^ 



(5.48) 
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The recursion relation described in Ref. [83] can then be written for odd n as 

/(n-l)/2 

^f^= E P2pkfZ,{2kfY]\-nT,{V-h' 

y k=i 

/(n-l)/2 

) = E ^2 [{2kfZ, {2kfY] I + nT,{V - h') , (5.49) 
y k=i 

and for even n as 

/ n/2 

^ = E ^2 [(2A; - ifZ, {2k - IfY] \ - nT,{V - h') 
yfc=i 

/ n/2 

) = E ^2 [(2A; - 1)2^, (2A; - 1)^^] | + nT^{V - h') . (5.50) 

If we let (5„ = 1 when n is even and 5^ = when n is odd, we can combine these two 
expressions [28]: 

^(n-l+<5„)/2 

Pn=[ E ^2 [{2k - 6nfZ, {2k - 6nYY] | - nT,{V - h' 

k=l 
'(n-l+<5„)/2 

En= [ E ^2 [(2A; - 5„)2Z, (2A; - 6nfY] ) + nT3(K - h^) . (5.51) 

fc=i 

In Ref. [28] it was shown that this recursive definition for P„ reproduces Eq. (5.24) in 
the large n limit. 

The backreaction of the multiple branes introduces corrections of the form n/N [83]. 
Ensuring that this ratio is small allows the continued use of the supergravity analysis. 
As we will see in the next section it will not be difficult to constrain this model when 

> ra. 



5.5.2. Bounds on the Tensor-Scalar Ratio for Finite n 



In the last section we introduced the multi-coincident brane model in the limit of a 
finite number of branes. In this section we will consider this model in the context of 
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tlie class of actions derived in Section 5.2, and show that current observational data 
can strongly constrain the ability of this model to produce an observable tensor signal. 
The last term appearing in the summation of P„ in Eq. (5.51) is 

pi-t = P2 [(n - ifZ, (n - ifY] . (5.52) 

This means that for all n, this term can be expressed in the form 

1 + 2(n - IfY - [2 + 3(n - lfY]{n - ifZRl 



>last 




(5.53) 



Inspection of Eqs. (5.48)-(5.51) implies that the relativistic limit is realised for any 
finite number of branes when {n—lYZR^ — )■ 1. In this case, the dominant contribution 
to the summations appearing in Eq. (5.51) will arise from the last term, Eq. (5.53). In 
the relativistic limit, therefore, the kinetic function appearing in the effective action 
simplifies to 



P = 2T, \^hyi + {n-lYY [l - '^'1 - nT, {V - h') 



where 



(5.54) 



" („ - D^yri {D • ^'-''^ 

2 

(n-l)W3^' ^^-^^^ 
and we have effectively imposed the relativistic condition 

X^^T3/i\ (5.57) 

in the numerator of Eq. (5.53). For the n = 2 and n = 3 cases, we have verified by direct 
calculation that when one calculates the speed of sound (2.81) and the non-linearity 
parameter (2.93) from the general expressions (5.48) and (5.49) and then imposes the 
relativistic limit (5.57), one arrives at the identical result by starting explicitly with 
Eq. (5.54). 

At this point we should consider the validity of the function in Eq. (5.54). The 
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recursive relation for P in (5.51) converges to tlie expression in Eq. (5.24) in tlie limit 
of large n [2n]. There must exist some value of n, beyond which P appears to resemble 
Eq. (5.24), rather than the approximate form proposed in Eq. (5.54). For a range of 
background solutions, numerical calculations suggest that the approximation is valid 
when n is less than 0(10) [S3]. 

As there are a large number of parameters in the theory, it is possible to find solutions 
where n ^ 10. However, a larger background flux would then be necessary, which would 
result in a situation where the conformal Calabi-Yau condition is no longer valid. In 
view of this, we focus on the sector of the theory where n < 10, which implies that the 
backreaction is under control and that the kinetic function is still of the required form. 

Eq. (5.54) is precisely of the form given by the general solution (5.8), where / = —1/2 
and^ 



M^) = -2T,hWl + {n- lyY, M^p) = . (5.59) 

We may, therefore, immediately conclude from Eq. (5.13) that /^l — —0.3/cl. More- 
over, since /3 ^ 0.9 in this scenario, Eqs. (5.10) and (5.12) reduce to 

CsPx ^ -1.3^1 + {n-iyYr^l . (5.60) 

We first consider the bound in Eq. (4.39). This applies at least for all UV scenarios. 
It follows after substitution of the relativistic limit (5.57) into the scalar perturbation 
amplitude, Eq. (2.86), that 

2 1 i/^ 1 

^^""^^wvrT(n-i)^r:a- ^^-^'^ 

Substituting the tensor-scalar ratio into Eq. (5.61) then results in a constraint on the 



^This is the case a 
fi{ip) becomes: 



18/5 or s = —1/3 in the analytic solution (A. 6) which after redefinition of the 



P 



h 8 - 4/2X1/3 _ 



/31 



3 • 



(5.58) 



This expression appears in a slightly different form to that in (5.53). However in deriving (5.53) 
we assumed the relativistic limit, which in turn imposes a non-trivial relation between X and 
Using this, and with a suitable redefinition of the functions, we can transform the above expression 
into the required form. 
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magnitude of the warp factor during observable inflation: 



Zl !:!^ (5 62) 

Eqs. (5.60) and (5.62) may now be substituted into the bound (4.39) to yield 

1100 [l + {n-in] 
< {AAf)t Vol(X5) ^ ^^^^^^ • ^^-^^^ 

It is clear that the parameter Y must be sufficiently large if the tensor perturbations 
are to be non- negligible. For the AdS^ x X5 throat, this parameter takes the constant 
value 

= {n - l)^Vol(X5) • ^^-^'^ 

As before, we choose natural field-theoretic values for the volume, Vol(X5) ~ tt^, and 
the string coupling, ^ 10^^, and further assume that {n — l^Y ^ 1. We again 
assume that the tensor-scalar ratio does not change significantly over the entire range 
of scales that are accessible to cosmological observation, which corresponds to AA/"* ~ 4. 
After substitution of the above values, therefore, the bound (5.63) simplifies to 

n<2.8xl0-i3^-^(ffl^ (5.65) 

As in Section 5.4, global tadpole cancellation constrains the magnitude of the back- 
ground charge N such that N < x/24. The maximal known value of the Euler number 
implies the upper limit of 

N < 75852 (5.66) 

for known solutions, although in principle higher values are possible. Imposing the 
WMAP5 bound > —151 in Eq. (5.65) and noting that n > 2 for consistency then 
implies an absolute upper limit on the tensor-scalar ratio: 

r, < 5 X 10"^ . (5.67) 

This limit is below the sensitivity of the Planck satellite (r > 0.02) [ ]. On the other 
hand, the projected sensitivity of future CMB polarisation experiments indicates that a 
background of primordial gravitational waves with > 10""^ should be observable [21, 
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188, 199]. In view of this, it is interesting to consider whether a detectable gravitational 
wave background could in principle be generated in this class of multi-brane inflationary 
models. We find from Eq. (5.65) that this would require 

n < 1 - 5.3 X IQ-^VNf^l < 1 - 0.014/^1 , (5.68) 

where the theoretic limit Eq. (5.66) for known compactifications has been imposed in 
the second inequality. We may deduce, therefore, that since we require n > 2 for 
consistency, a detectable tensor signal will require — /^l > 70. This implies that an 
observation of the tensors should also be accompanied by a sufficiently large — and 
detectable — non-Gaussianity. In other words, this class of models could be ruled out 
if tensors are observed in the absence of any non-Gaussianity. On the other hand, 
the current limit of > —151 implies that n < 3 is required for the tensors to be 
observable. Consequently, if tensor perturbations are detected, this would rule out all 
models with n > 4 or, alternatively, would require presently unknown configurations 
with N exceeding the bound (5.66). 

In the above analysis we assumed that the string coupling took the value Qs ^ 10~^. 
For the AdS^ x X5 throat, the bound (5.63) depends proportionally on gs and can 
therefore be weakened by allowing for larger values of the string coupling. For example, 
increasing this parameter by a factor of 4to Qs — 0.04 (so that it is still in the perturba- 
tive regime) relaxes the limit on the number of branes for the tensors to be detectable 
to < 5. Similarly, considering a smaller value for the volume of the Einstein manifold 
X5 will also weaken the upper limit. 

Let us re-iterate that this limit on n is well within the regime of validity for the theory, 
which we have argued is self-consistent for n < 10. Moreover, since the constraint (5.68) 
arises using the absolute maximal bound on the known Euler characteristics, it suggests 
that in realistic scenarios n will always be much smaller than this. Indeed, one could 
argue that only the n = 2 and n = 3 theories are likely to be valid over a large 
distribution of the fiux landscape. 

We must also ensure that our approximation (n — l)^y ^ 1 is valid for consistency. 
For the parameter values we have chosen this requires that QsN ^ (n — 1)^ and this is 
satisfied if the condition (5.68) holds. Note also that we require N ^ n for the super- 
gravity approximation to be under control and for backreaction effects to be negligible. 
This is also satisfied when Eq. (5.68) holds. 
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Finally, it should be emphasised that the derivation of the bound in Eq. (4.39) under- 
estimates the Planck mass by assuming that the volume of the throat is much smaller 
than the volume of the compactified Calabi-Yau three-fold. It is likely, therefore, that 
the actual constraint on r would be much stronger. Consequently, although the bound 
(5.68) does marginally allow for detectable tensors if n is sufficiently small, in practice 
this constraint would be further tightened by a more complete calculation. Nonethe- 
less, our analysis does not necessarily rule out these models as viable candidates for 
inflation. Rather, it suggests that it will be difficult to construct a working model that 
results in a detectable tensor signal. 

5.6. Discussion 

The relativistic DBI brane scenario represents an attractive, string-inspired realisation 
of the inflationary scenario. In Chapter 4 we showed that cosmological data has placed 
very strong constraints on the simplest UV models based on a single D3-brane. The 
strength of these constraints follows from field-theoretic upper limits on the tensor- 
scalar ratio, r, which in turn arise because the effective DBI action satisfies special 
algebraic properties. This provides motivation for considering generalisations of the 
scenario, in particular to multi-brane configurations. 

In this chapter we have identified a phenomenological class of effective actions for 
which the constraints on r are relaxed, if significant (and detectable) non-Gaussian 
curvature perturbations are generated during infiation. We have provided approximate 
and exact derivations of this class of models which coincide in the relativistic limit. 
It would be interesting to investigate whether the effective action (5.8) with values of 
/ 7^ —1/2 arises in string-inspired settings or elsewhere. 

In Section 5.3 we introduced the coincident n-brane model of Thomas & Ward [194]. 
We examined the predictions of this model in two limits, arbitrarily large n and small 
finite n. The large n model is similar to the single brane case. Using the results of 
Section 4.4, we showed that it is strongly constrained by current observations. 

The finite n model is of more theoretical interest as it exhibits the non-Abelian nature 
of the scenario. In Ref. [83] a recursive approach was derived to calculate the pressure 
and energy densities for n > 2 models using the n = 2 results. In the relativistic limit, 
these finite n models are included in the class of actions derived in Section 5.2 which 
relax the bounds on r. The backreaction of these models can also be kept well under 
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control. 

We proceeded to consider the question of whether the upper hmits on r could be 
relaxed to such an extent that a background of primordial gravitational waves might be 
detectable in future CMB experiments. The vast majority of string-inspired inflationary 
models that have been proposed to date generate an unobservable tensor background. 
We found that a detectable signal is possible, in principle, for typical string-theoretic 
parameter values if the number of coincident branes is either 2 or 3. This is consistent 
with known F-theory configurations and current WMAP5 limits on the non-Gaussianity. 
Furthermore, we found that the level of non-Gaussianity must satisfy — /^l ^ 70 if 
such configurations are to generate a detectable tensor signal. This is well within the 
projected sensitivity of the Planck satellite ['^^<^]. 

Our analysis invoked an AdS^ x X5 warped throat geometry. However, we made 
no assumptions regarding the form of the infiaton potential, other than imposing the 
implicit requirement that the universe underwent a phase of quasi-exponential expan- 
sion. In this sense, therefore, we have yet to explicitly establish that these inflationary 
models will be able to generate a measurable tensor signal. Nonetheless, since such a 
detection would provide a unique observational window into high energy physics, our 
results provide strong motivation for considering the cosmological consequences of these 
multi-brane configurations when specific choices for the infiaton potential are made. In 
particular, it would be interesting to employ the techniques developed in to identify the 
regions of parameter space that are consistent with current cosmological observations. 

In Part I of this thesis we have concentrated on using analytic techniques to constrain 
string-inspired inflationary models. In Part II, numerical techniques will be developed 
with the goal of constraining inflationary models using second-order perturbation the- 
ory. 
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6. Cosmological Perturbations 



6.1. Introduction 

Cosmological perturbation theory is an essential tool for the analysis of cosmological 
models. It will become more so as the quantity and quality of observational data 
continues to improve. With the recent launch of the Planck satellite, the WMAP 
mission reaching its eighth year, and a host of other new experiments, we will have 
access to more information about the early universe than ever before [104, 158]. 

To distinguish between theoretical models, it is necessary to go beyond the standard 
statistical analyses that have been so successful in recent years. As a result, much 
interest has been focused on non-Gaussianity as a new tool to classify and test models 
of the early universe. Perturbation theory beyond first order will be required to make 
the best possible use of the data. In Chapter 2 cosmological perturbations at first 
order were introduced. In Part II of this thesis, we outline an important step in the 
understanding of perturbation theory beyond first order, demonstrating that second 
order perturbations are readily amenable to numerical calculation, even on small and 
intermediate scales inside the horizon. 

Inflationary model building has for the past few years focused on meeting the re- 
quirements of first order perturbation theory, namely that the power spectra of scalar 
and tensor perturbations, as defined in Eqs. (2.62) and (2.67), should match those 
observed in the CMB. Inflationary models are classified and tested according to their 
predictions for the scalar power spectrum, scalar spectral index and tensor-scalar ratio. 
An important observable that arises at second order is the non-Gaussianity parameter 
/nl- As described in Section 2.6, this parameter is not yet well constrained by observa- 
tional data in comparison with the other quantities. In Part I, however, it was shown 
that /nl can already be used to rule out models with particularly strong non-Gaussian 
signatures. 

There are two main approaches to studying non-Gaussianity and higher order effects. 
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The first uses non-linear theory and a gradient expansion in various forms, either ex- 
plicitly, e.g. Refs. [164, 168], or through the AAf formalism, e.g. Refs. [112, 126, 128, 
171, 172, 191, 192]. 

However, a gradient expansion approach is restricted and can only be applied on 
scales much larger than the particle horizon. The second approach uses cosmological 
perturbation theory developed by Bardeen [15] and extends it to second order, e.g. 
Refs. [2, 16, 17, 27, -32, 60, 62, 127, 130, 1-35, 146, 152, 15-3, 17-3, 196]^ This approach 
works on all scales, but can be more complex in comparison to the AM formalism. The 
two methods lead to identical results on large scales [132]. We will follow the Bardeen 
approach here. 

In Section 2.6 the first order perturbations of the inflaton field were taken to be purely 
Gaussian. It is therefore necessary to go to second order if we are to understand and 
estimate the non-Gaussian contribution of any inflationary model (for a recent review 
see Ref. [136]). Deriving the equations of motion is not trivial at second order and 
only recently was the Klein-Gordon equation for scalar fields derived in closed form, 
taking into account the metric backreaction [133]. This allows for the first time a direct 
and complete computation of the second order perturbation, in contrast with previous 
attempts which have focused only on certain terms in the expression, for example 
Ref. [63]. 

In this chapter the equations of motion for first and second order field perturbations 
are described. These form the basis of the numerical calculation undertaken in Chap- 
ters 7 and 8. Chapter 7 describes the numerical implementation of the calculation, 
including the initial conditions used and the computational requirements. We outline 
the numerical steps taken in the system and examine the current constraints on the cal- 
culation. The calculation is based on the slow roll version of the second order equation, 
but solves the full non-slow roll equations for the background and first order systems. 
We present the results of the calculation in Chapter 8, including a comparison of the 
second order scalar field perturbation calculated for specific inflationary potentials. 

The models tested in this calculation are single field models with a canonical action. 
Significant second order corrections are expected only in models with a non-canonical 
action or multiple fields, or when slow roll is violated. Numerical simulations will be 
particularly useful in analysing models with these characteristics. Section 8.3 discusses 
planned future work to extend the current numerical system to deal with models beyond 

^For an extensive list of references and a recent review on these issues see Ref. [130]. 
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the standard single field, slow roll inflation. 

In Section 6.2 of this chapter the Klein-Gordon equations for flrst and second or- 
der perturbations are introduced. These will be the central governing equations of the 
numerical calculation. They are initially written in terms of the metric perturbations 
and then described in closed form, i.e., in terms of the fleld perturbations alone. In 
Section 6.2.3, the second order equation is written in a slow roll approximation. Sec- 
tion 6.3 describes the observable quantities that can be calculated from second order 
fleld perturbations. In Section 6.4 we discuss the results of this chapter. 

6.2. Perturbation Equations 

In this section we will briefly review the derivation of the flrst and second order perturba- 
tion equations in the uniform curvature gauge and describe the slow roll approximation 
that will be used. There are many reviews on the subject of cosmological perturbation 
theory, and here we will follow Ref. [1 >!']. The full closed Klein-Gordon equation for 
second order perturbations was recently derived in Ref. [133] and the results of that 
paper will be outlined below. 

6.2.1. Second Order Perturbations 

In Section 2.3 cosmological perturbations of a single scalar fleld were introduced at flrst 
order. The methods adopted in that section can be extended at second order to flnd 
gauge invariant quantities. Recall that scalar quantities such as the inflaton fleld, 93, 
can be split into an homogeneous background, ipo, and inhomogeneous perturbations. 
Up to second order ip becomes 

ifi{ri, x") = Mv) + SMV, x') + ^Sifi2iv, x') . (6.1) 

The metric tensor gf^j, must also be perturbed up to second order. In Eq. (2.33) the 
vector and tensor perturbations were included at flrst order. Here we consider only the 
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scalar metric perturbations [136]: 



-a" (1 + 201 + 02) 



9oo 



goi = a { Bi + -B2 



Qij = a' [(1 - 2V^i - ^2) + 2Ei,i,- + E2. 



(6.2) 



where 5ij is the flat background metric, 0i and 02 are the lapse functions at flrst and 
second order, ipi and -02 are the curvature perturbations, and i?2, Ei and i?2 are the 
scalar perturbations describing the shear. In addition to the flrst order transformation 
vector (2.36), there is now a second order transformation vector 



^2^ = (a2,/32;) 



(6.3) 



where the spatial vector part of the transformation has been ignored. 

As before, we can write down how a second order scalar quantity (such as 8(^2) will 
be transformed [132]: 



(6.4) 



where a tilde (~) denotes a transformed quantity. The metric curvature perturbation 
transformation at flrst order is straightforward, 1^1 = 1^1 — Hai, but at second order it 
becomes more complicated [io6]: 



V.2 = V^2 - 'Ha2 - jXl + '-V-^r^i, , 



1, 
4 



(6.5) 



where Xij is given by 



Xij =2 



a 



+ 4 



+ 2 (5iiaij + 5yai,,) + AUai (^uj + - 2ai,iaij + 2^1^,^^/ 



(6.6) 
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and Bii and Cuj were defined in Eq. (2.34). 

We will work in the uniform curvature gauge where spatial 3-hypersurfaces are flat. 
This implies that 

^1 = ^2 = ^1 = ^2 = 0. (6.7) 

These relations can be used at first and then at second order to define gauge invariant 
variables. It follows from Section 2.3 that the first order transformation variables in the 
flat gauge satisfy ai = ipi/H and /3i = — -Ei. At second order, for the transformation of 
scalar quantities, as in Eq. (6.4), we require only a2- This is found by using Eq. (6.5) 
to have the form 

where the flrst order gauge variables have been substituted into Xij. 

The Sasaki- Mukhanov variable, i.e., the fleld perturbation on uniform curvature hy- 
persurfaces [145, 170], was given at flrst order in Eq. (2.45) as 

6^1 = 6^1 + . (6.9) 
At second order the Sasaki-Mukhanov variable becomes more complicated [132, 135]: 



n 



where ^igat = —{Ei^i + Fn). From now on we will drop the tildes and only refer to 
variables calculated in the flat gauge. 

The potential of the scalar fleld can also be separated into homogeneous and per- 
turbed sectors: 

V{^) = Vo + 6Vi + ^6V2, (6.11) 

5Vi = V^5ipi, (6.12) 
SV2 = V^^5ipl + V^Sip2. (6.13) 

Finally, the Klein- Gordon equations describe the evolution of the scalar fleld and 
are found by requiring the perturbed energy-momentum tensor T^,^ to obey the energy 
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conservation equation V f^T^^ = (see for example Ref. [132]). For the background 
field, (fo, the Klein-Gordon equation is 



Sip'; + 2-HVi + 2aVv0i - V^Sipi - ifi'oV^Bi - ^[,0; + a^V^^Sifii = , (6.15) 



where all the variables are now in the flat gauge. 

In order to write the Klein-Gordon equations in closed form, the Einstein field equa- 
tions (2.5) are also required at first and second order. These are not reproduced here, 
but are presented for example in Section II B of Ref. [133]. 

6.2.2. Fourier Transform 

In general, the dynamics of the scalar field becomes clearer in Fourier space. However, 
terms at second order of the form {5(pi{x'^)) require the use of the convolution theorem 
(see for example Ref. [-iJl]). Following Refs. [ i i] and [13.j] we will write 5(p{k'^) for the 
Fourier component of 5ip{x^) such that 



ipl + 2H<^o + aV^ = . 



(6.14) 



The first order equation is 



and the second order version is given by 



(6.16) 




(6.17) 
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where is the comoving wavenumber. In Fourier space, the closed form of the first 
order Klein-Gordon equation then transforms into 



2 



n 



n 



6^i{k') = 0. (6.18) 



The second order equation requires a careful consideration of terms that are quadratic 
in the first order perturbation. In particular, we require convolutions of the form 



f{x')g{x' 



1 



(2vr) 



d'qd'p6%k'-p'-q')fip')g{q' 



(6.19) 



For convenience we will group together those terms with gradients of 5ipi{x^) and denote 
them by F. The full closed form, second order Klein-Gordon equation in Fourier space 
is then given by [133] 



5ip'^{k') + 2HV2(F) + k^5ip2ik 

,2 



+ a" 



+ 



(27r)3 
167rG 
U 

871G 



n 



n 



5ip2{k^ 



n 



[Q6^\{p')6vM) + ^',o?V^^5^,{p')5v,{q')\ 



(/^'o [2a^V^^'^5^^{p')5^i{q') + ^'^Q5^i{p')5^i{q')\ 



+ F(5^i(F),Vi(fcO) = 0, 



SttG 



5^i{p')5^i{q' 



(6.20) 



where we have defined the parameter Q = a^^STcGVoipQ/T-L + V^p) for convenience. The 
F term contains gradients of S(pi in real space and therefore the convolution integrals 
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include additional factors of k and q. The form of F is given by [133] 



n J q^ 

\ n J n 



n 



A;2 



Q / 47rG\ ^ piq^p^q"" + p^g ^ 

n \ n J k^q" 



+ 



+ 



2 / 



A;2 



\ n J n 



PiqPmq 
p2q2 



+ 



+ 



n 



SvrG 



piq^+p^ 2 , , i 

— p — Q ^vivp 



q^ + Piq^ 
k^ 



47rG'\ yki ( „p^Pj 



n J k^ 



2^ {Q6Mp') + ^'o^^'iip')) QSMq" 



P 



(6.21) 



6.2.3. Slow Roll Approximation 

In order to establish the viability of a numerical calculation of the evolution of second 
order perturbations from the Klein-Gordon equation, Chapters 7 and 8 will be limited 
to the framework of the slow roll approximation. This involves taking 



Mi: 

2a2 



(6.22) 
(6.23) 



such that Q ~ and ~ (87rG/3)a^Vo. In Chapter 2 the slow roll parameter eh 
was defined in Eq. (2.20). In this chapter and the rest of Part II, a different slow roll 
parameter will be used, denoted by Eh and defined in Refs. [] •) '] and [ ]. This new 
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parameter is the square-root of Eh and is given by 

£H = v^^ = yi^. (6.24) 

The second slow roll parameter is still rju = Sh — ^'h/'^'^^h- Following Ref. [I -)'!] we 
will implement the slow roll approximation by keeping terms up to and including 0{e^fj) 
and terms which are OIehVh)- Within this approximation the second order equation 
(6.20) simplifies dramatically, and with the F term included it reduces to 

dip'iik') + 2n6ip'^ik') + kH^2{k') + (aV^^ - 247rG((^'o)2) 5<^2(fc') 

^ /" J3„ j3„ r3o i i J\\ „2 fjr , SttG . 



167rG 2 I t r c- / ?\ f- / i 

+ - ^) M(P-)^^;(,')' 

= 0. (6.25) 

The numerical simulation described in Chapter 7 will solve the slow roll version of the 
second order equation given above, Eq. (6.25), together with the complete first order 
equation (6.18) and background equation (6.14). 



6.3. Observable Quantities 

Cosmological perturbations at second order are becoming increasingly important now 
that statistical quantities beyond the power spectrum and spectral index are being 
investigated. Observations, however, do not tell us anything about the infiaton field 
directly. In this section the second order perturbations described above will be related 
to observable quantities in order to demonstrate how a numerical calculation could be 
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employed in the near future to gain further insight into the nature of the field that 
drives inflation. 

The temperature fluctuations observed in the CMB can be directly related to the 
curvature perturbation TZ. In Section 2.3, TZ was defined at first order in terms of 
5ipi. When the second order contribution is included the total comoving curvature 
perturbation is defined as 

7^ = 7^l + ^7^2 . (6.26) 

The first order term is related to the inflaton perturbation in the flat gauge by TZi = 
TiSipi/ipQ. The second order part includes terms quadratic in 6ipi and so in Fourier 
space requires convolutions. We are interested in the value of TZ after horizon crossing 
for the calculation of and a determination of the non-Gaussianity produced during 
inflation. This allows us to neglect gradient terms in real space or terms proportional 
to k in Fourier space. In this limit the real space expression for TZ2 is [132] 

n,{v, = ^S^2 - 2-^ Vi^^i + M f - ^' - 2^'l • (6.27) 

^0 Wo) Wi)) \ "Po J 

Using the background evolution equation (6.14) and transforming to Fourier space 
implies that Eq. (6.27) can be written as 



^0 {^^r J 



27TT2'^¥'i(^'P')Vi(^,gO 

2^2^ + "^V^ + ^SnG)a'Vo ) 5^i(r/,p^)5^i(r^,g^) 



(6.28) 



Once the numerical calculation has been carried out at first and second order as de- 
scribed in Chapter 7 this quantity can be evaluated after horizon crossing. 

In Chapter 2 the non-Gaussianity parameter defined in terms of TZ in 
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Eq. (2.89). Writing Eq. (2.89) in Fourier space using Eq. (6.26) implies that 



7^(F) = 7^l(F) + ^/^-^^ / d.q'nMl'R-^ik^-q' 



1 



(27r)= 



(6.29) 



where () denotes the expectation value. A good approximation of the local non- 
Gaussianity produced is then given by 



5, 
6' 



/nl = ^T^2{k' 



(27r)= 



dg3 7^l(g*)7^l(F -g^ 



(2vr) 



dq^ni{q')ni{k'-q') 



. (6.30) 



Calculating 6(p2 and 7^2 therefore provides direct insight into the behaviour and pro- 
duction of the non-Gaussianity parameter f^^. 

To go beyond the local shape of the non-Gaussianity it is necessary to calculate 
the full bispectrum of the perturbations. In practice the bispectrum of the curvature 
perturbation on uniform density hypersurfaces, (, is used in setting observational limits. 
At first order this is simply related to the comoving curvature perturbation by (i = 
—TZi. At second order the relationship is more complicated. For large scales outside 
the horizon, ^2 can be related to the field perturbations in real space using [132] 



n 



4-3 



n 



(6.31) 



In Fourier space this again introduces a convolution integral of the first order pertur- 
bations. 

The bispectrum of ( is given by 



(C(ki)C(k2)C(k3)) = (27r)35(ki + k2 + k,)B{k,, k^, k,) 



(6.32) 



where translation invariance introduces the delta function. The k dependence of the 
bispectrum is usually separated from an overall amplitude factor and considered as a 
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shape function F{ki, k2, k^). The bispectrum is then of the form [13, 118] 

(C(ki)C(k2)C(k3)) = A(27r)35(ki + k2 + k,)Fik,, k^, k,) , (6.33) 

and for a particular shape function F the best estimator for A when the non-Gaussianity 
is small is given by [^J] 



. ^ Eic. C(ki)C(k2)C(k3)F(fci, fc3)/««; 

Ek,F(A:i,A:2,A:3)V««) ' ^ ^ 



where 0"^^ is the variance of the mode and the sums run over all the triangles in k space 
subject. If ks is chosen to be the longest side of the triangle then the triangle inequahty 
enforces 

k3<ki + k2. (6.35) 

Eq. (6.34) provides a blueprint for how to evaluate the bispectrum in terms of a partic- 
ular given shape. To compare a primordial bispectrum with the observed temperature 
bispectrum from the CMB it is necessary to construct the spherical harmonics of the 
bispectrum and use transfer functions to relate the primordial values with the observed 
values. We have not carried out this procedure in this thesis. However, one of the goals 
of our future work is to undertake such a comparison of the numerically generated 
bispectrum with the observed quantity. 

As mentioned above, the shape most often used in comparisons with observations is 
the local shape given by the ansatz in Eq. (2.89). The expression for Fiocai is [13, 105] 

Fiocai(A:i, k2, k,) = 2iV/^°£ (i + + i) ' (6-36) 

where the spectrum has been taken to be scale invariant and is a normalisation 
factor. 

This is not the only shape that has been considered and, as we have seen, non- 
canonical inflationary actions generate a bispectrum which is peaked when the magni- 
tudes of the k modes are approximately equal. The form of F for the equilateral case 
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is [105] 



1,3 

2 



3 J 



^2 



+ 



^1^2 ^3 



+ 5 perms 



(6.37) 



A third form has been proposed which is nearly orthogonal to the other two shapes 
[177]. These shapes all have the property that they are separable functions of each 
ki or can be constructed from these separable functions. This property eases analytic 
calculations but clearly does not hold for generic shapes. There has been a proposal to 
define the shape of the bispectrum in terms of a set of basis vector shapes [61, 118]. This 
would remove the need for only separable shapes to be considered and would allow for a 
straightforward analysis of the bispectrum from its primordial value up to the observed 
bispectrum in the CMB. 

We have seen that the second order scalar perturbation is not the direct observable 
quantity of interest. The bispectrum of the curvature perturbations, which contains a 
contribution from the second order non-linear part, can be compared with observations 
either by use of various shape functions or through a full analysis with transfer of the 
primordial values. A future goal of our work is to compare the bispectrum obtained 
numerically with that from observations. 



6.4. Discussion 

In this chapter, the equations of motion for a single scalar field up to second order in 
cosmological perturbations have been introduced. The second order gauge transforma- 
tion has been discussed and the transformation components determined for the uniform 
curvature gauge. In Chapter 2 first order classical perturbations were quantised in the 
Minkowski spacetime limit and normalised using the Wronskian condition in Eq. (2.54). 
This constraint also fixes the quantisation for other orders of the perturbation, including 
6<^2 [ ]■ 

The perturbation equations are better understood in Fourier space, although the cost 
of adopting this approach is the need to employ convolution integrals of the first order 
perturbations. When written in Fourier space, the second order Klein-Gordon equation 
can be described entirely in terms of the field perturbations and background quantities. 
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Eq. (6.20), first derived in Ref. [133], is valid on all scales inside and outside the 
horizon. When a particular slow roll approximation is made this equation simplifies 
to that found in Eq. (6.25). This slow roll version of the equation will be the central 
governing equation of the numerical calculation described in the next chapter. 

Finding numerical solutions of Eq. (6.25) is the first step towards solving the full equa- 
tion (6.20) for a single field and ultimately the multi-field equation given in Ref. [133]. 
Understanding cosmological perturbations beyond linear order is critical if higher order 
statistical effects are to be accurately calculated. Section 6.3 outlined the connection 
between 5ip2 and observable quantities such as the comoving curvature perturbation and 
the non-Gaussianity of the perturbations. Going beyond single field, slow roll models, 
non- linear effects become more important. In Chapters 7 and 8 the first step is taken 
towards calculating higher order perturbations for these models. 



7. 



Numerical System and 
Implementation 



7.1 



Introduction 



Our goal in Part II of this thesis is to show that, just as at first order, a direct numerical 
calculation of the second order perturbations of a scalar field system is achievable. 
In this chapter the implementation of this system is outlined. The structure of the 
numerical system follows the work done at first order by Martin & Ringeval [137, 165] 
and previously by Salopek et al. [ ]. 

The most important difference between an analytic and numerical treatment of the 
equations presented in Chapter 6 is the requirement to specify a finite numerical range 
of a finite number of k modes to be calculated. The upper cutoff in k, which marks 
the smallest scale considered, is well motivated by the difficulty in observing primordial 
perturbations at very small scales. On the other hand, we also need to specify the largest 
scale (or smallest k) that we will consider. Analytically, this is often taken to be the 
size of the universe, with k = being the equivalent mode. One immediate problem 
with this specification, however, is that the Bunch-Davies vacuum initial conditions 
diverge. The standard approach to this problem is to implement a cutoff at large scales 
beyond which the amplitude of perturbations vanishes. This is a pragmatic approach, 
but recently there has been some evidence that a sharp cutoff similar to this could be 
responsible for the lack of power at large scales in the WMAP data [95, 124, 187, 189]. 

The main concern is that the k range covers most, if not all, of the modes ob- 
served to date in the CMB. The WMAP team rely for their main results in Ref. [104] 
on i multipoles in the range i G [3,1000], which corresponds approximately^ to k G 
[0.92 X 10^*^°, 3.1 X 10'^^] MpL = [3.5 x 10-^ 0.12] Mpc'^ We will consider three dif- 

^The approximate conversion for £ is £ ~ and a Megaparsec is given in Planck units as IMpc^^ ~ 



108 



7.2: Numerical Equations 



109 



ferent ranges of k modes when producing the results in Chapter 8, all of which contain 
the WMAP pivot scale ^wmap = 0.002Mpc-\ The choice of k range is flexible with 
the only numerical constraint being that the number of modes at second order should 
be equal to 2' + 1, where / is a positive integer. This enables faster integration using 
the Romberg method, as explained below. 

In Chapter 6 the Klein-Gordon equations of motion were described for the background 
field, and the first and second order field perturbations. These form the basis of the 
numerical calculation in this chapter. In Section 7.2, these equations are rewritten in a 
form more amenable to numerical work. The time variable is changed from conformal 
time, 7], to the number of e-foldings. A/". The convolution terms which are present in 
the Fourier space equations are then written in terms of spherical polar coordinates and 
split into smaller units. 

Four inflationary potentials were chosen in order to test the numerical calculation. 
These are defined in Section 7.2.1 and the steps taken to establish the values of the 
required parameters are outlined. The initial conditions used for the first order pertur- 
bations are the Bunch-Davies vacuum conditions, as specified in Section 2.3. At second 
order the perturbations are set to be identically zero at the beginning of the simulation, 
as explained in Section 7.2.2. This section also describes the method and timing of the 
initialisation of the variables. 

In Section 7.3 the numerical method is discussed. The calculation can be split into 
four stages, each of which is described in depth, along with the logistical constraints and 
software environment. The numerical code is tested in Section 7.4 by comparing the 
computed value with an analytic result where this is possible. The choice of parameters 
in the calculation is determined by the results of this comparison. In Section 7.5, the 
results of this chapter are summarised and discussed. 

7.2. Numerical Equations 

The Klein-Gordon equations in Chapter 6 are not appropriate for a numerical calcu- 
lation and in this section we rearrange them into a more suitable form. This involves 
a change of time coordinate and grouping of terms into smaller units for calculation. 
The second order slow roll equation (6.25) can be further simplified by performing the 



2.6247 X 10-"MpL. 
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p integral and changing to spherical polar coordinates q,6 and u, where q = |q|. The 
d^g integral then becomes 



q'^dq / sin 9d6 / du . 
Jo Jo 



(7.1) 



For each k mode equation we take the 6 = 0,u = axis in the direction of k\ so that 
the angle between /c* and g* is 6, and the scalar product is g^/c* = qk cos 6. The argument 
of each S(pi or 6(p'-^^ term depends on 6 through — g*| = ^/k"^ + q'^ — 2kq cos 6 and so 
must remain inside the 6 integral. There is no u dependence in 6(pi with this choice of 
axes, so the last integral is straightforward to evaluate. 

In the slow roll case there are only four different 6 dependent terms, here labelled 
A-V: 



A{k\q') 
B{k\q^) 
C{k\q') 
V{k\q') 



sm{e)5ipi{k' - q')de , 
cos{9)sm{9)Sifii{k' - q')d9 , 



sm{9)Sip[{k' - q')d9 



cos{9) sm{9)6^[{k' - q')d9 . 



(7.2) 



When written in terms of the variables defined in Eqs. (7.2), the slow roll equation 
(6.25) becomes: 

Sif'^ik') + 2H5^'^{k') + kH^2{k') + (aV^^ - 247rG(^'o)2) 5^2{k') + S{k') = , (7.3) 

where S{k'^) is the source term which will be determined before the second order system 
is run: 



S{k' 



dg <a^V^^^q^6^i{q')A{k\q') 



3a%^g2 + ^g^ + 2fcV) A{k\q' 



q^Cik\q^)+(^2-^^ kqV{k\q^ 



9 g2 

- + — 
2 A;2 



kq^B{k\q' 



6^i{q' 



(7.4) 
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The full set of equations which must be evolved are then Eq. (6.14) for the background, 
Eq. (6.18) for the first order perturbations and Eqs. (7.3) and (7.4) for the second order 
and source terms. 

A more appropriate time variable for the numerical simulation is the number of e- 
foldings (2.17). We employ 

M = log(a/ainit) (7.5) 

as our time variable instead of conformal time. This is measured from Oinit, the value of 
a at the beginning of the simulation. We will use a dagger^ to denote differentiation 
with respect to JV. Derivatives with respect to conformal time, 1], and coordinate 
time, t, are then given by 

d7] dr] dAf dJ\f' ^ ' 

d_^dTi(W_d_^ _d_ 

dt dt dfj dAf dAf' ^ ■ ^ 

respectively, where H = dlna/dt and = aH. 

If a is set to be unity at the present epoch, we can calculate a^^it once the background 
run is complete and the number of e-foldings of inflation has been determined. The 
value of a at the end of inflation, Oend, is calculated by connecting it with gq (see for 
example Eq. (3.19) in Ref. [U I] or Eq. (7) in Ref. [155]). The relation between them is 
given by 

= , I'-oj 

ao Ctreh ^eq O-o 

where Oreh and are the values of a at the end of reheating and matter-radiation 
equality. Using the relationship between energy densities and the scale factor relevant 
to the matter and radiation dominated eras, together with the Friedmann equation 
(2.10), we can write 

Oend \ ^ 1 / H^Yid \ , 1 , / -^^reh \ , ^ i / -^^eq \ , i / '-^cq 



where -ffend, -f^reh and i^eq are the values of H at the end of inflation, at the end of 
reheating and at matter-radiation equality respectively. The value of a^q is taken to be 
4.15 X 10-^{Q„^h^)-^ and Heq = 4.63 x IQ-^^fi^/i^MpL [55, 155]. The mean value for 



^This should not be confused with the use of ^ as Hermitian conjugate, which is confined to Section 2.3. 
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ilmh^ determined by WMAP5 + BAO + SN measurements is fi^/i^ = 0.1369 [104]. 
With these values and taking oq = 1, the scale factor at the end of inflation is given by 



Q^end 



-72 



end 



PL 



reh 



M, 



PL 



(7.10) 



In Chapter 8, it is assumed that reheating occurs instantaneously at the end of 
inflation such that i^reh = -f^end- This gives a^nd — 10~^^ and approximately 65 e- 
foldings from the end of inflation to the present. It also fixes the horizon crossing time 
of the WMAP pivot scale, fcwMAP = 0.002Mpc~\ to be about 60 e-foldings before the 
end of inflation. Because the Hubble parameter is not kept fixed during the numerical 
calculation the number of e-foldings between horizon crossing and the end of inflation 
will depend on the form of the inflationary potential and the evolution of H. 

The background and first order equations, written in terms of the new time variable 
M, are 

and 



0, 



(7.11) 



3 + ^1 6^lik^) + 



+ 



8nG 



Sifii{k') = 0, (7.12) 



respectively. The corresponding second order perturbation equation takes the form 



+ 



V. 



247rG((/^S)' ) bip2{k') + S{k') = , (7.13) 
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with the source term given by 



+ 



8nG 
(oH) 



l + Y^]kq'B{k\q^) 



■ 2' 
5^i{q') 



q^C{k\q^)+[2-^-]kqV{k\q^) 



(7.14) 



where 



C{k\q') = —C{k'-q' 



V{k\q^) = —V{k^-q^) 

an .iQ 



sm{9)6^\{k' - q')d9 , 
cos{e) sm{e)5ipl{k' - q')de 



(7.15) 



The arguments of 6ipi and 6ip\ in the A-T> terms require special consideration. To 
compute the integrals, 9 is sampled at 



N0 = 2^ + l 



(7.16) 



points in the range [0, tt] (for some / G Z"*" to allow Romberg integration) and the 
magnitude of k^ — is found using 



|F - q'\ = ^k^ + q^ -2kqcos{e). 



(7.17) 



While 5ipi{k^) = 6ipi{k), the value of — g*| is at most 2/cinax, where k,q E [/cmin, ^max]- 
This means that to calculate the source term for the k range described we require that 
6ipi and 6ip\ be known in the range [0,2A;inax]- In Section 7.3, we will show that this 
first order upper bound does not significantly affect performance. On the other hand, 
— g*| can also drop below the lower cutoff of calculated k modes. As discussed above 
a sharp cutoff will be implemented and 6y:>i{k) = used for the values below fcmin- 
When the spacing of the discrete k values, Ak, is approximately /cmin the cutoff affects 
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only the k = q modes and is only significant close to kmm- Section 7.4 describes how 
the accuracy is affected by changing Ak and other parameters. Without extrapolating 
outside our computed k range it appears to be very difficult to avoid taking 6{pi = 
for a small number of k values below fcmin- 

The value of — g*| will not in general coincide with the computed k values of Sifi. 
We use linear interpolation between the neighbouring k values to estimate 6ipi at these 
points. We leave to future work the implementation of a more accurate and numerically 
more intensive interpolation scheme. 

Throughout the discussion above we have not specified any particular inflationary 
potential, V. Indeed the numerical code can use any reasonable single field potential 
provided that it drives a period of inflationary expansion in the e-folding range being 
simulated. In the next section the four potentials which have been tested are outlined. 

7.2.1. Potentials and Parameters 



Potential Parameter Value 

m 6.32 X lO^'^MpL 

A 1.55 X 10-13 

10 

a 3.82 X lO-i^Mp^L 

mo 1.74 X lO-^MpL 

Table 7.1.: The parameter values for the four potentials, chosen so that 'P-^^(A;wmap) is 
in agreement with the WMAP5 value. 



1 2 2 



a(p3 



To demonstrate the numerical calculation four different single field, slow roll poten- 
tials were chosen. These are not intended to represent an exhaustive selection, but they 
do exhibit an interesting variety of behaviours. The potentials used are: 

1. V{!f) = ^m'^if'^. This is the original chaotic inflation model which is still in good 
agreement with the observational data [4]. 

2. V{ip) = |A(y9^. Although increasingly in tension with observations this is a stan- 
dard large field model. 

2 

3. V{ip) = aip^i. This potential with an unusual fractional index is the effective 



7.2: Numerical Equations 



115 



potential resulting from the monodromy inflation model of D4 branes, where 
observable tensor modes are possible [1, 

4. V{{p) = Uo + ^m^Lp'^. This is a contrived toy model which requires inflation to be 
terminated by hand. We will set inflation to end when ^ 8. By taking a value 
of t/o = 5 X 10^^°MpL a blue spectrum (n^ > 1) can then be obtained [104, 120]. 



xlO" 



xlO 




(a) Viip) = imV' 



xlO- 





15 10 

V3/7l/pL 



(b) Viip) = ix^^ 



XlO-'' 




14 12 

(d) Vi^) = C/o + im2^2 



(C) V{ip) = aip3 y^j , y^i - I 2 

Figure 7.1.: Plots of the four different potentials investigated. 



In Figures 7.1 and 7.2 the potentials are plotted over the course of their evolution. 
Throughout the rest of this chapter we will use the quadratic model to demonstrate 
the calculation unless otherwise stated. In Chapter 8 the results for each potential will 
be compared. 

For each of the chosen potentials there is one free parameter that needs to be deter- 
mined. We choose the parameters m, A, a and mg so that calculated for each model 



7.2: Numerical Equations 



116 
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Knd - A/" 

(b) The potentials in terms ol J\f for the last 70 e-foldings of inflation. 



Figure 7.2.: Two comparisons of the four potentials. 
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is in agreement with the WMAP5 value at the pivot scale /cwmap = 0.002Mpc^ ~ 
5.25 X 10~^°MpL. The dependence of 'P-^^(A;wmap) on each of the parameters can be 
seen in Figure 7.3. Requiring P-l^^(fcwMAp) = 2.457 x 10~^ gives the values shown in 
Table 7.1. Here we have chosen the lower of the two possible values of mo shown in 
Figure 7.3d. 




6.2 6.3 



6.4 6.5 



(a) Fixing m for V{(p) = ^m?(p'^. 
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(c) Fixing a for V{(p) = cr(/?3 . 
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(b) Fixing A for V{lp) = lX(p'^. 
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(d) Fixing m for V{(p) = Uq + \ra'^Lp'^ 



Figure 7.3.: Parameter values for the different potentials are chosen by requiring con- 
sistency with the WMAP5 normalisation of the first order power spectrum. 



7.2.2. Initial Conditions 

The background system requires initial conditions for yjQ and ip\. These initial condi- 
tions and the range of e-foldings to be simulated must be selected with the choice of 
potential in mind. Not only must the e-folding range include an inflationary period, but 
the k modes to be calculated at first and second order must begin inside the horizon. 
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For example, the initial value ip^ = ISMpL for the ^m^ip'^ model gives the background 
evolution described below and shown in Figure 7.4. As the evolution quickly reaches the 
attractor solution, the choice for ip^ is not particularly important; changing the initial 
value adds or subtracts a small number of e-foldings of evolution before the modes are 
initialised [137, 165]. 

The initial conditions are set for each k mode a few e-foldings before horizon crossing. 
This follows Salopek et al. [iOy] and is justified on the basis that the mode is sufficiently 
far inside the horizon for the Minkowski limit to be taken. This initial time, A/initl^), 
is calculated to be when 

—^ = 50. (7.18) 

(a-H)linit 

The range of e-foldings being used must include the starting point for all k modes, 
but the parameter on the right hand side, here chosen to be 50, can be changed if 
needed. We use the small wavelength solution of the first order equations described in 
Section 2.3 as the initial conditions [169], with 

^V^llinit = 1= , (7.19) 

a Jlk 



a J2k V aH 



where conformal time r] can be calculated from 



j (W/aH ~ -{aH{l - 4))"^ , (7.21) 



when Eh changes slowly. For example /cwmap is initialised about 65 e-foldings before 
the end of infiation and crosses the horizon about 5 e-foldings later. We also use these 
formulae in the calculation of the source term in Eq. (7.14) to determine the value of 
6(pi for a k mode before its numerical evolution has begun. 

We are interested in the production of second order effects by the evolution of the 
the Gaussian first order modes and we make no assumptions about the existence of 
second order perturbations before the simulation begins. Therefore, we set the initial 
condition for each second order perturbation mode to be 6ip2 = 0, 6ipl = at the time 
when the corresponding first order perturbation is initialised. One argument in favour 
of this choice of initial conditions is that far in the past the perturbations are assumed 
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to be Gaussian and therefore the second (and higher) order perturbations would be 
identically zero. 

A numerical solution for the second order perturbation equation will contain a homo- 
geneous solution and a particular solution. The homogeneous part of the solution of the 
slow roll equation, Eq. (7.13), can be calculated analytically as done in Appendix A. 3. 
On their own the initial conditions we have chosen above do not remove this homoge- 
neous solution from the result for 5ip2 in general. In order to do this, and keep only 
the particular solution to the equation, it is necessary to ensure that the homogeneous 
solution is the trivial (0, 0) solution throughout the evolution. We have not attempted 
to do this in this thesis but it is an important issue for further study in the future. 
Approaches to removing the homogeneous part of the solution include calculating a 
semi-analytic value for the second order initial conditions which equals the particular 
solution or numerically trying to select the trivial homogeneous solution by introducing 
a ramping function to the source term in Eq. (7.14). In summary, the results quoted in 
Chapter 8 for the second order perturbations include both a homogeneous and partic- 
ular solution. Extraction of either of these parts from the full result remains an issue 
for future study. 

7.3. Implementation 

The current implementation of the code is mainly in the Python'^ programming language 
(with compiled Cython components) and uses the Numerical and Scientific Python 
modules for their strong compiled array support [ ^] . The core of the model compu- 
tation is a Runge-Kutta fourth order method (see, for example, Eq. (25.5.10) in [1]). 
Following Refs. [137] and [ : ], the numerical calculation proceeds in four stages. The 
background equation (7.11), rewritten as two first order equations, is evolved from the 
specified initial state until some end time required to be after the end of the inflation- 
ary regime. The end of inflation occurs when d^a/dt^ is no longer positive and the 
parameter e'\j = Sh = —HyH becomes greater than or equal to unity (see Figure 7.4). 
This specifies a new end time for the first order run, although the simulation can run 
beyond the strict end of inflation if required. For the V{ip) = Uq + ^m^ip^ model, the 
end of inflation is set by hand to remove the need for a second, inflation-terminating 



■^Python website: http : //www . python . org 
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field. The initial conditions for the first order system are then set as above. 




Figure 7.4.: The end of inflation is determined by calculating when e\j = sh = 
—W/H = 1 (red dashed line). On the horizontal axis, A/" is the num- 
ber of e-foldings from the start of the simulation. 

The system of ordinary differential equations for the first order perturbations in 
Eq. (7.12) is integrated using a standard fourth order Runge-Kutta method. A fixed 
time step method is used in order to simplify the construction of the second order source 
term. This is also necessary since it is not known a priori which time steps would be 
required at second order if an adaptive time step system were used. The first order 
equations are separable in terms of k and so it is straightforward to run the system in 
parallel and collate the results at the end. However, as will be discussed below, the 
first order calculation is not computationally expensive in comparison with the other 
stages and only takes of the order of a few minutes for around 8000 time steps with 
AA^ = 0.01 and 1025 k modes. 

Once the first order system has been solved, the source term for the second order 
system must be calculated. As the real space equation for the source involves terms 
quadratic in the first order perturbation, it is necessary to perform a convolution in 
Fourier space, as shown in Eq. (7.3). We do not transform back into real space due 
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to the presence of both gradient operators and their inverses. Instead, the slow roll 
version of the source term integrand was used, although it is worth remarking that the 
method can also be applied to the full equation. This stage is computationally the 
most intensive, and can be run in parallel since the calculation at each time step is 
independent of the others. The nature of the convolution integral and the dependence 
of the first order perturbation on the absolute value of its arguments requires that twice 
as many k modes are calculated at first order than are desired at second order. As the 
first order calculation is computationally cheaper than the source term integration, this 
does not significantly lower the possible resolution in fc-space, which is still limited by 
the source term computation time. Once the integrand is determined, it is fed into a 
Romberg integration scheme. As for 6, which was discretised by Ng points in Eq. (7.2), 
this requires that the number of k modes is 

Nk = 2' + 1, (7.22) 

for some*^ / G Z+ . This requirement can be relaxed by opting for a less accurate and 
somewhat slower standard quadrature routine. 

The second order system is finally run with the source term and other necessary 
data being read as required from the memory or disk. The Runge-Kutta method 
calculates half time steps for each required point. For example, if y{xn) is known 
and y{xn+i) = y{xn + h) is required (for step size /i), the method will calculate the 
derivatives of y at y{xn), y{xn + h/2) and y(x„ + h). As we need to specify the source 
term at every calculated time step, the requested time step for the second order method 
must be twice that used at first order. This decreases the accuracy of the method, but 
does not require the use of splines and interpolation techniques to determine background 
and first order variables between time steps. 

The second order system is similar in run time to the first order system. However, 
the source integration is more complex and involves the integration of x Ng values 
at each time step. When = 1025 and Ng = 513, the first order evolution lasts 
around 200 seconds. The source calculation, on the other hand, takes approximately 
200 seconds for each time step. Each of the four terms A — V is approximately 16 
gigabytes in size at each time step for these values of A^^ and Ng. However, only the 
integrated result is stored for use in the second order run. This is approximately 16 



^The number of discretised k modes Nk does not need to be equal to Ng. 
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kilobytes in size for each time step. Results for each stage are stored in the open HDF5 
standard [9, 79], which can deal efficiently with large files, is very portable and allows 
for data analysis independent of the Python/Numpy programming environment. 

The full calculation contains around 8000 time steps, making the source term cal- 
culation approximately 470 hours long. Each time step is independent of the others, 
however, so parallelisation of the system is straightforward. The results in Chapter 8 
were obtained on the Virgo Cluster in the Astronomy Unit at Queen Mary, University 
of London. The code was run on ten nodes, each containing four Opteron cores with 
a clock speed of 1994Mhz. With this configuration the run time of the source term 
calculation is reduced to under twelve hours. 

7.4. Code Tests 

The numerical code has been tested in a variety of controlled circumstances in order 
to quantify the effects of different parameter options. In particular, it is important to 
establish whether the values specified for the number of discretised 6''s, Nq, the size 
of the spacing of the discretised k modes, Afc, and the range of k values significantly 
impact on the results. The sections of the code that solve ODEs are straightforward 
and follow standard algorithms. 

As mentioned above, the WMAP results [ ] use observations in the range k G 
[0.92 xlO-'^o, 3.1 X 10-58] j^^PL = [3.5xlO-^0.12]Mpc-^ We will consider three different 
k ranges both in our results and the tests of the code^: 

Ki = [1.9 X 10"^ 0.039] Mpc^i , Ak = 3.8 x lO^^Mpc^^ , 

K2 = [5.71 X 10-^ 0.12] Mpc"^ , AA; = 1.2 X lO^^Mpc"^ , 

K3 = [9.52 X 10"^ 0.39] Mpc"^ , Ak = 3.8 x lO^^Mpc"^ . (7.23) 

The first, Ki, has a very fine resolution but covers only a small portion of the WMAP 

^Thc k ranges in AfpL are: 

Ki = [0.5 X 10"^\ 1.0245 X 10"^^] A/pL , Ak = I x lO^^^MpL , 

K2 = [1.5 X 10"'^\ 3.0735 X 10-58] MpL , A/c = 3 x lO^^^MpL , 

K3 = [0.25 X 10"^", 1.02425 x 10""] A/pl , Ak = 1 x 10-^°ArpL . 
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range. The next, K2, is closest to the WMAP range and still has quite a fine resolution. 
The final range, K^, has a larger k mode step size, Ak = 1 x 10~^°MpL = 3.8 x 
10-^Mpc-\ and covers a greater range than the others. It extends to much smaller 
scales than WMAP can observe. 
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Figure 7.5.: The analytic solution of given in Eq. (7.26) for k G Ki. The value of a 
is set as 2.7 x 10^^. 



The central calculation in the code is of the convolution of the perturbations for 
the source term, Eq. (7.14). The first of the 9 dependent terms in Eq. (7.2), A, can 
be convolved analytically for certain smooth choices of 6(pi{k). Taking 6(pi{k) to be 
similar in form to the initial conditions (7.19) gives 6(pi{k) oc l/\/k with proportionality 
constant a. If denotes the following integral of the A term: 

dq'6Mq')6^,{k' - q') = 2n dq q'6^,{q')A{k\q') , (7.24) 

^min 

then substituting 6ipi{k) = a/y/k implies that 

J^(A;) = 27ra2 / dqq^ dO {k^ + q^ - 2kqcose)-^^^ sinO . (7.25) 
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(a) The relative error for different Nq, the 
number of discretised Os, keeping the other 
parameters fixed and using the range. 
The upper blue line [Ng — 129) and mid- 
dle green line {Ng = 257) have relative er- 
rors at least an order of magnitude larger 
than the lower red line [Ng — 513). 




(b) The relative error for the three different 
k ranges Ki, K2, (starting from the 
left). The parameter Afc is set equal to 1 x 
10-6iMpL,3 X 10-'^U/pL,l X 10-60 Af, 
respectively. 
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Figure 7.6.: Comparison of relative errors in for different Ng and k ranges. 



This lias the analytic solution 



'iSk 



3F 



log 



+ Vk n 

+ log 



TT 



H arctan . 

2 \vk — kit 



Vk 



Vkmin + k + a/A; 

mir 



(Sk -\- 8k^^^ (^^k -\- A^max \^~knmx k 



Ukk^ 



k ~\~ A^niax ~l~ \/^niax A? 



+ \/kmm (S/c + S/Cj^jjj) (^\^k -\- fcmin + ^/k ~kn 
+ Mfc/Cuiin (^^/k^^\-~k^^ — ^/k — /Cmin^ J' • 



(7.26) 



The k dependence of can be seen in Figure 7.5. We have tested our code against 
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this analytic solution for various combinations of k ranges and Nq. The relative error 

[analytic — calculated] , , 

I analytic I 

is small for all the tested cases, but certain combinations of parameters turn out to 
be more accurate than others. The relative errors of all the following results are not 
affected by the choice of a so we will keep its numerical value fixed throughout as 
2.7 X 10^^. 
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Figure 7.7.: The relative error in the integral for different values of Ak. The other 
parameters are fixed: fcmin = 1 x 10~^°Mpl, A^'^ = 1025 and No = 513. The 
value of Ak is less than fcmin for the upper blue line {Ak = 1 x 10"^^Mpl) 
and the middle green line [Ak = 3 x 10"^^Mpl). These have relative errors 
at least an order of magnitude larger than the lower red line for which 
Ak = k^i^ = 1 X 10-6°MpL. 

We first tested the effect of changing Ng, the number of samples of the 9 range [0, vr]. 
Figure 7.6a plots these results for the k range K3 with Ak = 1 x lO^^^AfpL. Only 
three values of Ng are shown for clarity. It can be seen that increasing Ng decreases 
the relative error when the other parameters are kept constant, as one might expect. 
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As mentioned above the choice of k range is especially important as the convolution 
of the terms depends strongly on the minimum and maximum values of this range. 
Indeed, this is clear from the analytic solution in Eq. (7.26). Figure 7.6b shows the 
difference in relative error for the three different k ranges described above with Ak = 
3.8 X 10-^ 1.2 X 10"^ and 3.8 x IQ-^Mpc"^ {Ak = 1 x 10-^\ 3 x 10-^\ 1 x IO-^^Mpl), 
respectively. The accuracy is similar in all three cases. 

Another important check is whether the resolution of the k range is fine enough. 
Varying Ak can not be done in isolation if the constraint (7.22) for A^^ is to be satisfied. 
For this test the end of the k range was changed with Ak but the other parameters 
were kept fixed at k^i^ = 1 x IQ-^^Mpl = 3.8 x 10-^Mpc~\ Nk = 1025 and Ng = 513. 
Figure 7.7 plots these results again for three indicative values. For Ak < A;min, there 
is a marked degradation in the accuracy of the method for the upper two lines. This 
is understandable as many interpolations of multiples of Ak below /Cmin will be set to 
zero. Once Ak is greater than fcmm, the relative error is insensitive to further increases 
in the value of Ak. (This is not shown in the figure.) 
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Figure 7.8.: The relative error in the integral for three different values of Ak. In 
contrast to Figure 7.7, kmin = 1 x IO'^^Mpl = 3.8 x lO^^Mpc"^ < AA; for 
each case. 
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The analytic solutions for the B, C and V terms are given in Appendix A. 4. The 
relative errors between the analytic and calculated values for Jg, and are shown 
in Figure 7.9 for the three final k ranges, with /3 = 10~^^. The errors for the Jg- and 
Ijy terms are very small, being of the order of 10~® and 10^^, respectively. The relative 
error for the term is larger, especially for small k values. However, the error is still 
below 0.08% for each of the Ki, K2 and ranges. 

It should be noted that these tests only show the relative errors in the computation 
of integrals of the four terms in Eq. (7.2). They do not represent errors for the full 
calculation. However, they do show that the accuracy is good compared with the 
analytic result. 

7.5. Discussion 

This chapter has described the implementation of the numerical calculation of second 
order cosmological perturbations. The Klein-Gordon equations in Chapter 6 are the 
central focus of this computation. In Section 7.2 these equations were rewritten using 
M as the time variable, a choice more suitable for numerical work. The convolution 
integrals in Eq. (6.25) can be expressed in spherical polar coordinates and split into four 
sub-terms A-V in Eqs. (7.2) and (7.15). Computing the source term Eq. (7.14), which 
is written using A-T), is the most complex and time consuming part of the calculation. 

To demonstrate the numerical code, four different potentials have been chosen and 
these were described in Section 7.2.1. One parameter for each potential was determined 
by comparing the calculated with the WMAP5 value. 

In Section 7.2.2 the initial conditions for the computed quantities were explained. 
Each k mode is initialised well inside the horizon using the Bunch-Davies vacuum 
conditions from Section 2.3. The second order perturbations are initially set to zero. 
This choice concentrates focus on the generation of second order effects during the 
observable period of inflation. 

The technical implementation of the code was discussed in Section 7.3. There are 
four stages in the procedure. First, the background fields are evolved over a specified 
time period. The end time of inflation is then determined by the condition Sh = ^ 
and the scale factor calculated for this time. With this information the initialisation 
of the first order modes can be performed. In the second stage, the first order pertur- 
bation equations are solved. These results are used in the third stage to calculate the 
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(c) The relative error in for each k range. 
Figure 7.9.: The relative errors in /g, and Ijj for each of the k ranges with /3 = 10 
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source term (7.14) at each time step and for each k value. Finally, the second order 
perturbation equations are solved using the source term results. 

To test the source term calculation the numerical results were compared to ana- 
lytic solutions in Section 7.4. Numerical parameters such as Ne and Nk were set by 
minimising the relative error between the two approaches. 

In Chapter 6 the evolution equations for second order perturbations were introduced. 
In this chapter the practical implementation of a numerical calculation of these pertur- 
bations was discussed. In Chapter 8 the results of this numerical calculation will be 
examined and the next steps towards an improved procedure will be described. 



8. Results and Future Work 



8.1. Introduction 

The main result of Part II of this thesis is the numerical integration of the Klein-Gordon 
equation of motion for second order scalar field perturbations, Eq. (6.25). The slow roll 
approximation of the source term for second order perturbations was employed, but 
the complete versions of the evolution equations were used for the background and 
first order perturbations. In this chapter the results of the numerical calculation will 
be presented. This represents the first step towards a full calculation of the Klein- 
Gordon equation at second order. In addition to the new results obtained, plans will 
be described for future work aimed at improving the numerical system and increasing 
its range of apphcabihty. 

As a proof of concept, the numerical system was tested with four different potentials, 
V{(p) = \m?(^'^, (^^^ and Uo + ^mlcp'^, and results computed across three different 

k ranges. As expected, the second order perturbation for a single, slowly rolling infiaton 
field that we have calculated is extremely small in comparison with the first order term. 
However, there are differences apparent between the potentials, which will be outlined 
in Section 8.2.3. 

We have listed the potential parameters m. A, a and mo in Table 7.1. These were 
found using the WMAP5 normahsation at ^wmap = 0.002Mpc-^ = 5.25 x IQ-^^MpL 
[104]. We have also outlined in Eq. (7.23) the three k ranges that have been used: 

Ki = [1.9 X 10"^ 0.039] Mpc"^ , Ak = 3.8 x lO^^Mpc^^ , 
K2 = [5.71 X 10"^ 0.12] Mpc"^ , Ak = 1.2 x lO^^Mpc^^ , 
/sTs = [9.52 X 10^^ 0.39] Mpc-\ AA; = 3.8 x lO^^Mpc"^ (8.1) 

Many of the results will be quoted for fcwMAP which lies in all three of these ranges. 
Given that the first order perturbations for the chosen potentials produce an almost 
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scale invariant power spectrum with no running, it is no surprise that the results from 
the three different k ranges are very similar. The second order source term is somewhat 
dependent on the lower bound of k (upper bound on size). This is also to be expected 
and in the scale invariant case a logarithmic divergence can be shown to exist [124]. 
We have implemented an arbitrary sharp cutoff at /cmm; below which 6(fi is taken to 
be zero. As mentioned in Chapter 7, there is some evidence to suggest that a similar 
cutoflF might be supported by the WMAP data [05, 1S7]. 

In Section 8.2, the numerical results for the computation described in Chapter 7 are 
presented. Comparisons of the results from the four different test potentials will be 
made in Section 8.2.3. Since this represents the first stage towards a full calculation 
of the source term, the next steps that will be required are outlined in Section 8.3. 
Finally, in Section 8.4 we discuss some of the consequences of our results. 

8.2. Results 

8.2.1. Results for V{ip) = \m^^'^ 

At first order the solutions obtained for the quadratic potential agree with previous 
work in Refs. [137, 165, 169]. Oscillations are damped until horizon crossing (when k = 
aH) after which the curvature perturbation becomes conserved. Figure 8.1 shows the 
evolution of the real and imaginary parts of the first order perturbation from when the 
initial conditions are set, at k/aH = 50, to just after horizon crossing. The horizontal 
axis for most of the following figures parametrises the number of e-foldings remaining 
until the end of inflation (ATend — A/"), instead of the time variable M employed in the 
calculations. 

Figure 8.2 shows the evolution of the second order perturbations for the scale /cwmap- 
As mentioned above, the overall amplitude of the second order perturbation is many 
orders of magnitude smaller than the corresponding first order one. The results given 
for 6(p2 in this chapter are for the full solution which includes a homogeneous and 
particular solution as described in Section 7.2.2. In Figures 8.1 and 8.2 the field values 
have been rescaled by A;^/^/(-\/27r) to allow for a better appreciation of the magnitude 
of the resulting power spectra. 

The source term S{k^) is calculated using Eq. (7.4) at each time step using the results 
of the first order and background simulations. This term drives the production of second 
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Figure 8.1.: The first order perturbation S^pi rescaled by k"^/"^ / {\/2'k) from tlie beginning 
of tlie simulation until around horizon crossing (red dot-dashed line). The 
real (blue) and imaginary (green dashed) parts of the perturbation are 
shown for the scale /cwmap- 
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Figure 8.2.: The real (blue line) and imaginary (green dashed) components of the second 
order perturbation 5v^2(^wmap) from the beginning of the simulation until 
around the time of horizon exit (red dot-dashed line). 



8.2: Results 



134 



order perturbations as shown in Eqs. (6.25) and (7.13). Figure 8.3a shows the absolute 
magnitude of the source term for a single k mode, /cwmap, for all time steps calculated. 
The source term is large at early times, and closely follows the form of the spectrum of 
the first order perturbations, as can be seen from Figure 8.3b. Figure 8.4 shows how 
the source term depends on the choice of k range. After horizon crossing, the source 
term is independent of the specific choice of Ki {i = 1,2,3). Before horizon crossing, 
however, there is a strict hierarchy with the smaller k ranges, Ki and K2, leading to 
smaller source contributions. As discussed in Section 7.4, Ak should be at least as 
large as k^m in order for the error to be reduced to a minimum. In Figure 8.5 the 
source term is plotted at three different values of k for the range Ki. As k increases, or 
equivalently the length scale decreases, the magnitude of the source term after horizon 
crossing decreases. 

It is informative to compare the magnitude of the source term with the other terms 
in the second order evolution equation (7.13). We denote these other terms by T: 

T(F) = (3 + Sdik') + (^) ' SMk') + - 247rG(v.S)^) ^Mk') ■ (8.2) 

Figure 8.6 then shows the absolute magnitude of both S and T. It is clear that for the 
scale /cwMAP the contributions to the source term are only of comparable magnitude 
during the early stages of the simulation. Figure 8.7 shows a comparison of |5'|/|T| 
for three different k values. The magnitude of S is closer to that of the rest of the 
terms for the larger k mode. A priori, the range of k modes where S will be large for a 
particular chosen potential is not known. However, once the relevant values of k have 
been determined, it may be possible to significantly reduce the time required for the 
simulation by restricting the calculation of S to those regions where it is important. 
Figure 8.7 shows that it is not possible to arbitrarily ignore the contribution of 5* either 
inside or outside the horizon. The full calculation on sub- and super-horizon scales is 
important for the evolution of the second order perturbation at different scales. 

In Figure 8.8 the value of 15*1 at the initialisation time for each k mode is shown. 
The initial magnitude of the source term is much smaller for larger values of k (smaller 
scales). Because the smaller k modes begin their evolution earlier, the relative difference 
in 15*1 is not as pronounced when measured at a single time step (see for example 
Figure 8.5). It should also be remembered that the magnitude of the other terms in 
the second order ODE is small for larger k modes as shown by the ratio I^I/ITI in 
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(a) Absolute magnitude of the source term. 
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(b) Power spectrum of first order scalar perturbations Vg^^ — ^^|(5(^ip. 

Figure 8.3.: Source term and first order power spectrum for tlie WMAP pivot scale 
Wmap- 
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Figure 8.4.: A comparison of the source term (7.14), for the scale ^wmap = 5.25 x 
10~^^Mpi, over the three different ranges Ki, K2 and K^, which were 
specified in Eq. (8.1). Before horizon crossing there is a significant difference 
in the amplitude of the source term for /cwmap- After horizon crossing, 
however, the magnitude of S is independent of the choice of K^. 
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Figure 8.5.: The source term (7.14) for three different k values in the Ki range, including 
the WMAP pivot scale, ^wmap = 5.25 x IQ-^^Mpl (middle green line). As 
the value of k increases or equivalently the scale decreases, the magnitude 
of the source term decreases. The calculation of the source term for each 
k value starts from the time step at which the corresponding first order 
perturbation is initialised, i.e., when k/aH = 50. 
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Figure 8.6.: The source term (lower blue line), as defined in Eq. (7.14), is compared 
with the T term (upper green line), as defined in Eq. (8.2), for /cwmap- The 
source term is of comparable magnitude at the beginning of the simulation. 
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Figure 8.7.: The quotient of the 5* term, Eq. (7.14), and the T term, Eq. (8.2), for 
three different k values in the range Ki, including the WMAP pivot scale, 
^WMAP = 5.25 X 10~^''MpL. For small values of k the source term is not 
comparable to the magnitude of T after horizon crossing. However, for 
larger k values (smaller scales) the two terms have comparable magnitude. 
It is, therefore, important to calculate the source term over the full range 
of e-foldings. 
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Figure 8.7. 
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Figure 8.8.: The absolute magnitude of the source term at the initial start time for each 
k value when k/aH = 50 deep inside the horizon. The results are for the 
range Ki. 

The source term can also be compared at different time steps over all the k values. 
In Figure 8.9 the green line shows l^l when all 6(fi modes have started to evolve. The 
lower red line illustrates 15*1 after all modes have exited the horizon, around 52 e-foldings 
before the end of inflation. 

8.2.2. Comparison of V{(p) = ^m^cp'^ Results with Analytic 
Solution 

In this section results for the quadratic model will be compared with an analytic solution 
for this model. However, an analytical result is difficult to obtain for the case of the full 
first order solution in terms of Hankel functions with the phase information included. 
The analytical solution we will use, therefore, is the non-interacting de Sitter space 
solution with the phase information ignored. The first order perturbations are then 
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Figure 8.9.: The absolute magnitude of the source term for all k values in the range 
Ki at two different time steps. The green line shows IS"! when all modes 
have been initialised. The lower red dashed line shows IS"! approximately 
52 e-foldings before the end of inflation, when all modes have exited the 
horizon. 
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given by 

and the derivative in terms of is 

6^\{v,k') = (l + ^)-^^^/k. (8.4) 

aV2k\ kr,J\ aHr, J amV2 

The source term is found using Eq. (7.14) and the values of the background quantities. 
The analytical solution of Eq. (7.14) for this choice of first order solution is given in 
Appendix A. 5 in Eqs. (A.35-A.42). 

In Figure 8.10 the analytical and calculated solutions are plotted for one timestep 
about 60 e-foldings before the end of inflation. At a single time step the correlation 
between the two solutions for S is very good. There is a deviation at small values of k 
due to the analytical solution getting rapidly smaller as k approaches zero. This is not 
replicated in the calculated version. However, this only strongly affects the result for 
the smallest values of k and for /cwmap for example the relative error of the calculated 
solution compared to the analytical solution is about 10~^. The relative error of the 
calculated solution is shown in Figure 8.11 for a single time step. 

The analytical and calculated values for the source term can also be compared for a 
single k value across a range of time steps. In Figure 8.12 the absolute magnitude of the 
source term for /cwmap is plotted for a range of a few e-foldings before horizon crossing. 
The analytical and calculated results are extremely similar and not distinguishable in 
the plot. The relative error of the calculated solution is around 10~^. Because the first 
order perturbations do not include any phase information the result is much smoother 
than the result generated using the full first order phase information. However, as the 
phase angle is a function of — g*| it cannot be trivially ignored in the computation 
of the full convolution integral. 

8.2.3. Comparison of Models 

All the results quoted so far have been for the quadratic potential. In this section the 
results for all four potentials will be compared using the K2 range. Figure 8.13 shows 
the power spectrum of first order curvature perturbations, V^_^, for each potential. The 
|m^(y9^, \X(p^ and acp^ models all clearly have a red spectrum with < 1. On the 
other hand, the Uq + \m'^ip'^ model has a blue spectrum {ris > 1) when Uq to chosen to 
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Figure 8.10.: A comparison of the analytical and calculated solutions for the source term 
at one time step, approximately 60 e-foldings before the end of inflation. 
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Figure 8.11.: The relative error of the calculated solution compared to the analytical 
solution for the source term. The error is shown for all k values at one 
time step, approximately 60 e-foldings before the end of inflation. 
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Figure 8.12.: A comparison of the analytical and calculated results for the source term 
of the fcwMAP mode before horizon crossing. 
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Figure 8.13.: Comparison of the power spectrum V-^^ for the four different models. 

The three models with potentials |m^y9^, |A<y9^ and a^p^ have red spectra 
{ug < 1) while the Uq + \m1ip'^ model has a blue spectrum (n^ > 1). 
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Potential Ug Us — I 

|mV^ 0.965 -0.035 

|A(/7^ 0.949 -0.051 

aipl 0.977 -0.023 

Uo + \ml^^ 1.002 0.002 

Table 8.1.: The spectral index for scalar perturbations for each of the four potentials 
used. These values are calculated for the /cwmap scale, five e-foldings after 
it crosses the horizon. The potential parameters are listed in Table 7.1. The 
value t/o = 5 X 10~^°MpL was chosen to ensure a blue spectrum {rig > 1). 



be 5 X 10~^''MpL, as specified in Section 7.2.1. The values of rig obtained for the four 
potentials are given in Table 8.1. 

The source term for each model is shown separately in Figure 8.14 for ^wmap using 
the K2 range^. Although these terms are qualitatively similar, differences between 
them are apparent. Figure 8.15 brings together the source terms at ^wmap to enable 
a direct comparison to be made. The fcwMAP mode begins at different times for the 
different models. Each result is therefore plotted in terms of the initialisation time for 
that mode. This change in duration is a consequence of allowing H to evolve during 
the calculation. 

The source term results for the quadratic and quartic potentials are very similar. 
Indeed, from horizon crossing to near the end of inflation the results appear to coincide. 
The |A(/9^ mode has a slightly longer duration and at late times is reduced in comparison 
with the |m^v9^ one. Figure 8.16 shows that at early times the relationship is more 
complicated with the |A<y9^ mode being larger for a significant period. 

2 

In the early stages the amplitude of the V{ip) = crips model is very similar to the other 
two results described above. After horizon crossing, however, there is a significant drop 
in the amplitude of S in comparison with the ^^m^ip^ and \\ip'^ models. This continues 
until late in the evolution when 15*1 increases swiftly to reach levels above the others. 
The duration of the mode in this model is shorter than in the other two models described 
so far. 

The fourth model, with potential V{({)) = Uq + \m1({P', is a contrived toy model. As 
described in Section 7.2.1, in order to perform the single field calculation, the end time 



"'^ These plots use a different k range to tlie ones comparing V{ip) — ^w?ip^ and |;A(p"' in Ref. [^4]. 
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(c) V{v) = (d) V{^) ^Uo + \ml^^ 



Figure 8.14.: Plots of the source term for the four different potentials studied. 
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Figure 8.15.: Comparison of the source term evolution for the four different models. 

After horizon crossing the magnitude of the source term is larger for the 
quadratic and quartic models than for the other two. Towards the end 
of the numerical calculation there is a marked increase in 15*1 for three of 
the models as Eh increases towards unity. The end time of inflation is 
specified by hand for the contrived toy model, so this effect is not seen. 
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Figure 8.16.: Comparison of the source term evolution for the four different models at 
early times. This figure highlights the early evolution of the four models 
shown in Figure 8.15. Before horizon crossing the magnitude of the source 
term is comparable for each model. After horizon crossing differences 
between the models become apparent. 
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of inflation must be specified by hand. In this simulation ~ 8 is taken as the end 
time. The potential is extremely flat in this region and the effect of this can be seen 
in the source term of the model. Before horizon crossing it is of comparable magnitude 
to the other terms. However, a steep decrease in IS"! ensures that it is a few orders 
of magnitude smaller than the other terms after horizon crossing. In contrast to the 
behaviour of the other models, the source term does not increase close to the end of 
inflation. This is due to the enforced end time cut-off which means that Eh does not 
become large. 

In this section we have described the results for four different single field potentials. 
As expected for single field slow roll models they exhibit similar properties. In the next 
section plans to extend the calculation to deal with more interesting models will be 
outlined. 

8.3. Future Directions 

There are many possible ways to improve the program outlined in Chapter 7. Chief 
amongst these is the implementation of the full second order source term given in 
Eqs. (6.20) and (6.21). As we have seen the slow roll approximation is very helpful 
in reducing the equations of motion to a manageable size. However, many interesting 
models break the assumptions of a slowly rolling field and to investigate these models 
it is necessary to use the full field equations. 

Models in which the field potential is not smooth due to the presence of a feature are 
particularly interesting examples of single field inflation for which slow roll is broken. 
As the derivatives of the potential can be large around the feature, these models must 
necessarily be handled without assumptions about the size of the slow roll parameters. 
In Ref. [3] a model with a step potential was proposed which takes the form 



where ips, c and d parametrise the location, height and width of the step feature. A 
bump model has also been proposed [12], with potential 




(8.5) 




(8.6) 
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where again iph, c and d parametrise the feature. At first order these models introduce 
noticeable differences in the scalar power spectrum. They are also known to be able 
to produce significant amounts of non-Gaussianity in shapes which are not similar to 
either the local or equilateral types described in Section 2.6 [11, 42]. 

It will also be important to go beyond slow roll in the multiple field case. To obtain 
analytic results, the study of multi-field models has often been restricted to those with 
either sum or product separable potentials. Even very simple models with two fields 
such as the double inflation model with the potential given by [184, 197] 



can violate slow roll when the fields (p and x ^^e close to equality. To go some way 
towards considering the full range of possible multi-field models with arbitrary infla- 
tionary potentials then requires that the full non-slow roll evolution equations are used. 

As far as the implementation of the code is concerned, the extension to the non-slow 
roll single field case is the next step. Although clearly more complicated than the slow 
roll case of Eq. (7.13), only three more 9 dependent terms need to be added to the A-V 
terms listed in Eq. (7.2). The four potentials considered above all result in slow roll 
inflation. Therefore, it is not expected that using the full source equation will result in 
an appreciably different outcome in these models until near the end of the inflationary 
phase. Once the field has stopped rolling slowly, new observable features are expected 
to arise, as is indeed the case at first order. 

Eqs. (6.20) and (6.21) must be written in terms of A/", with the 6 dependent terms 
grouped together, in order to set up the numerical system completely at second order. 
The main equation becomes 



^ 2 2,^ 22 

2 2 ^ 



(8.7) 
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V^^ + 9.TiG[2iplV^ + MG{ipl) VA 5<^2(A;0 + 5fuii(A:0=0, (8.8) 
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where the full source equation is given by 
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The Ffuu term in Eq. (8.9) requires the use of three further 9 integrals in addition to 
those presented in Eq. (7.2). These take the form 



S{k\q')= I cos\e) sm{e)6^i{k' - q')de , 
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(8.10) 



It is worth noting that the term sin {e)/\k'-q'\^ tends to zero in the limit where k = q 
and ^ — )■ 0. The Ff^u term can now be written in terms oi T and Q: 



full 



MG 1 
(27r)2 {aHY 



dq q^ 



+ (87rG)Q^S(^ + |^) Vi(g') 



+ 



-2Q(87rG)(^S)2|5^I(gO 



B{k\q' 



8.3: Future Directions 



154 



+ 



2 + (87rG)(^Sf 



4 



1 



+ 



2aH 



1 



) 



) 



+ 2Q%,{q^) + (2- - I) ^t(aiJ) Vi(gO 



q \ q 




(8.11) 



Eqs. (8.9) and (8.11) are clearly more complicated than the slow roll versions used in 
Chapter 7. The numerical complexity is not significantly greater, however, once the 
three terms in Eq. (8.10) have been calculated. The running time of the full calculation 
will clearly be a significant constraint. 

With this in mind, the performance of the numerical simulation could be improved 
by analysing the most time consuming processes and investigating what optimisations 
might be implemented. The current, perhaps inelegant, procedure will allow any perfor- 
mance improvements to be benchmarked for accuracy as well as for speed. As discussed 
above, A^^ was set to 1025 for the test runs. This provides good coverage of the WMAP 
k range, but it is not clear whether it sufficiently approximates the integral to infinity 
for the source term. Logistical factors, including the running time and memory usage of 
the code, restrict the choice of A^^. By optimising the routines for reduced memory and 
increased speed it is hoped that the range of scales can be extended and the resolution 
enhanced. 

Beyond these considerations, the next significant step is to implement a multi-field 
version of the system. This would allow the investigation of models that inherently 
produce large second order perturbations. In Ref. [133] the second order Klein-Gordon 
equation for multiple fields was presented and upgrading the simulation to use these 
equations should be a straight-forward (if lengthy) process. Extending the current data- 
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structures and routines to a fixed number of extra fieids wiii increase ttie numerical 
complexity and the run-time of the code. 

For example, let us suppose that the second order perturbations of two scalar fields, 
(p and X, are to be calculated. Let V denote the potential and Vo its background value. 
As the coding environment we have used can easily handle arrays of variables, it is 
useful to write the equations in vector form. The following definitions will be used: 



/ Sipi 
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(8.12) 
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In conformal time the Priedmann equation becomes 



Q(<^o')^<^o' + aVo^ , (8.14) 



where denotes the transpose of </5. The background vector equation of motion is 
given by 

V?o" + aT^V'o' + a Vi = , (8.15) 

where is the zero vector. The first order vector equation takes the form 



(8.16) 



where 1 is the identity matrix. 

We will outline the second order vector equation using the slow roll approximation. 
In the multi-field case there are many more slow roll parameters than in the single field 
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scenario. Extending the definition of Eh in Eq. (6.24) to two fields gives 



(8.17) 
(8.18) 



There are now four ?7//-type parameters corresponding to the different combinations of 
second derivatives of V . These can be written together in matrix form as 



(8.19) 



where I,J = (p, x- The magnitude of rju is only small in the adiabatic direction, so 
terms including rjjj are included when making the slow roll approximation [ ]. 
The second order, slow roll, vector equation for the perturbations is given by 

6cp2'\k') + 2H8ifi2'{k') + {ei + a'V2 - 2AnGipo' {^o'f) Scp2{k') 

+ S{k')=0, (8.20) 

where the slow roll source term equation is 
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Following the method of Section 7.2, the d^p integral is evaluated and the d^q integral 
is written in spherical polar coordinates. The 6 dependent terms, which are equivalent 
to Eq. (7.2) in the single field case, are given by 
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C{k\q') 
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(8.23) 



The first two equations are not vector equations but are needed for the explicit matrix 
term in Eq. (8.21). We rewrite that equation with these definitions to obtain: 
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This expression for S reduces to Eq. (7.4) when only one field is considered. At least in 
the slow roll case, the multi-field source term equation is not considerably more complex 
than the single field one. The extra numerical complexity arises from the calculation 
of the new 6 dependent terms in Eq. (8.23). 

8.4. Discussion 

Part II of this thesis has described the numerical solution of the evolution equations 
for second order scalar perturbations. The closed form of the Klein-Gordon equation 
(6.25) has been employed for the first time. We have shown that direct calculation of the 
field perturbations beyond first order using perturbation theory is readily achievable, 
although it is non-trivial. 

This first demonstration has been limited to considering the slow roll approximation 
of the source term in Eq. (6.25) which is quadratic in first order perturbations. Slow 
roll has not been imposed on the background or first order equations. Four different 
potentials were used to demonstrate the capabilities of the system. The singularity at 
k = 0, which arises as larger and larger scales are considered, is avoided by implementing 
a cutoff at small wavenumbers below fcmm- This is a pragmatic choice necessary for the 
calculation. It is also necessary to specify a maximum value of k. This choice is dictated 
by computational resources and with reference to observationally relevant scales. In this 
demonstration, k ranges have been used which are comparable with the scales observed 
by the WMAP satellite. By comparing the analytical results of the convolution integral 
with the numerical calculation, values of the parameters Ng, Nk and Ak were chosen in 
such a way that the numerical error was minimised. The convolution scheme that has 
been implemented works best when Ak > fcmin- 

We have seen explicitly that the second order calculations for the chosen potentials 
can be completed once the cut-off for k^i^i is imposed. As expected for these potentials 
the magnitude of second order perturbations is extremely suppressed in the slowly 
rolling regime, in comparison with the first order amplitude. We have also shown that 
the evolution of the source term during the infiationary regime can be readily calculated. 

By computing the perturbations to second order, we have direct access to the non- 
Gaussianity of 6ip. When used to investigate models that predict a large non-linearity 
parameter, /nl, this technique could yield greater insight into the formation and devel- 
opment of the non-Gaussian contributions by studying the effects of the different terms 
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in the source equation (7.14). It was shown recently that /nl can be calculated directly 
from the field equations [149, 175], instead of using the standard method based on a 
Lagrangian formalism [130]. The method presented here will therefore eventually allow 
a full numerical calculation of /nl to be made. 

Our numerical code evolves the second order perturbation itself and gives an insight 
into how this field behaves through the full course of the infiationary era. This is in 
contrast to other approaches which only consider the result for the three point function 
of the field, or alternatively of the curvature perturbation. The computational system 
handles perturbations with scales both inside and outside the horizon. Any effects of 
horizon crossing are visible and no assumptions need to be made about the form of the 
solution inside the horizon. 

The numerical code, when developed with the full equation, will not require any 
simplifying assumptions about the form of the potential used. This allows models 
which are not amenable to analytic analysis to be examined. Examples of models 
which require consideration beyond the slow roll approximation include single field 
models with a step or other feature in the potential, and multi-field double infiation 
models where the field values are roughly equal. 

The code we have developed is also applicable in other physical circumstances. Be- 
yond scalar perturbations the form of the source term is similar in other interesting 
cosmological physics. The generation and evolution of non-Gaussian curvature per- 
turbations is, of course, directly related to the behaviour of the second order scalars 
as has been described in Section 2.6 and Section 6.3. Investigating and classifying 
non-Gaussian signatures for infiationary models is the main goal of our future work. 

The generation of vorticity in a cosmological setting has physical parallels with the 
equations we have studied. This second order effect arises through the vector perturba- 
tions which we have not considered in this thesis. Vorticity in the early universe could 
also lead to the generation of primordial magnetic fields, an area which is of increasing 
interest [29, 18]. The wave equations for tensor mode perturbations also exhibit the 
same form as the scalar equations with a source term at second order. The code we 
have developed could be modified to examine the behaviour of gravitational waves in 
the early universe at second order. 

In summary, we have demonstrated that numerically solving the closed Klein-Gordon 
equation for second order perturbations is possible. The slow roll version of the source 
term was used in the calculation, but as described in Section 8.3, the extension of 
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the system to include the full source term is achievable. The analytic and numerical 
solutions for the convolution terms were compared directly and found to be in good 
agreement. The models used have been shown to have negligible second order pertur- 
bations in line with known analytic results. 
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9. Conclusion and Discussion 



As the number of viable cosmological models increases, the need to constrain them 
becomes more important. At the same time, the quantity and quality of observational 
data continue to improve. There now exists the opportunity to go beyond linear sta- 
tistical analyses and confront the predictions of models with observational data from 
the non-linear regime. In this thesis both analytic and numerical methods have been 
developed to constrain inflationary models. 

The framework used in this thesis is the Friedmann-Robertson- Walker universe, re- 
viewed in Chapter 2. During the accelerated expansion of the inflationary period, 
quantum fluctuations seeded energy density variations, which in turn gave rise to the 
diverse structure of the present universe. First order cosmological perturbation theory 
is necessary to describe the evolution of these fluctuations. The main observable quan- 
tities can be calculated at horizon crossing by using the Bunch-Davies vacuum initial 
conditions. The departure of the perturbations from a purely Gaussian random field is 
parametrised by /nl, which is described in two limits, local and equilateral. In addition 
to canonical actions, we also introduced non-canonical models, for which the speed of 
sound of the perturbations plays a crucial role. When the sound speed is small, the 
amplitude of /^l for these models is large. 

In Part I, analytic methods were developed to constrain string theory inspired non- 
canonical inflationary models. The Dirac-Born-Infeld scenario was outlined in Chap- 
ter 3. In this model, a D3-brane propagates in a six-dimensional warped throat. The 
radial position of the brane from the tip of the throat assumes the role of the inflaton 
field. The non-canonical nature of the DEI action restricts the kinetic energy of this 
field no matter how steep the potential. This allows an inflationary period of sufficient 
duration to occur. This model has been widely regarded as a very promising realisation 
of an inflationary model in a string-theory context. 

In Ref. ['0], Baumann & McAllister used the Lyth bound [123] to limit the tensor- 
scalar ratio. Their analysis was based on the conservative assumption that the brane 
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could not propagate further than the full length of the throat. In Chapter 4, we showed 
that this bound can be tightened by applying it over the portion of the throat through 
which the brane passes during the directly observable stage of inflation. Restricting 
the field variation to be over these approximately four e-foldings constrains r to be less 
than 10""^ for standard parameter values. 

The most optimistic estimates of advances in experimental techniques and data anal- 
ysis, including foreground reduction, indicate that observations of r > 10~^ might be 
achievable in the future [21, 190]. Therefore, the new bound in Eq. (4.12) immediately 
rules out the observation of a tensor mode signal from this model. 

In addition to this, we also derived a lower bound on r in Eq. (4.20). This depends 
on observable quantities, namely the scalar spectral index and the equilateral non- 
Gaussianity. Saturating the WMAP5 observational limit on and taking the best 
fit value for n^, we found that the most conservative lower limit is r > 0.005. This is 
clearly incompatible with the previously derived upper bound. Therefore, for standard 
parameter values, the D3-brane DBI scenario is not viable. Numerical simulations by 
Peiris et al. in Ref. [15 1] have demonstrated the lower bound, in the relativistic limit, 
using the Hamiltonian flow approach. 

In Section 4.4, a phenomenological approach was taken to easing the upper bound 
on the tensor-scalar ratio. By considering a DBI-type action with unspecified field 
functions, /j, we showed that the generalised lower and upper bounds can be consistent 
if the product of /a and fs is sufficiently large on observable scales. This provides 
a guide to the types of models which could evade the inconsistency of the bounds on 
r. For more general models with a non-canonical action, a bound on r which relates 
the geometry of the throat, the number of e-foldings of observable inflation, and the 
derivatives of the action has been derived in Eq. (4.39). This bound, although it does 
not in general relate to observational quantities can be used when the details of a 
particular physical model are known. 

The discovery of the incompatible bounds on r for DBI inflation has had a noticeable 
impact on the research community, spurring interest in finding models which evade 
these bounds. Many such models have been proposed with varying degrees of success. In 
Section 4.5 these were categorised according to whether they featured single or multiple 
fields, and single or multiple branes. Some of these models are still constrained by the 
bounds on r but not to the same extent as the standard DBI scenario. For example 
the parameter space of the models with wrapped brane configurations is still extremely 
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limited by the observational values from WMAP5 [ ']. For other models an analysis in 
terms of the bounds derived in this thesis has yet to be undertaken. As the observational 
limits on and r continue to improve, an important step in ensuring the validity of 
DBI based models is to check whether equivalent bounds to those derived here exist, 
and whether they can be met for any significant proportion of the parameter space. 

General actions which relax the bounds on r, were derived in Chapter 5. We found 
a class of actions similar in form to that of the DBI model. However, instead of a 
square-root in the kinetic term, the index of the main term of these actions depends 
on the constant of proportionality between A and the sound speed of inflaton fluctu- 
ations. The upper bound on r can be derived for these general actions and when the 
index is below the critical value of 1/2 (corresponding to the standard DBI scenario), 
this bound is significantly relaxed. When new models are proposed in the future, our 
phenomenological derivation of this family of actions will allow those models for which 
the bound on r is less stringent to be easily identified. 

One such model is the single field, multi-coincident brane scenario of Thomas & Ward 
[194]. When n D3-branes propagate in a throat, the non-Abelian interactions between 
the branes result in major departures from the single brane case. This is in contrast to 
the non-interacting branes model, in which the total action is simply the sum of copies 
of the single brane action. 

In the limit of a large number of branes being coincident, the effective action is similar 
to n times the single brane model, with the addition of a fuzzy potential term. Indeed 
this model is of the type considered phenomeno logically in Section 4.4 and the bounds 
derived in that section can be applied. These bounds on r can be somewhat relaxed 
when this potential is large, but the model is still strongly constrained by current 
observations. For an AdSr, throat, standard parameter choices limit the number of 
allowed branes to be less than 150, at which point the assumption of arbitrarily large 
n becomes questionable. 

More promising is the finite n limit of the coincident brane model. The non-Abelian 
nature of the interactions leads, in this case, to a recursive relation for the n-brane 
action in terms of the n = 2 one. In Section 5.5, we showed that the action for finite 
n is one of the class of bound-evading actions described above. This identification is 
possible because the last term of the recursive sum dominates in the relativistic limit. 
This approximation is valid at least when n < 10 and the backreaction of the multiple 
branes is kept well under control for this range of n. 
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Although the bounds on r are eased for this model, we showed that observations 
strongly constrain the possibility of an observable tensor signal being generated. If 
an observable tensor-scalar ratio is considered to be r > 10~^, then only the two or 
three brane cases are capable of producing such a signal. This bound on n depends 
on the WMAP5 limit on f^. If, as expected, the observational limits on tighten 
considerably in the future, the possibility of an observable tensor signal from the multi- 
coincident brane model could be ruled out. 

On the other hand, the choice of r > 10~^ as the threshold of an observable signal 
is very optimistic. If foreground removal techniques and the signal-to-noise ratios of 
future experiments cannot reach this threshold, and instead reach r > 10~^, no number 
of branes will be able to produce an observable tensor signal when combined with 
the current limits on the non-Gaussianity. There will then be little possibility of a 
distinguishing observational signature for these coincident brane models. 

In Part II of this thesis, numerical methods were used to test inflationary models up 
to second order in cosmological perturbation theory. The Klein-Gordon equation at 
second order was derived in Ref. [lii] for the multi- field case. In Chapter 6, second 
order gauge transformations were outlined and the Klein-Gordon equation reproduced 
for a single scalar field model. In contrast to the AM approach, this equation is valid 
on all scales, both inside and outside the horizon. 

In Fourier space, the second order Klein-Gordon equation (6.20) contains a convolu- 
tion term of the first order scalar field perturbation. For this first demonstration of the 
numerical system, a slow roll approximation of the full equation was used. 

Calculating the second order scalar field perturbations provides the possibility of 
a unique insight into the generation and evolution of non-linear contributions to the 
scalar curvature perturbation. One advantage of using the inflaton field equations is 
that we can directly investigate how the perturbations are generated. If, instead, we 
integrated the evolution equation for a derived observable quantity, there would be a 
degree of separation from the physical origins of this process. Indeed, there is, as yet, 
no known evolution equation for the main observable quantity, the comoving curvature 
perturbation, at second order. Using cosmological perturbation theory also provides 
control over the calculation. The domain of applicability of the perturbative expansion 
is well defined and the resultant equations are certain to be valid in this domain. 

The main observable quantity is not however the second order scalar perturbation, 
but rather the departure from Gaussianity in the CMB temperature map, parametrised 
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by the amplitude of the bispectrum of the perturbations. In Section 6.3 we outlined 
how /nl could be calculated from the numerically found 5ip2 both for the local type and 
more generally using the bispectrum of the uniform density curvature perturbation. As 
the observational limits on /nl are tightened over the course of the remaining WMAP 
releases and future Planck data, the importance of comparing the predictions for /nl of 
inflationary models with the observed values will only increase. In this thesis we have 
not computed /nl for the models we have considered, but this is an important future 
step that will be undertaken. 

The long term aim of this continuing project is to analyse multi-field, non-slow roll 
models, in which non-Gaussian effects are expected to play an important role. As a 
step towards this goal, we described the implementation of a numerical calculation of 
the single field, slow roll, second order equation in Chapter 7. The construction and 
evaluation of the convolved source term in Eq. (7.4) proved to be the most numerically 
complex step required. 

To allow a numerical calculation, an energy scale cutoff must be implemented. We 
used a sharp cutoff at small wavenumbers below which the perturbations were taken to 
be identically zero. Another cutoff at small scales (or equivalently large wavenumbers) 
was dictated by practical considerations of calculation size and computation time. 

We defined, in Eq. (7.2), four 9 dependent integrals into which the convolution term 
can be decomposed. By comparison with an analytic solution for a particular smooth 
choice of the first order perturbation, an estimate was made of the relative error present 
in the integration of each term. The number of discrete wavenumber values, their 
spacing, and the number of discrete 6 values were chosen to minimise the error in one 
of the integrals. From these parameter values, three different finite ranges of discrete 
values of the wavenumber were defined. These all contain the WMAP pivot scale at 
^WMAP = 0.002Mpc~^ and cover the WMAP observed scales to varying degrees. Despite 
the k ranges having being chosen to minimise the relative error in the integral of only 
one of the 9 dependent terms, the integrals of the three other terms also display small 
relative errors for these ranges. The analytic solutions which have been found will form 
an important part of the verification of any future modifications to the numerical code. 

The execution of the code is in four stages, building on previous calculations of first 
order perturbations in Refs. [137, 165, 169]. To begin, the background equations are 
solved and the end time of inflation is fixed. The initial conditions for the first order 
perturbation can then be set and solutions found for the evolution equations. Despite 
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the large volume of calculations required at each time step, the easily parallelisable 
nature of the source term calculations allows the run time of the third stage to be 
reduced significantly. The final stage of the calculation uses the source term results to 
solve the second order perturbation equations. 

The initial conditions for the second order perturbations are taken to be 5ip2 = 
and 5(y92 = as described in Section 7.2.2. For this choice of initial conditions the 
homogeneous part of the solution of the second order equation is zero at all times. As 
the perturbations are supposed to become more Gaussian the further back in time they 
are considered, in the limit of the far past the second order perturbations should be 
zero. It remains to be investigated whether the choice of initialisation time is sufficiently 
far in the past for this assumption to be accurate. At first order it is known that the 
perturbations are well approximated by the Bunch-Davies vacuum initial conditions 
even just a few e-foldings before horizon crossing. However, this choice of initialisation 
time may not be the most appropriate for the second order perturbations. In future 
work it would be worth considering whether the analytic Green's function solution for 
5ip2 at very early times could be integrated until the numerical initialisation time and 
used as the initial condition for the perturbation. 

To test the code, four different, large field, monomial potentials were used. These were 
the standard quadratic and quartic potentials, a fractional index potential derived from 
the monodromy string infiation model and a toy model in which infiation is stopped by 
hand and a blue spectrum is produced. Each potential depends on a single parameter, 
which was fixed by comparing the resultant scalar curvature power spectrum with 
the WMAP5 normalisation. The slow roll approximation can be applied to all four 
potentials. These potentials are not meant to represent an exhaustive survey of single 
field slow roll models but are sufficiently different to exhibit different power spectra and 
second order source terms. 

We presented the results for each potential in Chapter 8. The first order results 
match those in Refs. [137, 165, 169]. The results of the source term calculation show 
that before horizon crossing, the source term amplitude decays rapidly for all four 
potentials. The amplitude changes less after horizon crossing, until later times when it 
increases as the slow roll approximation breaks down. 

The four different potentials have similar amplitudes before horizon crossing but 
reach different values after horizon crossing. The differences in the slow roll parameters 
for each potential are compared with the source term values in Appendix A. 6. The slow 
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roll parameters do not appear to be directly related to the amplitudes of the source 
terms, at least in a linear fashion. 

The choice of wavenumber range affects the amplitude of the source term, as expected, 
due to the implementation of a sharp cutoff at large scales. This dependence is only 
apparent, however, before horizon crossing. For a particular range, the magnitude of 
the source term decreases as wavenumber increases. However, the ratio of the source 
term to the other terms in the Klein-Gordon equation increases with wavenumber. 

As expected, the amplitude of the second order scalar perturbations is much smaller 
than that of the first order ones. After the generation of the second order perturbations 
at early times, their evolution is that of a damped harmonic oscillator similar to the 
first order evolution. 

We have shown that the magnitude of the source term can be important throughout 
the full evolution and that it is not sufficient to calculate this term only for modes 
either entirely inside or outside the horizon, i.e., taking a short or long wavelength 
approximation respectively. We have been able to access both of these regimes, by 
solving the evolution equations of the inflaton field perturbation. This is in contrast to 
other approaches which could have been taken, for example using the AJV formalism, 
which is only applicable in the large scale limit. 

The construction of any numerical code involves a considerable commitment of time 
and resources so it is important to understand why such an endeavour has been un- 
dertaken. The numerical calculation of first order cosmological perturbations is an 
invaluable part of the cosmologist's toolkit. It allows analytic predictions of inflation- 
ary models to be confirmed where these exist, but also generates predictions where no 
analytic solution is possible. Another important use is to test predictions based on the 
slow roll approximation against the full evolution equations. We have taken the first 
step towards upgrading this standard numerical calculation to include second order 
scalar perturbations. The ability to check the predictions of inflationary models at sec- 
ond order will be a powerful tool to constrain these models and check the consistency 
of any analytic assumptions that have been made. Solving the inflaton field equations 
provides the most direct access to the non-linear effects that the increase in available 
statistics have made observationally important. 

We have presented the first numerical calculation of the Klein-Gordon equation for 
second order scalar perturbations which was derived in Ref. [133]. Although we have 
restricted ourselves in this thesis to the single field, slow roll version of the second 
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order equation, the expertise gained and the lessons learned in the development of the 
numerical system will be of significant assistance when the next steps towards a full 
multi-field calculation are taken. 

In the past, a numerical calculation on the scale we have achieved would have been 
the preserve of dedicated super-computing facilities. We have demonstrated that a 
calculation of this scope is now possible using relatively modest local resources. If 
the computational power available increases, the practical limits on the resolution and 
extent of the k ranges will ease. Further improvements in the efficiency of the code will 
also loosen these constraints. 

Another consideration in the development of a numerical system is the possibility of 
code re-use. One of our future goals is to develop our code into a numerical toolkit 
which can be applied to a variety of physical situations. The equations of motion of the 
inflaton scalar field are similar in form to the governing equations of other important 
cosmological phenomena. Therefore, it should be possible to adapt the numerical sys- 
tem we have constructed and apply it to other areas of interest. The form of the second 
order equation and source term are similar to those applicable in the evolution of tensor 
perturbations and the generation of vorticity in the early universe. The fiexibility of 
the numerical system we have developed will be a positive factor in any attempt to 
apply our code to these physical systems. 

The numerical calculation described is the first step towards a system capable of 
handling the multi-field, non-slow roll models for which non-linear contributions are 
important. In Section 8.3, the next steps towards this goal were outlined. Continuing 
our work by calculating the second order perturbations for both the non-slow roll, single 
field case and the slow roll, multi-field case will pave the way for the eventual calculation 
of the non-slow roll, multi-field equations. 

The full single field, non-slow roll, second order equation can be treated using the 
method already described. Three more 6 dependent terms are necessary to compute 
the full convolution integral in this case. It will be important to find analytic solutions 
for these three extra terms, as already done for the terms in the slow roll case, in order 
to gauge the effectiveness of the extended code. When the extension to non slow-roll 
models is complete, it will be possible to investigate models with a step or other feature 
in their potential. These models can exhibit large amounts of non-Gaussianity produced 
around the feature with a shape dependence that is more general than that of the local 
and equilateral forms. 
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The Klein-Gordon equation for the multi-field case introduces further complexity. 
We plan to expand the numerical system to encompass two or three scalar fields. The 
differences between single and multi-field models are already apparent for these cases. 
We have presented the slow roll source term equation for multiple fields in vector no- 
tation. The definitions of the four 9 dependent terms used in the single field, slow roll 
model were also extended to the multi-field case. Beyond the slow roll approximation, 
the full multi-field equation should be treatable in a similar manner to the single field 
case, by introducing further 9 dependent terms. 

To conclude, it is worth reiterating our opening remarks. Cosmology has moved from 
being a theorists' playground to a genuine scientific discipline. Infiationary models can 
now be strongly tested by observations and the next generation of experiments will place 
even tighter limits on the viable parameter space of such models. In this thesis, analytic 
arguments have constrained string theory inspired infiationary models and numerical 
methods have paved the way to calculating higher order cosmological perturbations. 



A. Appendix 



The following materials supplement the calculations and discussions in the main thesis. 

A.l. Analytic Solution of Generalised Sound Speed 
Relation 

Eq. (5.3) can be analytically solved in full generality without imposing the limits (5.5) 
on the derivatives of the kinetic function. This allows us to determine the most general 
class of models where the non-linearity parameter satisfies the condition oc 1/Cg at 
leading order. 

In general Eq. (5.3) takes the form 

(2 - a)PxP,xx + 4XP2^^ = '^XP^xPxxx (A.l) 



and this reduces to 



,2 , 3(2 - «) T 



aT,x = (6 - a)T' + ^ ^ , (A.2) 

where T = P^xx/P,x- Eq. (A.2) can be transformed into the linear equation 

C,,, + . ^ ,A.3) 

Za A a 

after the change of variables U = 1/T and the general solution to Eq. (A. 3) is given by 

P,xx 1 



Px X[/2(^)X(-6)/2"-2] ■ 

Integrating a second time implies that 



(A.4) 



',x = -/i(^)(l-/2(v^)X-)'^^'^^ , (A.5) 
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where s = {a — 6) /{2a) and we have redefined the arbitrary integration functions fi{ip)- 
Finally Eq. (A. 5) can be formally integrated in terms of a hypergeometric function 

P = -f,X2FJ--,-^;l--J,X-^) , (A.6) 
\ s 2s s J 

which represents the most general solution for this class of models. Note that we have 
set the remaining constant of integration to zero to ensure that the kinetic function 
vanishes in the limit of zero velocity. In fact this expression admits many different 
classes of solution, arising as limits of the expansion of the hypergeometric function. 



A. 2. Generalised BM bound for Finite n Models 

For completeness we should also consider the BM bound (4.33) for the finite n multi- 
coincident brane models. This is given by 

42 



< -^V1 + (^-1)'>^/n1 , (A.7) 
and in the case of an AdS:^ x X5 throat simplifies to 

-*<-W— (A.8) 

Comparing the limits in Eqs. (5.65) and (A.8) implies that the bound (4.39) is stronger 
than the corresponding BM bound (4.33) if 

n > 1 - 5.5 X lO-i^iV^/WJr/^l , (A.9) 
and this condition is always satisfied if 

- 5.5 X IO-^^AT^/Xs/nI < 1 • (A.IO) 

Moreover, the bound (A. 10) will itself be satisfied for all values of and if it is 
satisfied when the limits = —151 and A^ = 75852 are imposed. Hence, we conclude 
that the bound (4.39) is stronger for Afcs < 75. In general, it is difficult to quantify 
the magnitude of N'cs without imposing further restrictions on the parameters of the 
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models and, in particular, on the functional form of the inflaton potential. However, 
if the ratio Sh/P,x remains approximately constant during the final stages of inflation, 
one would anticipate that ^ 60. Nevertheless, if <^ 75852, the bound (A. 9) will 
only be violated for n < 3 if Ncs ^ 60. 

A. 3. Discussion of Homogeneous Solution for 
Second Order Equation 

The homogeneous equation for the second order perturbations is 

5^p'i{r], k') + 2-HV2(^, k') + [A;2 + aV^^ - 2ATiG{^'Qf] 5ip2{v, k') = . (A.ll) 

During slow roll, with the slow roll variables Eh and r/n defined in Chapter 2, this 
becomes 

V2 + 2^V2 + [k^ + + sh)] Sy^2 = 0. (A.12) 

If we let u = a6ip2, this equation can be rewritten as 

u" + [k"^ + n^{3r]H - 2eH - 2)]u = 0. (A. 13) 

When Eh is small, the conformal time rj is given by 

so we can rewrite Eq. (A. 13) as 



u" + 



2 , 1 3r]H - 2eH - 2 



k' + 



M = 0. (A.15) 



If the derivatives are taken in terms of {—rj) instead of rj this is in the form of a Bessel 
equation with solutions in terms of Hankel functions given by 

r(l,2) 



u^,2 = ^/^Hl'''\-kri), (A. 16) 



A. 3: Discussion of Homogeneous Solution for Second Order Equation 



174 



where ' are the Hankel functions (Bessel functions of the third kind), and z/ is given 
by 

.2 = ^'^^ + ^ ■ (A.17) 

4(1 - 2eH) ^ ' 

The full solution for u is then 

Mfuu = Ci^H';^\-kri) + C2V^H^J\-kTi) , (A.18) 

where Ci, C2 G C. When the (real) argument of the Hankel functions goes to +00 they 
have the following asymptotic form [ ]: 

i/(i)(^)^./Xe^(-f-f), (A.19) 

V TIZ 

/Ae-(-f-f). (A.20) 

V 7CZ 

So at early times when rj — > ~oo and —kr] — )■ +cxd we have the following expressions 
for u: 

= .[X (c,e-'ikv+^^-^+f) + c^e+'ikri+^^'^+m ^ (A.21) 
\ nk 



k.[T (_C7ie-*('=^+t'^+i) + C2e+*('^+t'^+^)) , (A.22) 
\ 7rk 

(A.23) 

where we have assumed that u is slowly varying far in the past, i.e., the derivatives of 
the slow roll parameters are very small. 

As explained in Section 7.2.2, the results given for 6(p2 are for the full solution in- 
cluding the homogeneous part. To remove the homogeneous part of the solution the 
initial conditions for the full 6(p2 should be chosen such that Ci = C2 = at all times. 



A. 4: Analytic Tests for B, C and V Terms 



175 



A.4. Analytic Tests for B, C and V Terms 

Analytic solutions can also be found for the B, C and T> terms. The B term integral, 
Ib, is given by 



Ib = 2tt dqq^5^,{q')B{k\q') 

^ ^min 

^27ra^ dqq^ dO {k^ + q'^ - 2kqcose)-^^* cosO sinO , (A.24) 



and has the following analytic solution when 5(pi{q) = a/y^: 



168A;2 



I 



63k^ 



log 



Vk 



y/k + /Cinin 



+ log 



— + arctan , 
2 Vy^^ 



(-65A;^ + 8kkl^) (yk + /c^ax + VA;^ax - A;) 



(65A;3 - S/c/c^iJ (VfcTfc" - ^/k - k^^ 



+ (-22A;2/,^i^ + 16A;^. J + ^min + \/k - k^^ | . (A.25) 



If, in addition to 5ipi{q) — a/^, we also take 

5^\{q) = 



a .dy/q 



(A.26) 
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then the C and T> terms can be integrated analytically. The integral of the C term is 



%= / d^q6ip,{q^)6ipl{k'-q') = 27r / dqq'6ip,{q%{k\q') 



= -27rQ;^ / dggz / [[k"^ + q^ - 2kqcos\ 

J femin J 



+ ^ {k'^ + q^ -2kq cose) ^ ) sinM^, 



(A.27) 



and the analytic solution is given by 



U - 1 



240pk 



15k^ 



log 



\/k + /Cmin n 

+ log 



— + arctan ( "'""'^ 
2 Vv^^^ 



Vk + k 

max 



(15A;3 + 136A;A;^3,) ( V^T^ + V^ma. - k^ 



+ {llSk^k^^ - 48A;^^) + k^^ - y/k^^ - k^ 
y/k~ (15/c^ + I36kkl^) (^^/kTk~ + \/k ^min^ 



(A.28) 



The integral of the V term is 



I^^2n dqq'Sip,{q')V{k\q') 



(A.29) 



—27ra^ I dqq^- 

^min *^ 



k'^ + q^ - 2kqcose) 



+ + (f - 2fcg cos e) * ) cos ^ sin OdO , 



(A.30) 
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and the analytic solution is 



-In - i- 



900/3k^ 



135^ 



log 



max k + y/k~\ I y/k + k 

max 

+ Vk 

n 

+ log 



TT 

— + arctan . 

2 V y/k — k„ 



Vk 



Vk + k 

mm n 



{-l%bk^ + 168^2^^^ - 32fc^,J {yk + k^^ - x/a"^ 



+ (70A; /Cmax + IS/^/Cjjjg^^) {^\/~k^r~kmax ~\~ V^^max ^ 



(_185A;^ + IQ^k^kl^ - 32A;^,J {^k + k^m - ^/k-k. 



+ (70fc3/ci„in + IQkkl^^) ( V/c + fc^in + 



(A.31) 



A. 5. Analytic Solution for Source Term 

Suppose the first order perturbations are given by the non-interacting de Sitter space 
solution such that 

5^i(r^,F) = ^fl-f , (A.32) 



and the derivative in terms of is 



2k \ kf] 



:V^. (A.33) 



2A; V kr]J \ aHr] J 

The analytic solution of Eq. (7.14) for this choice of first order solution can be written 
in terms of four integrals of the A-T) terms: 



{ J4 + Jb + Jg- + } , 



(A.34) 
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where 



^ (aH) 



2„2 



2k^q 



min 



min 



{any 



9 _ g2 
2 A;2 



kq)6^W)V{k\q'). 



5^,{q')A{k\q^), 
(A.35) 
(A.36) 

(A.37) 

(A.38) 



The analytic solution for J4 is given by 



J A = 



V 



( 



240A; arctan 



k — k 



{rfk"^ - m-qk - 24) 



- 120A;7r (77^^;^ - m-qk - 24) - 80v^ 



- Uk 



kmax k -\- y/k + /Cmax^ ^ + Sfcjj^g^-jj- ^\//jmax ^ ^/k + /^max^ 
+ 48i ^A^max + /Cmax 'x/^max + ^ ^V'^max ^ + + /Smax^ ^ f] 

+ 72 ^^//c + /Cmax A/^max ^ 



-80v^ 



3 ( yk A^min + 
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+ 240A; {rfk'' + 24) log (2\fk) - 240k {rfk^ + 24) log (2 (^A" + ^A"^)) 
- 240A; {rfe + 24) log (2 [y/k^ + ) 

+ 240A; {rfk'^ + 24) log (^2 a/A; 

min 

+ 3v^ (31777^A;^ - imOi-qk - 560) A;^ 



+ 45 {ri^k'^ - Uirik - 16) arctan | ^/ ) A;2 + 45 (r^^/c^ + 8ir/fc + 16) log hy/k) k^ 



- 45 + 8ir//c + 16) log (2 + \/Kr..-k)) 

- 45 (ry^/c^ _ siry/c + 16) log (2 + + /c^ax)) A;^ 

+ 45 (77^A;2 - SiTyA; + 16) log (2 (^^/k~ + V^- + Aw)) - y (^^^^ - 12ir;A; - 16) ttA;^ 

max k y/k + A^max) A^ 

+ 1077(?7A;^ax + 12i) (v^ max k + •y^A; + A^max) A^ 
+ 8 (r/^fc^ax + lO^^^/Cmax + 30) (^\/A;max - k - a/T+I") A;^ 
- 16A;inax (llr7^A;^ax - 20i7?A;i„ax - lO) (^V^max - A; + V^aT+A") A; 

+ 128A;jjj3j(. (fj kj^g^ — birjkmecK ~ 5) ^V^A^max ~ k — \J k + A;niax) ^ 

- '^\J k^J\lhrf' i^k- A^min + a/A; + A;inin) A;^ 
+ 10?7 (^?7A;i„in (^V^ - A;min - + ^min) " 6i (^3 V A: - A;min + 2 VA: + A;i„in) ) k^ 

+ 8 {^\Jk- A;min + + ^min) /^min " 5^?? (^\J k - Al^in " 2 VATTaJ") A;min 

+ 30 (^a/A^TAw - - ^min) ^ A;^ 
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16/Crni 



^11?7 (^y/k A^uiin y/k + /Cmin^ ^min lOiTj (^y/k /Cmin 2-\//c + /Cmin^ ^min 

+ 10 

+ 64A;jjjj^ ^2?7 ^ ^/k — /Cmin + \/A; + /Cmin^ ^min ^ ^//c — /Cmin — 2 + k^^^i^ ^min 

+ 10 ( - ^min - + ^min) ^ ^ | • (A.39) 



The analytic solution for is 



- -^M L.oo'^L2u I 29400 (4^^fc^ + I5ivk + 120) arctan ^ ' ^ 



aH)^^" 2 2822A00r]^k\ ^ ' ' > \\ k - k^^u 

+ 29400 {4:ri^k^ - hli-qk - 120) log {2VPj k^ 

- 29400 {4r]^k^ - 51tr]k - 120) log (2 (v^A" + ^A"^)) fc^ 

- 29400 {4:ri^k^ + 51ir]k - 120) log (2 (^y/k^ + a/^T^)) k^ 
+ 29400 {47]^e + 51i7]k - 120) log (2 (^y^k~+ ^/k + k~}j^ 

- 14700 {Arfk'^ + ISiTyA; + 120) T^k^ + 280VA;inax(A; + k^^) ^4207^^^;^ 
+ 577(20077^:^3^ + 303i)A;2 + (-416772^:^3^ + 6414i77A;^a^ + 11592) ^ 
2Ak^^ {STj^kl^, - mrikr,^ - 126) k + 48A;Lx (Sr^'^^x + 53^7;A;„,a^ + 84) 



- 280^/A~(A~^|^420772A;^ - 577(20077A;„,a^ + 303i)A;3 
+ (-416772A;Lx + 6414i7;A;„,a. + 11592) ^ + 24A;n,ax {^rj'kl^ - Slir^k^^ - 126) k 
+ 48A;^^ (877'A;^3, + 53i77A;^a. + 84) 
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- 2S0^/k~^/k + k~ymr^^k^ + 5rj{200rik^in + S03i)k^ 

+ {-AWrj^kl^ + 6AUirik^,n + 11592) k^ - 2Ak^i^ (Sry^A:^.^ - 87irik^^ - 126) k 
+ {Sv'kl,^ + 53i7^k^^ + 84) 

- 280^/{k^^k~n)k~(^20r)'^k'^ + 57^(1083^ - 200r)k^in)k^ 
- 2 (2087^2^^.^ + 2547ir7A;n,in + 5796) + 24A;„iin {Sv^k^^^ + 67irik^in + 126) k 

+ ASkli, {Sv^kl^ - 7Sir]k^i^ - 84) 



^'^ 105 {27{)rfk'^ + 2765ir;A; + 4956) arctan ^ ■ / ^""'^ 



(0/^)2 28224007^2^3 1 V / ' / ' j \\k-k 

+ 105 (27077^A;^ - 37i5irjk - 4956) log (2Vk^ k^ 

- 105 (27077^/^2 _ 3745i77A; - 4956) log (2 (^>/k~ + /c 

- 105 {270ri\'' + 37A5irik - 4956) log ^2 T^T^)) k 

+ 105 {27Qrfk^ + 3745i?7A; - 4956) log (^2 (^^/k~ + ^A^fAw)) A;^ 
105 

- — (27077^ A;^ + 2765i77A; + 4956) ttA;^ 



- Vk^^^K^^^^^(^S350ri^k^ + 525ri{36rik^^ - 749i)A;^ 

+ 70 (2 1 67^2^^^ _ 3745,^A;,nax - 7434) k' - 40A;„ax {3516ri^kl^^ - 133irik^,^ + 8673) k' 
- A8ki,, (13607/2A:Lx - 18655z77A;^a. - 22442) ^ 
+ 64A;^ax (60077^ A;^3^ - 5110i77A;inax - 5733) k 

+ 2hQkt,^ (3007;2fc^^ + 1855z7;A;^ax + 2646) 
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+ y/kraUk + kms^) 1^2835077' A;^ - 525r){36r)kr,^ - lAQi)^ 

+ 70 (2167^2^^^ _ 3745,^^^^ _ 7434) k^ + AQk^^ {3516rj^kl^ - ISSirik^^ + 8673) k^ 
-A8kl^ {mOrj'kl^^ - 18655ir)k^^ - 22442) e-QAk^^, {(^OOrj'kl^ - BUOirjk^^ - 5733) k 

] 



+ 256k^^, (SOOt^^A;^^ + 1855i7]k^^ + 2646) 



- y/k~^k + krain ^2835077^ A;^ - 52577(3677A;min - 749i)A; 

+ 70 {216ri^k'^^ - 37A5irikrnin - 7434) k^ + AOkram (SSier^^A;^;,, - ISSzT^/e^in + 8673) k^ 
-48A;^in {UGOri'^kl^^ - 18655irikrain - 22442) k'^-Gik^^^ {GOOrj'^kl^^ - dllQirjkram - 5733) k 

+ 256A;^i, (300r/'A;^i, + 1855ir;A;„,in + 2646) 

- V{k - ^min)A;min (^8350r}'^k^ + 525r]{36r}kr^in + 553i)k^ 

+ 70 (216?7^/ef„i^, + 2765irik^,in + 7434) A;^ - AOk^-.^ (3516??^A4i„ - 9247ir7A;n,i„ - 8673) A;^ 
-48A;^i, (l360r;2A;^i„ + 15155z77A;^in + 22442) A;2+64A;^.^ (eOOr;^^^.^ ^ 3710^r7A;^i„ + 5733) k 

+ 256A;^i„ (300r7'A;^i„ - 2555irik^i^ - 2646) j | . (A.40) 

The analytic solution for is 



5 



+960ir)k-80p^r)'^ {rfk"^ - 12ir]k - 24) +20^77 {-3irfk^ - 32rfk'^ + 96i77A; + 192) +1920 
15 / 

+ — A;7r 9ri^k^ - mirfk^ - BQOrj'^k'^ + 960i77A; - 80/3^7]^ {ri^k"^ - Uirjk - 24) 



2 



+ 20/37; {-3irfk^ - 327fk'^ + 96i77A; + 192) + 1920 
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+ 38A0i{k - k^^fk^^T}^ + 80 {69k^ - 178k^^k + ISAkl^) rf 
+ 400/3^ ((3A;2 - Uk^^k + 8Cax) + m{k - k^^)rj - 72) 77^ + 19200i(A; - k^^)ri 



+ 320/3 



12z(A; - k^^fk^^rf + {21^ - Q2k^^k + 56A;^^) if 



+ 120i(A; - k^^)r] - 180 



77 - 28800 



+ \/A;^ax(A; + A;„a.) (l5A;^ - 10A;„,a.A;=^ - 248A;^3,A;' - 336A;^3,A; - 128A;iJ 77^ 

+ 3840iA;^ax(A; + k^^^frf + 80 (69A;' + llSk^^k + 184A;^3,) 77^ 
+ 400^' ((3A;' + 14A;^axA: + 9>kl^) rf - A%i{k + A;^a.)77 - 72) 77' 



19200i(A; + A;„iax)77 + 320/3 



12iA;^a.(A; + /^ma^j'ry^ + {2lk^ + 62A;^a.A; + 56A;^^) rf 



l2Qi{k + A;^a.)77 - 180 



77 - 28800 



+ ^/k~^/k + k~[ 9 {-Ihk" + 10A;„,inA;=^ + 248A;^i„A;2 + 336A;^,A; + 128A;^iJ 77^ 



- 3840iA;„,in(A; + k^i^)W - 80 {mk'' + n9>k^ir.k + 184Ci„) rj' 

- 400/3' ((3A;' + 14A;„,inA; + %kl,^) rf - m{k + k^^)r] - 72) 77 



+ 19200i(A; + krain)r] + 320/3 



- I2ik^^{k + k^i^frf - {2ie + 62k^ink + 56kl,^) rf 



+ 120i(A; + A;inin)77 + 180 



77 + 28800 j 
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+ 60z (15A;3 - 54A;^in^2 ^ g^2^,^^ ^ ^^g^s^, j ^3 _ gQ ^39^2 _ 3g^^.^^ ^ 104^;^^^ ^^ 
+ 400/3^ ((3A;2 - 14A;^i„A; + 8A;^iJ 77^ + 12i(A; + 2k^,^)r) + 72) 77^ 



+ 4800i(A; + 2k^■,^)^ + 20i/3 



3 (15A;3 - 54A;„,i„A;' + 8A;^i,A; + IGA;^;,) 77^ 



?7 + 28800 I 



+ 32i (3A;2 + 4A;„ii„A; + 8^;^;^) ^^ + 480(A; + 2A;inin)?7 - 2880i 
+ \hk {9r]^k^ - mri'^k'^ - 320/37? {rj'^k'^ - 12) + SOP^r]'^ {rj^k'^ + 24) + 1920) log 
- Ihk I Qt/^/c^ - 40077^^2 - 320/37/ (77^/c^ - 12) 



+ 80/3277' (77'^;' + 24) + 1920 j log (2 



- 15A; |^977^A;^ - 40077^^;' - 320^77 {ri^k'^ - 12) 
+ 80^^77' (77' A;' + 24) + 1920 j log (2 + v^^+^)) 



+ Ibk \^r]^k^ - Amrfk'^ - 320/3?/ [rfk"^ - 12) 
+ 80/3V(^'A;2 + 24)+ 1920 jlog(2(y^+^/^Tfc~)) I. (A.41) 



The analytic solution for J~ is 



+ 20/3^77' (277'A;' + 13i7/A; - 12) + 10/377 (2777'A;' + ^2i^k - 48) - 240 
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+ 2520A;7r ^^ri'^k^ + 230r]^k^ + 260ir)k + 20p^r)^ {2ri^k^ + ISirjk - 12) 



+ 10/3r/ (27r/2fc2 + 52ir}k - 48) - 240 ) + ^^^^+^^-^-- ( _ 35^2 | ^^^^2^, 

kri 



2 7.2 
max 



+ 192irik^^ + 6A/37]{3i7]k^,, + 1) + 64 j A;" - 5r; |^ - 2M7i-'ki^ - llli^'ki,, + 2656r;A;„,a. 
+ 9mi3''rf{2,r]k^^ - i) + {-^^^iv^k^^ + 5536r)k^^ - 1920i) - 960ij k^ 
- 6 1^ - 58877^ A;^^ + ISOBirj^ki^ + 9965r)^kl^ + 15520ir)k^^ 
+20PY {^Wkl^. + 77Givkms. + 252) +5/377 ( 261i77^A;^3,+2173772A;^3,+6208i77A;„,a.+2016 



+5040 |/c2+4/c^,.| 84r^^/cl,-330^77'/c^^-3275r;2fc^^+7065^77/Cmax- 15/3^77^1 20r]'ki^ 



- mir]k^^^ + 1008 j - 5/37/ {miri^kl^ + 7157/^^^^ - 2W7)k^^ + 6048) - 15120 ) k 
+ ^kl^ - 84r/'Cax + 330z773/c^^ - 145077^^;^^ + 2385z77A;^a. + l^^W {^^kl 
+ l^9ir]k^^ + 378 I + 10^77 (33^77=^ A;^^, - 11577^^;^^ + Allir^k^^ + 1134) + 5670 
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+ 16A;-^V^ma. - k 3577" (lll?7'Cax + 192i77A;^ax + 64.pr){3ir)k^^ + 1) + 64) 



+ 5r/ 1^29 V/c;^^ + TlTir^^/c^,, - 2656r;fc^a. - ^m'n\^nkr.^ - i) 



3u3 
max 



+ il3ri {7nrfkl^^^ + bbmr]k^^ + 1920) + 960i j A;^ + 6 1^ - SSSr^^/c^^^ + ISOSir^^fc 

+ m^fyrfkl^ + 15520ir7A;^ax + 20^^77^ (45772^,-^^^ + imr^k^^ + 252) 
+ 5^77 (261i773A;^^ + 2173772^;^^ + 6208i77A;^ax + 2016) + 5040 j ^ + Ak^^ i^W^t 
- ?>mrfkl^^ - 3275ri^kl^ + imbir]k^^ - Ib^^rf {2Qrfkl^ - mir]k^^ + 1008) 



5/377 (66^77^^;^^^ + 7157^2 A;^^^ - 2826i77A;max + 6048) - 1512oj k 

\^ - m'kL. + ^^Ozfj'ki^ - UbOrj^ki^ + 2385z7;A;^ax + 15^^ (^Orj'k, 

+ 159i77A;^ax + 378^ + IO/377 (33^772^:^3^ - 11577^^;^^ + 477ir}k,^^ + 1134) + 5670 

+ 16A;-f/V^ + ^min 1^3577' (lll?7'A;^in + 192i77A;^i„ + 64^77(3i77A;^in + 1) + 64) k^ 

- 577 ^294772A;^„ + mirfkl,^ - 265677A;^in - 960^'77'(377A;^in - i) 

+ il3r] (7177/2 + 5536i7/A;min + 1920) + 960i j /c^ + 6 ^ - 5887;^A;^i„ 
+ 1305^77=^ A;^„ + 996577^^;^;, + 15520^77^;^^ + 20/3^77^ (45772A;^i, + 776i77A;^i„ + 252) 
+ 5^77 (261i772A;^i„ + 2 173772 A;^i^ + 6208^77 A;^in + 2016) + 5040 ) k'' 
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- A71ir]k^in + 1008^ - b^r] {66ir]^kl^^ + 715r?^A;^i„ - 2826ir)k^^ + 6048) - 15120 j k 

+ 159ir]kr,i^ + 378^ + lO^r? (aair^^A;^^ - IISt^^A;^^^ + 477ir)k^i^ + 1134) + 5670 
- 16A;-f/V^-^min 1^3577' {niv^k^^ + mirjk^^ + 6A^r){3ir]k^^ + 1) + 64) k' 
+ 1077 1 14777=^A;^i„ - 1104i772Cin + 155277^;^!^ + 480^2772(377A;^i„ - i) 



3i,3 
min 



+ 16/377 (-69^772^;^;^ + 18777/c„iin - 60i) - 480i J A;^ + 6 1^ - dSSrj'^k^^ + A80ir]^k, 
+ 8165772Cin + 14720i77A;^i„ - 20/3^ (4577'Cin - 736^77^;^;^ - 252) 
+ 5/377 (96ir7^A;^in + 1453772A;^i„ + 5888i77A;i„in + 2016) + 5040 j k"^ 

+ 4A;min 1^8477^ + 120i77='A;^i„ - 2675r]^kl^ + 6165ir]k^n + (^^Orj^kl^^^ 

+ Allirjk^i^ - 1008^ + BiPrj (2477^ A;^;^ + 475ir]'^kl^^ + 2A66r)k,^,, + 6048i) - 15120 j k 
+ 8Ci„ (sAr)^ki^ + 120i77^A;^i„ + 205077'A;^i, - 3285i77A;^in + 15^'77' (^Ori'^kl^ 
- 219ir)kmin - 378 ) + lOprj {l2i7]^ki^ + 2Zhrfkl^^^ - 657ir]k^in - 1134) - 5670 



A. 5: Analytic Solution for Source Term 



188 



- 5040A; ^ - drj'^k^ - Ahirfk^ - Ib^rj^k'^ - 3A0ir)k + 20^^r)'^ {2r)'^k'^ - 17ir]k + 12) 

+ 5/377 {-^irfk^ - 22ri^k'^ - 136ir)k + 96) + 240^ log (2Vk^ 
+ 5040k - 9r)^k^ - Abirf^k^ - IbOrj^k'^ - SiOirjk + 20P^r]^ {2r{^k'^ - nirjk + 12) 
+ 5/377 {-Qiv^k^ - 22rfk'^ - ISGirjk + 96) + 240^ log (2 + V^Wx- A;)) 

+ 5040A; - Qrj'^k^ + Abirfk^ - IhOrfk'^ + MOir]k + 20^^77^ {2ri^k'^ + lli-qk + 12) 
+ 5^77 (9i77^A;^ - 22772^;^ + 136i77A; + 96) + 240^ log (2 [^/k^ + ^/k + K^^ 

- 5040A; - 9r}^k^ + Abirfk^ - IbOrfk'^ + 340i77A; + 20/3^77^ (2772^;^ + 17^77^; + 12) 

+ 5/377 (9i77^^^ - 2277^^;^ + 136^77^; + 96) + 240 j log (2 (^^/k~ + V^ + /Cmin)) 
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+ 20/3^7^2 {iri^k'^ + Ibi-qk + 120) + IQ^-q {Qirfk^ - hrfk'^ + mi-qk + 480) 
+ 2400 j arctan ^ j A;^ + 315 ^9r;^A;^ + lOir^^A;^ + 2Q{)rfk'^ - 1020ir;A; 

+ 20/3^7^2 (4r;2/c2 - Sliryfc - 120) + lO^ry {irfk^ + 2>7ifk'^ - 2mr]k - 480) 



- 2400 j log (2\fPj k^ - 315 {^-q'^k^ + mrfk^ + 29077^^:2 - 1020^77^; 
+ 20/3^77^ (477^/^2 - 51i77A; - 120) + 10/37? (i7;^A;^ + 2,7rfk^ - 2QAir]k - 480) 



- 2400 j log [2 (^v^A~ + V/c„,ax- fcjj A;'^ - 315 |^97;^A;^ - 10i77=^A;^ 

+ 2907/=^ + 102077/A; + 20/3^7?^ {Arfk'^ + 5177/A; - 120) + 10/37/^ - 777^ A;^ + ?,7r{^k'^ 

+ 204i77A; - 480^ - 2400 j log (2 + ^/kTk~^^ k^ + 315 ^977^ A;^ - 10777^ A;^ 

+ 29077^A;2 + 1020777A; + 20^^77^ (4772A;^ + 51777A; - 120) 
+ 10^77 (-777^ A;^ + 37ri'^k'^ + 2mr]k - 480) - 2400 j log (2 (yk~ + a/A^TAw)) A;^ 

- 977^ A;^ + 60777^ A;^ - 13077^A;^ + 300777A; + 20^^77^ (477^A;^ + 15777A; + 120) 

+ 10/377 (6777^ A;^ - 577^ A;^ + 60777A; + 480) + 2400 j ttA;^ 



315 



+ \/K^.Vk + k^A 28357j^k^ - Qmri\bl3rj - Sik^^^rj + 5)k^ + Utj^ I8OO/3-7; 



2„2 



- U28ki^r]^ + 27mkme^r) + 5/3{542ir]k^^^ + 3201)7? + 14205J A;' 
+ 477 ( 236477=^ A;^3^ + 6620777' A;^^^ - 946577A;^a^ + 75/3^77' (20077A;^ax + 303?) 
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+ 5/377 (l324i77^A;^^ + n07r]k^^ + 9090i) + 22725i j 
- 24 1 - 139277^fc^^ + 2500^7^3 fc^^ + ISOlOr^^A;^^^ - 16035ir]k^^ + b^'^r)'^ { 208ri'^k 



-3207ir]k^^^-5796j +10/3r/ (250^772^^3^ + UObri^kl,^ - 3207irik^^^ - 5796) -28980 j A;^ 

+ 160/c^ax (^iv'kt.^ - SSirj'ki^ - 597r)'kl^ + 783irjk^^ - 9(3'rj' (^8rj'ki^^ 
- S7ir)k^^ - 126^ + ^77 (-88i77=^A;^,, - 66977^^;^^ + 1566i77A;„,a. + 2268) + 1134^ k 
+ 320Ca. - '^^V'kt^ + SSirj'ki^ - SA8r,'kl,^ + A77ir)k^^ 
+ 9PY {Sv'ki^ + ^Sivkms. + 84) + 2^77 (^Uirj'kl^ - ^^Wkl 

— •)"» 

+ ^A~v^A~^<J - 283577^^^ - 63077^(5^/377 + 3A;„,a.7; + bi)k^ - Urj" [ 18^13^ 



,2 

''max 



- 1428A;^^7;2 + 2710iA;„,axr7 + 5/3 (542^7; A;n,ax + 3201)7; + 14205 j k'' 

+ 477 |^2364772A;^3^ + 6620i772A;^3^ - 946577A;i„ax + 75/3^77^(20077 A;max + 303i) 

+5^77 (1324^77^ A;^3^ + 110777A;max + 9090i)+22725i j k^+2A {^I392r]^ki^^+2bmirfkl^^ 
+ 13010772A;^^ - 16035i77A;„,ax + 5^277' (208772A;^^ - 3207i77A;^ax - 5796) 
+ 10/37/ (250i7/2A;j^^ + 14057/^ A;^^^ - 3207i7/A;max - 5796) - 28980 j 

+ leO/c^ax 1^2477^/0^^ - 88z7;='A;^,, - ^97rfkl^ + 783z7;A;^ax - '^I^W {^rfkl,^ - 87ir]k^.. 

- 126 ) + ^77 (-88^77^^;^^ - 66977^^:^3^ + 1566i77A;^a^ + 2268) + 1134 ) A; 
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+ 1428A;^i„772 - 1660ik^,,rj + 5/3{-332ir)k^^ - 1761)7? - 10605 j k'^-Arjy- 2364r]^kl^^ 
+ 13A80ir]''kli^+3M65r]k^i^+7bpY{'2QQvkn,in-l083i)+5P7](2696ir]^kl^^^^ 



- 32490i^ - 81225i j k^ - 24 ^392rj'^k^^ - lOOOirj^k'^^ - 10930r]'^ k^^^^ + 12735ir]k^ 



+ 5/3Y (208772A;^i„ + 2547ir)kmin + 5796) + 10/3r) - lOOirj'^kii^ 
- 989ri^kl^ + 2547i7yA;^in + 5796^ + 28980 j k^ + mk^i^ (^4:ri^kiin + 32ir^^A;^i 
-453r)^kii^ + 603ir}k,^n + 9P^r]^ (8r7^A;^in + Q7ivk^m + 126) +Pr}(32ir}^kli^-381r)^k 



2 

min 



+ 1206ir^A;inin + 2268^ +1134j k + 320ki^ ( 24r)^k^^ + 32ir)^kii^ + 492r)^kii^-657ir]k,, 



+ 9PW {^rfkl^ - 73irik^in - 84) + 2Pri\^16iri^ki, 

+ Vk~Vk + k~(^ - 2835r)^k^ + 630r)%5i^r] + 3k^^r] + hi)k^ - Urf (^8mpW 
- 1428A;^i„r?2 + 27mk^s.^r] + hp{h42ir]k^i^ + 3201)77 + 14205^ k'^ - Ar) (236Ar)^k^ 



min 



+ 6620^77^^^.^ - 946577A;„,in + 75/3^772 (20077A;^i„ + 303i) + 5^77 ( 1324i77'A;^i„ + 110777A;„ 
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+ 9090zJ + 22725iJ /c^ + 24 | - 1392r]^k^-^^ + 2500i7]^ki^^ + ISOlOr^^fc^^i^ - 16035ir]k^ 



+ 5/3^7]^ (20877^^;^.^ - 3207i7]k^i^ - 5796) + 10l3r]\250ir]^ki^^ + 14057/^^;^;^ 
-3207ir/fcj^in-5796 ) -28980 | k'^-lGOk^i^ I 24ri^k^i^-88iri^kl-^-597r]'^kl^^+783ir]kr, 



„2 

''min 



- 9PY {8v^ki,^ - 87tvk^,^ - 126) + f3r]{^- 88ir]'ki,^ - m9rj^kl 
+ lbmir]k^^+22Q8\+ll3A\k-32Qkl^(-2Ar]^k^^,^+88i^^ 



+ 9PW {8rfkl,^ + 53i77A;^in + 84) 
+ 2137] [AAirfkl^^ - 1381^^^,^ + A77ir]k^^ + 756) + 756^ j | . (A.42) 

A. 6. Discussion of properties of source term for 
different potentials 

The evolution of the source term for the four potentials has been discussed in Sec- 
tion 8.2.3, with particular emphasis on the evolution after horizon crossing as shown 
in Figure 8.15. Here the differences apparent at early times, shown in Figure 8.16 are 
commented on. 

At early times the first order perturbations are still very close to the Bunch-Davies 
initial conditions as outlined in Section 7.2.2. In particular the perturbations are highly 
oscillatory with phase exp(— A;?]), where rj is the conformal time. When eh is small this 
is given by 

It is therefore instructive to plot the slow roll parameter eh for the four potentials at 
these early times, as has been done in Figures A.l and A. 2. For completeness the other 
slow roll parameter rjn defined in Eq. (2.27) has been plotted in Figures A. 3 and A. 4. 

Figures A.l and A. 3 show Eh and rju for the four different models. Figures A. 2 and 
A. 4 show the early stages of the evolution as in Fig 8.16. 

As is clear from these figures the change in the slow roll parameters is not easily 
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Figure A.I.: The value of Eh for the four potentials. 
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Figure A. 2.: The value of Eh for the four potentials at early times. 
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Figure A. 3.: The value of 7]h for the four potentials. 

related to the differences in the profiles of the four potentials in Figure 8.16. In partic- 
ular, although Eh and rjn are quite different for the quadratic and quartic models, the 
magnitude of S after horizon crossing for these models is very similar. 

At the earliest stages of the calculation of S, one or two e-foldings after the ini- 
tialisation of the first order perturbation, there appear to be small oscillations which 
affect the models in different ways. The highly oscillatory initial conditions, combined 
with the small but appreciable differences in Eh and tih contribute to this effect. In 
Figure A. 5 the real part of the phase of the initial condition for 6(pi is plotted just 
after initialisation for the four potentials. The small differences in phase for each model 
combined with the sharp cutoff at large and small k values could explain the variations 
in 15*1 at early times as seen in Figure 8.16. 




V{lp) = |rnV' 
V{ip) = |A(^4 
V{ip) = aipi 
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Figure A. 4.: The value of tih for the four potentials at early times. 




Figure A. 5.: The real part of the phase in the Bunch Davies initial conditions for the 
four different potentials at early times. 
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